【noi.ac #889】背包

题目

题目描述
有 nn 个物品排成一行，每个物品有一个价值，所有的价值都是小于 13 的正整数。

有 qq 次操作。每次操作输入 l,r(1≤l≤r≤n)l,r(1≤l≤r≤n)，你需要对于第 ll 个至第 rr 个物品：

对于 k=0,1,2,…,12k=0,1,2,…,12 输出：从这些物品中任选若干个（可以是0个），总价值 mod13=kmod13=k 的方案数，对 2 取模；
将第 ll 个物品的价值增加 1（如果原价值是 12 则将其价值修改为 1）；
将第 rr 个物品的价值增加 1（如果原价值是 12 则将其价值修改为 1）。
输入格式
从标准输入读入数据。

第一行输入两个整数 n,qn,q。

接下来一行，输入一个字符串，依次表示 nn 个物品的价值。这里用 a 表示 10，b 表示 11，c 表示 12。

接下来 qq 行，每行输入两个整数 l,rl,r。

输出格式
输出到标准输出。

对于每次操作，输出一个长度为 13 的 01 串，依次表示 k=0,1,2,…,12k=0,1,2,…,12 时的答案。

样例
输入
10 10
a2a7b44c48
4 7
6 9
4 10
9 10
8 9
3 5
9 9
1 7
2 5
2 4
输出
0010011110010
0000110110000
0001011001101
1100010001000
1000011000001
1000101111110
1000000100000
1010001111000
0110001100000
1101000101111
子任务
最终测试数据中，所有的初始价值与所有的询问都是独立均匀随机生成的。

测试点 n q
1 100 100
2 1000 1000
3 10000 10000
4 125000 100000
5 250000 100000
6 500000 100000
7 1000000 100000
8 2000000 100000
9 4000000 100000
10 8000000 100000

思路

直接用压位线段树来维护区间的答案。两个区间合并的方式是可以枚举两边模2为1的位置来更新答案，然后套用线段树维护即可。

代码

#include<bits/stdc++.h> 
using namespace std;
const int M=1<<20,N=8e6+77;
int calc(int x,int y)
{
	int z=0;
	for(int i=0;y;i++,y>>=1)if(y&1)z^=x<<i;
	return (z>>13)^(z&8191);
}
int b[M],c[M],f[M],m,n,q,s[N],t;char a[N];
void query(int l,int r)
{
	int x=1;
	for(l+=m-1,r+=m+1;l^r^1;l>>=1,r>>=1)
	{
		if(~l&1)x=calc(x,f[l^1]);
		if(r&1)x=calc(x,f[r^1]);
	}
	for(int i=0;i<13;i++)printf("%d",x&1),x>>=1;
}
inline void add(int u){for(u|=m,f[u]=f[u]>>12?3:f[u]<<1^3;u>>=1;f[u]=calc(f[u<<1],f[u<<1|1]));}
int main()
{
	scanf("%d%d%s",&n,&q,a+1);
	for(int i=1;i<=q;i++)scanf("%d%d",b+i,c+i),s[++t]=b[i],s[++t]=b[i]+1,s[++t]=c[i],s[++t]=c[i]+1;
	s[++t]=1,sort(s+1,s+t+1);
	for(int i=m=t=1;i<=(q<<2|1);i++)if(s[i]!=s[t])s[++t]=s[i];
	if(s[t]>n)t--;
	for(int i=1;i<=q;i++)b[i]=lower_bound(s+1,s+t+1,b[i])-s,c[i]=lower_bound(s+1,s+t+1,c[i])-s;
	for(s[t+1]=0;m<t+2;)m<<=1;
	for(int i=1,j=0;i<=n;f[m|j]=calc(f[m|j],1<<(a[i]<60?a[i]^48:a[i]-87)|1),i++)if(s[j+1]==i)f[m|++j]=1;
	for(int i=m-1;i;i--)f[i]=calc(f[i<<1],f[i<<1|1]);
	for(int i=1;i<=q;i++)query(b[i],c[i]),puts(""),add(b[i]),add(c[i]);
	return 0;
}