计蒜客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=,f=;char ch=getchar();
         while(ch<''||ch>''){if('-'==ch)f=-;ch=getchar();}
         while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
         return x*f;
     }
     int dat[];
     int n,N;
     int sum(int k){
         int ret=;
         while(k<=N){
             ret+=dat[k];
             k+=(k&-k);
         }
         return ret;
     }
     void add(int k){
         while(k>=){
             dat[k]+=;
             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[],a[];
     pdd b[];
     pair<double,int> c[];
     int d[];
     void solve(){
         K=::K,A=::A,B=::B;
         n=::n;
         for(int i=;i<=n;i++){
             k[i]=read();a[i]=read();
         }
         for(int i=;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+,b+n+);
         for(int i=;i<=n;i++){
             c[i]=make_pair(b[i].second,i);
         }
         sort(c+,c+n+);
         N=;
         for(int i=;i<=n;i++){
             if(==i||c[i].first!=c[i-].first){
                 N++;
             }
             d[c[i].second]=N;
         }
         ll ans=;
         for(int i=;i<=n;i++){
             ans+=sum(d[i]);
             add(d[i]);
         }
         printf("%lld\n",ans);
     }

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

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