计蒜客NOIP2017提高组模拟赛（三）day2-直线的交点

传送门

简单几何+逆序对

发现当两条直线甲乙与平板的交点在上面甲在较左的位置，那么下面甲在较右的位置就可以相交

然后把上面的位置排下序，下面离散化+树状数组即可

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<queue>
#include<vector>
#define INF 0x7f7f7f7f
#define pii pair<int,int>
#define pdd pair<double,double>
#define ll long long
#define MAXN 
using namespace std;
int n;
int K,A,B;
namespace solve1
{
    int read(){
        int x=0,f=1;char ch=getchar();
        while(ch<'0'||ch>'9'){if('-'==ch)f=-1;ch=getchar();}
        while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
        return x*f;
    }
    int dat[100005];
    int n,N;
    int sum(int k){
        int ret=0;
        while(k<=N){
            ret+=dat[k];
            k+=(k&-k);
        }
        return ret;
    }
    void add(int k){
        while(k>=1){
            dat[k]+=1;
            k-=(k&-k);
        }
    }
    pdd Node(int k1,int a1,int k2,int a2){
        double x=1.0*(a2-a1)/(k1-k2);
        double y=1.0*(k1*a2-k2*a1)/(k1-k2);
        return make_pair(x,y);
    }
    int K,A,B;    
    int k[100005],a[100005];
    pdd b[100005];
    pair<double,int> c[100005];
    int d[100005];
    void solve(){
        K=::K,A=::A,B=::B;
        n=::n;
        for(int i=1;i<=n;i++){
            k[i]=read();a[i]=read();
        }
        for(int i=1;i<=n;i++){
            b[i]=make_pair(Node(k[i],a[i],K,A).first,Node(k[i],a[i],K,B).first);
        }
        sort(b+1,b+n+1);
        for(int i=1;i<=n;i++){
            c[i]=make_pair(b[i].second,i);
        }
        sort(c+1,c+n+1);
        N=0;
        for(int i=1;i<=n;i++){
            if(1==i||c[i].first!=c[i-1].first){
                N++;
            }
            d[c[i].second]=N;
        }
        ll ans=0;
        for(int i=1;i<=n;i++){
            ans+=sum(d[i]);
            add(d[i]);
        }
        printf("%lld
",ans);
    }

}
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if('-'==ch)f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

int main()
{
//    freopen("data.in","r",stdin);
    K=read();A=read();B=read();
    n=read();
    solve1::solve();
    return 0;
}