Codeforces 286E Ladies' Shop

大致题意：

让你构造尽量少的数，

使得从这些数中， 选出若干个数（可以重复选同一个数，和<=m）它们的和组成的集合，

与输入的n个数的集合完全相同

n,m<=1e6:

本题做法：

首先我们能轻易得出一个结论，

构造的这些数，一定出现在给定的n个数之中。

这个结论是显然的，这里不再赘述。

于是我们得到一个新的题意：

从n个数中选出若干个数，使得它们能够组成的和的集合与原序列完全相同。

以下两个算法都基于这个结论：

算法1（来自cwy的一个tle做法）：

考虑贪心。从小到大枚举输入的n个数，

如果当前数无法用之前选中的数的和来表示，

那么答案中一定包含这个数，即我们需要选中这个数。

那么，我们要维护一个背包，维护当前已选中的数能够构成的<=m的和。

时间复杂度nm。

算法2：

我们考虑，由于要从这n个数中选择若干个数，

那么未被选择的数，一定能用被选择的数的和来表示，所以才能被删去。

我们来考虑一种做法：

对于原序列中任意一个a,b（a和b可以相等），

若a+b存在，那么a+b应当被删去。

若a+b<=m，但a+b不存在，那么无解（这是显然的）

我们来证明该做法的正确性。

无解的情况很好判断，以下证明的前提是该序列有解。

那么有解的前提是，对于序列中任意一个a,b，若a+b<=m，那么a+b存在于序列中（这是显然的）

如果对于原序列中一个数s(s<=m)，该数能够用序列中某一些数的和来表示，

设s=a0+a1+a2+....+an=a0+(a1+a2+..+an) 。由于有解，a0, a1+a2+..+an也存在于序列中。

那么该数就能够用序列中某两个数的和来表示。

所以判断s是否应当被删去的方法是，判断s是否能表示成a+b，使得a,b存在于序列中。

说了这么多，感觉这个结论还是很显然。

那么关键来了，怎么实现这个算法？

我们需要把原序列中任意两个数的和全部标记起来，然后一一从原序列中删去。

这就是fft模板题了吧。

时间复杂度mlogm

代码：

#include <cmath>
#include <cstdio>
#include <complex>
#include <iostream>
using namespace std;
#define ref(i,x,y)for(int i=x;i<=y;++i)
int read()
{
    char c=getchar();int d=0,f=1;
    for(;c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
    for(;c>='0'&&c<='9';d=d*10+c-48,c=getchar());
    return d*f;
}
const int N=1<<21;
const double pi=acos(-1);
int rev[N],sz,SZ,n,m,a[N];
bool ans[N];
typedef complex<double> com;
class array{
public:
    com a[N];
    void fft(int tp)
    {
        ref(i,0,SZ-1)if(rev[i]<i)swap(a[rev[i]],a[i]);
        for(int i=2;i<=SZ;i<<=1)
        {
            int I=i>>1;com w(cos(pi/I),tp*sin(pi/I));
            for(int j=0;j<SZ;j+=i)
            {
                com W(1,0);
                for(int k=j;k<j+I;++k,W=W*w)
                {
                    com A=a[k],B=a[k+I];
                    a[k]=A+B*W;a[k+I]=A-B*W;
                }
            }
        }
        if(tp==-1)ref(i,0,SZ-1)a[i]/=SZ;
    }
}A;
int main()
{
    n=read(),m=read();
    ref(i,1,n)a[i]=read();
    while((1<<sz)<m)++sz;
    ++sz;SZ=1<<sz;
    ref(i,0,SZ-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(sz-1));
    ref(i,1,n)A.a[a[i]]=com(1,0),ans[a[i]]=1;
    A.fft(1);
    ref(i,0,SZ-1)A.a[i]=A.a[i]*A.a[i];
    A.fft(-1);
    ref(i,1,m)if((int)(A.a[i].real()+0.5)>=1)
    {
        if(!ans[i])
        {
            cout<<"NO"<<endl;
            return 0;
        }
        ans[i]=0;
    }
    cout<<"YES"<<endl;
    int sum=0;
    ref(i,1,m)if(ans[i])sum++;
    printf("%d
",sum);
    ref(i,1,m)if(ans[i])printf("%d ",i);
    printf("
");
}