CF632E Thief in a Shop
Posted troverld
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了CF632E Thief in a Shop相关的知识,希望对你有一定的参考价值。
XI.CF632E Thief in a Shop
听说这题NTT常见的模数都会被卡?orzorz。
一看就是一个完全背包的样式。当然咯,FFT题当然应该用FFT做呀。我们构建(g(x)=[xin ext{商品价值的集合}])。计算(g^k),答案即为(g^k(x))所有有值的位置。
代码:
#include<bits/stdc++.h>
using namespace std;
const double pi=acos(-1);
int n,m,rev[1<<20],lim=1<<20,lg=20,qwq[1<<20],bs[1<<20];
struct cp{
double x,y;
cp(double u=0,double v=0){x=u,y=v;}
friend cp operator +(const cp &u,const cp &v){return cp(u.x+v.x,u.y+v.y);}
friend cp operator -(const cp &u,const cp &v){return cp(u.x-v.x,u.y-v.y);}
friend cp operator *(const cp &u,const cp &v){return cp(u.x*v.x-u.y*v.y,u.x*v.y+u.y*v.x);}
}f[1<<20],g[1<<20];
inline void FFT(cp *a,int tp){
for(register int i=0;i<lim;i++)if(i<rev[i])swap(a[i],a[rev[i]]);
for(register int md=1;md<lim;md<<=1){
cp rt=cp(cos(pi/md),tp*sin(pi/md));
for(register int stp=md<<1,pos=0;pos<lim;pos+=stp){
cp w=cp(1,0);
for(register int i=0;i<md;i++,w=w*rt){
cp x=a[pos+i],y=w*a[pos+md+i];
a[pos+i]=x+y;
a[pos+md+i]=x-y;
}
}
}
}
inline void mul(int *a,int *b,int *c){
for(register int i=0;i<lim;i++)f[i]=cp(a[i],0);
for(register int i=0;i<lim;i++)g[i]=cp(b[i],0);
FFT(f,1),FFT(g,1);
for(register int i=0;i<lim;i++)f[i]=f[i]*g[i];
FFT(f,-1);
for(register int i=0;i<lim;i++){
c[i]=(int)(f[i].x/lim+0.5);
c[i]=(c[i]>=1);
}
}
inline void ksm(int x){
if(x==1){memcpy(qwq,bs,sizeof(bs));return;}
ksm(x>>1);
mul(qwq,qwq,qwq);
if(x&1)mul(qwq,bs,qwq);
}
int main(){
for(register int i=0;i<lim;i++)rev[i]=(rev[i>>1]>>1)|((i&1)<<(lg-1));
scanf("%d%d",&n,&m);
for(register int i=1,x;i<=n;i++)scanf("%d",&x),bs[x]=true;
ksm(m);
for(register int i=0;i<lim;i++)if(qwq[i])printf("%d ",i);
return 0;
}
以上是关于CF632E Thief in a Shop的主要内容,如果未能解决你的问题,请参考以下文章
CodeForces 632E Thief in a Shop
CodeForces - 632E Thief in a Shop 完全背包
CodeForces 632E Thief in a Shop
Educational Codeforces Round 9 E. Thief in a Shop NTT
codeforces Educational Codeforces Round 9 E - Thief in a Shop