BZOJ 2244 拦截导弹 cdq分治

2244: [SDOI2011]拦截导弹

Time Limit: 30 Sec   Memory Limit: 512 MBSec   Special Judge
Submit: 877   Solved: 303
[ Submit][ Status][ Discuss]

Description

某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度、并且能够拦截任意速度的导弹，但是以后每一发炮弹都不能高于前一发的高度，其拦截的导弹的飞行速度也不能大于前一发。某天，雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段，所以只有一套系统，因此有可能不能拦截所有的导弹。

在不能拦截所有的导弹的情况下，我们当然要选择使国家损失最小、也就是拦截导弹的数量最多的方案。但是拦截导弹数量的最多的方案有可能有多个，如果有多个最优方案，那么我们会随机选取一个作为最终的拦截导弹行动蓝图。

我方间谍已经获取了所有敌军导弹的高度和速度，你的任务是计算出在执行上述决策时，每枚导弹被拦截掉的概率。

Input

第一行包含一个正整数n，表示敌军导弹数量；

下面 行按顺序给出了敌军所有导弹信息：

第i+1行包含2个正整数h_i和v_i，分别表示第 枚导弹的高度和速度。

Output

输出包含两行。

第一行为一个正整数，表示最多能拦截掉的导弹数量；

第二行包含n个0到1之间的实数，第i个数字表示第i枚导弹被拦截掉的概率（你可以保留任意多位有效数字）。

Sample Input

4

3 30

4 40

6 60

3 30

Sample Output

2

0.33333 0.33333 0.33333 1.00000

【数据规模和约定】


对于100%的数据，1≤n≤5*104， 1≤hi ，vi≤109；

均匀分布着约30%的数据，所有vi均相等。

均匀分布着约50%的数据，满足1≤hi ，vi≤1000。

HINT

鸣谢kac提供sj程序！

Source

第一轮day2

第一个问题，显然求个二维的LIS就可以。可以转化为三维偏序问题，用cdq分治问题求解。

第二问，实际上就是问某个导弹在多少个最长的LIS序列当中出现过。

观察发现，可以把下标，速度，高度作为一个三元组，求LIS长度就是求三元偏序当中的最值问题，求LIS数量就是三元偏序求和。这样，就可以用cdq分治求解。对于每一层，维护一个树状数组，来得到答案。对于每个位置，正反求两次长度和数量，分别代表此元素向前、向后取最长LIS和此时的LIS数量。最后全部汇总即可。

#include <cstdio>
#include <iostream>
#include <string.h>
#include <algorithm>
#define mem0(a) memset(a,0,sizeof(a))
#define meminf(a) memset(a,0x3f,sizeof(a))
using namespace std;
typedef double db;
const int maxn=100005,inf=0x3f3f3f3f;
int a[maxn],b[maxn],c[maxn];
int num,n;

struct Tree{
	int l,t;
	db g;
};
Tree f[maxn];

struct query {
	int p,a,b,id,l[2];
	db g[2];
	bool operator< (const query &x) const {
		return id<x.id;
	}
};
query q[maxn],tem[maxn];

bool cmp(int x,int y) {
	return q[x].a<q[y].a || (q[x].a==q[y].a&&q[x].id<q[y].id);
}

int lowbit(int x) {
	return x&(-x);
}

Tree gettree(int tt) {
	Tree ans;
	ans.l=0;ans.g=0.0;
	for (int t=tt;t;t-=lowbit(t)) {
		if (f[t].t==num) {
			if (f[t].l>ans.l) {
				ans.g=f[t].g;ans.l=f[t].l;
			} else if (f[t].l==ans.l) ans.g+=f[t].g;
		}
	}
	return ans;
}

void update(int tt,int l,db g) {
	int t=tt;
	for (int t=tt;t<=n;t+=lowbit(t)) {
		if (f[t].t!=num) {
			f[t].t=num;
		    f[t].g=g;f[t].l=l;
		} else {
			if (f[t].l==l) f[t].g+=g; else 
			    if (l>f[t].l) f[t].g=g,f[t].l=l;
		}
	}
}

void cdq(int l,int r,int p) {
	if (l==r) {
		if (q[l].l[p] < 1) 
			q[l].l[p] = 1, q[l].g[p] = 1.0;
		return;
	}
	int mid=(l+r)/2;
	memcpy(tem+l,q+l,sizeof(query)*(r-l+1));
	int h=l,t=mid+1;
    for (int i=l;i<=r;i++)  (tem[i].p<=mid?q[h++]:q[t++])=tem[i];
	cdq(l,mid,p);
	num++;
	h=l,t=mid+1;
	for (int i=mid+1;i<=r;i++)
    {
        while (h<=mid&&q[h].id<q[i].id) 
		    update(q[h].b,q[h].l[p],q[h].g[p]),h++;
        Tree now=gettree(q[i].b);
        if (!now.l) continue;
        if (now.l+1>q[i].l[p])  q[i].l[p]=now.l+1,q[i].g[p]=now.g;
            else if (now.l+1==q[i].l[p]) q[i].g[p]+=now.g;
    }
	cdq(mid+1,r,p);
	merge(q+l,q+mid+1,q+mid+1,q+r+1,tem+l);
	memcpy(q+l,tem+l,sizeof(query)*(r-l+1));
}

int main() {
	int i,j;
	num=0;
	scanf("%d",&n);
	mem0(f);
	for (i=1;i<=n;i++) {
		scanf("%d%d",&a[i],&b[i]);
		q[i].a=a[i];q[i].b=b[i];q[i].id=c[i]=i;
		q[i].l[0]=q[i].l[1]=q[i].g[0]=q[i].g[1]=0;
	}
	sort(a+1,a+n+1);sort(b+1,b+n+1);
	int t1=unique(a+1,a+n+1)-a-1,t2=unique(b+1,b+n+1)-b-1;
	for (i=1;i<=n;i++) {
		q[i].a=t1-(lower_bound(a+1,a+t1+1,q[i].a)-a)+1;
		q[i].b=t2-(lower_bound(b+1,b+t2+1,q[i].b)-b)+1;
	}
	sort(c + 1, c + n + 1, cmp);
	for (i = 1; i <= n; i++) q[c[i]].p = i;
	cdq(1,n,0);
	for (i=1;i<=n;i++) {
		q[i].p=n-q[i].p+1;
		q[i].id=n-q[i].id+1;
		q[i].a=t1-q[i].a+1;
		q[i].b=t2-q[i].b+1;
	}
	reverse(q+1,q+n+1);
	cdq(1,n,1);
	reverse(q+1,q+n+1);
/*	for (i=1;i<=n;i++) {
		cout << q[i].l[0] << ' ' << q[i].l[1] << ' ';
		cout << q[i].g[0] << ' ' << q[i].g[1] << endl;
	}*/
	int ml=-1;db sums=0.0;
	for (i=1;i<=n;i++) {
		if (q[i].l[0]+q[i].l[1]>ml) {
			ml=q[i].l[0]+q[i].l[1];
			sums=0.0;
		}
		if (q[i].l[1]==1&&q[i].l[0]+q[i].l[1]==ml) 
		    sums+=q[i].g[0]*q[i].g[1];
	}
	printf("%d\n",ml-1);
	for (i=1;i<=n;i++) {
		if (q[i].l[0]+q[i].l[1]==ml) {
			printf("%.5f",q[i].g[0]*q[i].g[1]/sums);
		} else printf("0.00000");
		if (i!=n) printf(" ");
	}
//	system("pause");
	return 0;
}

/*
4
3 30
4 40
6 60
3 30
*/