Codeforces956D. Contact ATC

$n leq 100000$个飞机在坐标轴上，给坐标给速度，坐标速度异号，还有一个风速在$[-w,w]$区间，$w$比最小的速度绝对值要小。由于风速不知道，所以问有多少对飞机可能在原点相遇。

思维定势：$frac{x_i}{v_i+v}=frac{x_j}{v_j+v}$，$v$是风速，然后推下去，会推到一个三维偏序。。

没有观察题目性质。这个时间是关于风速单调而连续的，所以只要风速最小和风速最大这两个东西求个逆序对就行了。

似乎卡精度，用了分数。

这种题要写题解感觉最近脑子有点锈。。有没有神犇愿意帮忙除个锈啊QAQ

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
//#include<math.h>
//#include<queue>
//#include<vector>
#include<algorithm>
//#include<iostream>
//#include<assert.h>
using namespace std;

int n,w;
#define maxn 200011

struct frac
{
    int a,b;
    bool operator < (const frac &x) const {return 1ll*a*x.b<1ll*b*x.a;}
    bool operator == (const frac &x) const {return 1ll*a*x.b==1ll*b*x.a;}
};

struct Poi{frac x,y; int z;}p[maxn];
bool cmpx(const Poi &a,const Poi &b) {return b.x<a.x || (a.x==b.x && a.y<b.y);}
frac lisa[maxn]; int li=0;

struct BIT
{
    int a[maxn],n;
    void clear(int m) {n=m;}
    void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]+=v;}
    int query(int x) {int ans=0; for (;x;x-=x&-x) ans+=a[x]; return ans;}
}t;

#define LL long long
int main()
{
    scanf("%d%d",&n,&w);
    for (int i=1,a,b;i<=n;i++)
    {
        scanf("%d%d",&a,&b);
        if (a<0) p[i].x=(frac){-a,b-w},p[i].y=(frac){-a,b+w};
        else p[i].x=(frac){a,w-b},p[i].y=(frac){a,-w-b};
        lisa[++li]=p[i].y;
    }
    sort(lisa+1,lisa+1+li);
    for (int i=1;i<=n;i++) p[i].z=lower_bound(lisa+1,lisa+1+li,p[i].y)-lisa;
    
    sort(p+1,p+1+n,cmpx);
    t.clear(n);
    LL ans=0;
    for (int i=1;i<=n;i++)
    {
        ans+=t.query(p[i].z);
        t.add(p[i].z,1);
    }
    printf("%lld
",ans);
    return 0;
}