排名【二分答案】【数据结构+模拟】

【一句话题意】给定n个人的初始位置，再给出每个人的速度，问在无限长的跑道上前m次“追赶事件”(一个人超过另一个人）n,m<=1e5

【分析】
这道题二分答案乱搞，先二分时间，再归并排序求逆序对，然后有个不强和另一个超强的剪枝。（弱）如果二分出的mid，求得的逆序对个数等于m，则跳出循环。（强）如果逆序对个数已经大于m了，则结束归并排序。

std是说，用链表模拟（或者说维护？）相对位置，再用优先队列维护下一次超越事件，好像就可以了。

【code】
考场上二分答案算法，没加剪枝果断超时。

#pragma GCC optimize(2)
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define lb long double
using namespace std;
const int maxn=2e5+1000;
const lb eps=1e-9;
int n,m;
template <typename T>inline void read(T &x){
	x=0;T fl=1;char tmp=getchar();
	while(tmp<'0'||tmp>'9') {if(tmp=='-')fl=-fl;tmp=getchar();}
	while(tmp>='0'&&tmp<='9') x=(x<<1)+(x<<3)+tmp-'0',tmp=getchar();
	x=x*fl;
}
struct node{int a,b,id;}p[maxn];
struct node1{int id;lb val;}p1[maxn],p2[maxn];
struct node2{int a,b;lb t;}ans[maxn];
lb a[maxn],b[maxn];
int cnt,tot;
void marge(int l,int r){
	if(l==r) return ;
	int mid=l+r>>1;
	marge(l,mid),marge(mid+1,r);
	if(cnt>m) return ;
	int i=l,j=mid+1,k=l;
	while(i<=mid&&j<=r){
		if(a[i]>a[j]) b[k++]=a[i++];
		else cnt+=(mid-i+1),b[k++]=a[j++];
	}
	while(i<=mid) b[k++]=a[i++];
	while(j<=r) b[k++]=a[j++];
	for(int i=l;i<=r;i++) a[i]=b[i];
}
void check(lb x){
	for(int i=1;i<=n;i++)
		a[i]=(lb)p[i].a+p[i].b*x;
	cnt=0,marge(1,n);//nlogn
	return ;
}
void smarge(int l,int r){
	if(l==r) return ;
	int mid=l+r>>1;
	smarge(l,mid),smarge(mid+1,r);
	int i=l,j=mid+1,k=l;
	while(i<=mid&&j<=r){
		if(p1[i].val>p1[j].val) p2[k++]=p1[i++];
		else {
			for(int e=i;e<=mid;e++)
				ans[++tot].a=p1[j].id,ans[tot].b=p1[e].id;
			p2[k++]=p1[j++];
		}
	}
	while(i<=mid) p2[k++]=p1[i++];
	while(j<=r) p2[k++]=p1[j++];
	for(int i=l;i<=r;i++) p1[i]=p2[i];
}
bool cmp1(const node x,const node y){return x.a>y.a;}
bool cmp2(const node x,const node y){return x.id<y.id;}
bool cmp3(const node2 x,const node2 y){
	if(fabs(x.t-y.t)<=eps){
		if(x.a==y.a) return x.b<y.b;
		else return x.a<y.a;
	}
	return x.t<y.t;
}
void solve(lb x){
	for(int i=1;i<=n;i++)
		p1[i].val=(lb)p[i].a+p[i].b*x,p1[i].id=p[i].id;
	cnt=0,smarge(1,n);
	sort(p+1,p+n+1,cmp2);
	for(int i=1;i<=m;i++)
		ans[i].t=(lb)(p[ans[i].b].a-p[ans[i].a].a)/(p[ans[i].a].b-p[ans[i].b].b);
	sort(ans+1,ans+m+1,cmp3);
	for(int i=1;i<=m;i++)
		printf("%d %d\n",ans[i].a,ans[i].b);
}
int main(){
	cin>>n;	for(int i=1;i<=n;i++) read(p[i].a),read(p[i].b),p[i].id=i; cin>>m;
	sort(p+1,p+n+1,cmp1);
	lb l=0,r=1e9,mid;
	while(cnt!=m){
		mid=(l+r)/2,check(mid);
		if(cnt==m) l=r=mid;
		else if(cnt>m) r=mid;
		else l=mid;
	}
	solve(r);
	return 0;
}