POJ2274

这真的是一道数据结构的好题。

题意是在一条直线上有n辆车,每辆车有一个初始位置x[i]和速度v[i],问最终(在无限时间后)一共会发生多少次超车事件(mod 1000000),以及输出这些事件(如果大于10000次输出前10000次)

对于第一问,很明显的求逆序对,如果满足x[i]<x[j]且v[i]>v[j],那么就会产生超车事件

所以树状数组即可

第二问则比较复杂,但我们可以发现每次超车发生时都是相邻的两车

所以我们发现这个可以用来维护

要注意一下的是当事件发生的时间相同时,要按编号大小进行输出(WA死在这里)

只要堆的 bool operator 重载一下即可

CODE

#include<cstdio>
#include<queue>
using namespace std;
typedef double DB;
const int N=250005,mod=1000000;
int tree[105],n,x[N],v[N],num[N],id[N],t,ans;
struct data
{
	int x,y;
	DB t;
	bool operator <(const data &s) const
	{
		if (s.t<t) return 1;
		if (s.t>t) return 0;
		return num[s.x]<num[x];
	}
};
priority_queue <data> heap;
inline char tc(void)
{
	static char fl[100000],*A=fl,*B=fl;
	return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
	x=0; char ch=tc();
	while (ch<'0'||ch>'9') ch=tc();
	while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline void write(int x)
{
	if (x/10) write(x/10);
	putchar(x%10+'0');
}
inline int lowbit(int x)
{
	return x&(-x);
}
inline void add(int x)
{
	while (x)
	{
		++tree[x];
		x-=lowbit(x);
	}
}
inline int get(int x)
{
	int res=0;
	while (x<100)
	{
		res=(res+tree[x])%mod;
		x+=lowbit(x);
	}
	return res;
}
inline void swap(int a,int b)
{
	int temp=x[a]; x[a]=x[b]; x[b]=temp;
	DB t=v[a]; v[a]=v[b]; v[b]=t;
	temp=id[a]; id[a]=id[b]; id[b]=temp;
}
int main()
{
	register int i;
	//freopen("CODE.in","r",stdin); //freopen("CODE.out","w",stdout);
	read(n); 
	for (i=1;i<=n;++i)
	{
		read(x[i]); read(v[i]); num[i]=id[i]=i;
		add(v[i]); ans=(ans+get(v[i]+1))%mod;
		if (i^1&&v[i-1]>v[i]) heap.push((data){i-1,i,(double)(x[i]-x[i-1])/(v[i-1]-v[i])});
	}
	write(ans); putchar('
'); 
	while (t<10000&&!heap.empty())
	{
		int X=heap.top().x,Y=heap.top().y;
		DB T=heap.top().t; heap.pop();
		if (num[X]+1!=num[Y]) continue;
		write(X); putchar(' '); write(Y); putchar('
'); ++t;
		int now1=num[X],now2=num[Y]; num[X]=now2; num[Y]=now1;
		swap(now1,now2);
		if (now1>1&&v[now1-1]>v[now1]) if (num[id[now1-1]]+1==num[id[now1]]) heap.push((data){id[now1-1],id[now1],(double)(x[now1]-x[now1-1])/(v[now1-1]-v[now1])});
		if (now2<n&&v[now2]>v[now2+1]) if (num[id[now2]]+1==num[id[now2+1]]) heap.push((data){id[now2],id[now2+1],(double)(x[now2+1]-x[now2])/(v[now2]-v[now2+1])});
	}
	return 0;
}
原文地址:https://www.cnblogs.com/cjjsb/p/8831774.html