1. 程式人生 > >CF924D Contact ATC(數學+BIT求逆序對)

CF924D Contact ATC(數學+BIT求逆序對)

這題妙呀。
我們設風速為w時,飛機i通過原點的時間為tx,
風速為-w時,飛機j通過原點的時間為ty。
則如果存在某個風速w,使得i和j能夠一個時間通過原點,需要滿足(txi-txj)*(tyi-tyj)<=0.這是什麼!這就是二維偏序!
所以我們算出tx,ty,求個二維偏序就好了。

不過這題大概卡精度,我們要手寫分數。
還要注意相等的時候誰在前面的問題。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <utility>
using namespace
std; #define ll long long #define inf 0x3f3f3f3f #define N 100010 inline char gc(){ static char buf[1<<16],*S,*T; if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;} return *S++; } inline ll read(){ ll x=0,f=1;char ch=gc(); while(ch<'0'||ch>'9'){if(ch=='-'
)f=-1;ch=gc();} while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc(); return x*f; } inline int gcd(int x,int y){return y?gcd(y,x%y):x;} int n,W,m,c[N];ll ans=0; struct frac{ int a,b; frac(int _a,int _b){int g=gcd(_a,_b);a=_a/g;b=_b/g;} frac(){} friend bool operator<(frac a,frac b){ return
(ll)a.a*b.b<(ll)b.a*a.b; }friend bool operator==(frac a,frac b){ return a.a==b.a&&a.b==b.b; } }aa[N]; struct hhh{ frac x,y;int pos; friend bool operator<(hhh a,hhh b){ return a.x==b.x?a.pos>b.pos:a.x<b.x; } }a[N]; inline void add(int x,int val){for(;x<=m;x+=x&-x) c[x]+=val;} inline int ask(int x){int res=0;for(;x;x-=x&-x) res+=c[x];return res;} int main(){ // freopen("a.in","r",stdin); n=read();W=read(); for(int i=1;i<=n;++i){ int x=read(),v=read(); a[i].x=frac(abs(x),abs(v+W));a[i].y=frac(abs(x),abs(v-W));aa[i]=a[i].y; }sort(aa+1,aa+n+1);m=unique(aa+1,aa+n+1)-aa-1; for(int i=1;i<=n;++i) a[i].pos=lower_bound(aa+1,aa+m+1,a[i].y)-aa; sort(a+1,a+n+1); for(int i=1;i<=n;++i){ ans+=ask(m-a[i].pos+1);add(m-a[i].pos+1,1); }printf("%I64d\n",ans); return 0; }