CC Killing Monsters 分块+DP(SoS)

题意:n个怪兽 第i只的血量为h[i].Q次操作,第i次操作为(x,y) 对所有满足i&x=i的下标i 其血量减少y.
n<=2^17 ,q<=2^18 .x,y,h[i]<=1e9. 每次操作后输出还有多少只怪兽其h[i]>0.

给出x  x&i==i的下标i 肯定为x的子集  例如x=011110 i= 0????0 ?可以为0或者1.
子集什么,试一试用SoS 假如下标i二进制表示为 101100 
那么对它有伤害的操作其x为 1?11?? i为x的子集 则i为1的位置必须置为1,为0的位置可以为1或者0 

此时F(mask,i)表示 包含mask的集合s有多少个(伤害之和) 并且s和mask的不同位置在前i个位置.

现在将操作分块,对每一块 可以用SOS求出 F[i]：这一块的操作 对下标i造成多少伤害. 
若F[i]>=h[i] 在暴力遍历该块 找出最后一个杀死i的操作. O(sqrt(q)*n + n*bit)

对每个操作记录其干掉多少个人 就能知道每个操作过后还活着的人数.

#include <bits/stdc++.h>
#define rep(i,a,b)	for(int i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int N=(1<<18)+1,BIT=18,inf=1e9+5,mask=(1<<17)-1;
int res[N],x[N],y[N],h[N],n,Q,f[N];
int p=0,block,num;
void SOS()
{
	memset(f,0,sizeof(f));
	for(int k=0;k<block&&p<Q;k++)
		f[x[p]]+=y[p],p++;
		
	for(int bit=0;bit<BIT;bit++)
		for(int k=n-1;k>=0;k--)
			if(!((k>>bit)&1))
				f[k]+=f[k^(1<<bit)],f[k]=min(f[k],inf);
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n;
	rep(i,0,n-1)	cin>>h[i];
	cin>>Q;
	rep(i,0,Q-1)	cin>>x[i]>>y[i],x[i]&=mask;
	block=1000,num=Q/block;
	if(Q%block)	
		num++;
	for(int i=0;i<num;i++)
	{
		SOS();//求每个i在这一块操作中收到的伤害
		for(int k=0;k<n;k++)
		{
			if(h[k]<=0)	continue;
			h[k]-=f[k];
			if(h[k]>0)	continue;
			ll sum=0;
			for(int j=i*block;j<(i+1)*block&&j<Q;j++)
			{
				if((k&x[j])==k)
				{
					sum+=y[j];
					if(sum>=h[k]+f[k])
					{
						res[j]++;
						break;		
					}
				}
			}		
		} 
	}	
	int ans=n;
	rep(i,0,Q-1)
	{
		ans-=res[i];
		cout<<ans<<"\n";
	}
	return 0;
}

まとめ

最开始想着对序列进行分块,然而这题在线对询问进行分块..(离线的莫队也是将询问分块)

SOS DP  简要回顾
给出序列a,求F[mask]：序列a中,下标为mask的子集的和

设S(mask,i) ：mask的子集 并且 子集和mask的不同位置只在前i个位置中产生. (从右往左编号[0..])
若mask的第i位为1 则当前子集第i位可以为0或者1.
S(mask,i)=S(mask,i-1)+S(mask^(2^i),i-1)