Description
给出 \(n\) 个 \(\leq m\) 且不同的数 \(a_1,\dots,a_n\)，现在要求从这 \(n\) 个数中选出最少的数字，满足这 \(n\) 个数字都可以由选出的数字组合成（就是做一个完全背包能做出来），并且任意组合出来的数字，只要不超过 \(m\)，就必须让这个数字在给出的 \(n\) 个数中。问是否可行，如果可行，请求出最少选多少数字。 \(n,m\leq 10^6\)。
Sol
先判断是否可行，再看哪些数可以省略。
求出 \(a\) 数组的生成函数，即构造多项式 \(F(x)=\sum f_i\cdot x^i\)。\(f_i\) 为 \(1\) 当且仅当 \(a\) 数组中出现 \(a_*=i\)
然后求出 \(G(x)=F^2(i)=\sum g_i\cdot x^i\)。如果 \(g_i>0\) 那就说明给出的这 \(n\) 个数可以合成 \(i\) 。
于是就得到了从原来的 \(n\) 个数中拿出 \(0\sim 2\) 个的结果。
然而最多拿出 \(m\) 个。
所以还要继续，用快速幂求得 \(f^m\)。如果多项式快速幂的话，复杂度 \(O(n\log^2n)\)，用多项式ln+多项式exp求的话，复杂度 \(O(n\log n)\)。但是多项式exp常数太大了！
事实上是有只做 \(1\) 次FFT的方法的。
显然如果 \(f_i>0\) 的话，\(g_i>0\)。
那我们只要保证满足 \(f_i=0,g_i>0,i\leq m\) 的 \(i\) 不存在就好了。
如果第一轮不存在这些不合法的，那接下来肯定也不存在。感性理解一下这就相当于构成了一个封闭的集合。
所以只做 \(1\) 次FFT就行了。
然后考虑一下哪些数可以省略
如果一个数 \(i\) 可以被其他数表示出来，那 \(g_i\) 一定 \(>2\)。所以 \(g_i=2\) 的 \(i\) 就是必选的。
时间复杂度 \(O(n\log n)\)。
Sol
#pragma GCC optimize(2)
#include<bits/stdc++.h>
using std::min;
using std::max;
using std::swap;
using std::vector;
typedef double db;
typedef long long ll;
#define pb(A) push_back(A)
#define pii std::pair<int,int>
#define all(A) A.begin(),A.end()
#define mp(A,B) std::make_pair(A,B)
const int N=4e6+5;
const int mod=998244353;
int lim,rev[N];
int n,m,a[N],b[N];
int ksm(int a,int b=mod-2,int ans=1){
while(b){
if(b&1) ans=1ll*ans*a%mod;
a=1ll*a*a%mod;b>>=1;
} return ans;
}
int getint(){
int X=0,w=0;char ch=getchar();
while(!isdigit(ch))w|=ch=='-',ch=getchar();
while( isdigit(ch))X=X*10+ch-48,ch=getchar();
if(w) return -X;return X;
}
void ntt(int *f,int g){
for(int i=1;i<lim;i++) if(i<rev[i]) swap(f[i],f[rev[i]]);
for(int mid=1;mid<lim;mid<<=1){
int tmp=ksm(g,(mod-1)/(mid<<1));
for(int R=mid<<1,j=0;j<lim;j+=R){
int w=1;
for(int k=0;k<mid;k++,w=1ll*w*tmp%mod){
int x=f[j+k],y=1ll*w*f[j+k+mid]%mod;
f[j+k]=(x+y)%mod,f[j+k+mid]=(mod+x-y)%mod;
}
}
} if(g>3)
for(int in=ksm(lim),i=0;i<lim;i++) f[i]=1ll*f[i]*in%mod;
}
signed main(){
n=getint(),m=getint();
for(int i=1;i<=n;i++){
int x=getint();
a[x]=b[x]=1;
}
lim=1;while(lim<=m+m) lim<<=1;
for(int i=1;i<lim;i++) rev[i]=(rev[i>>1]>>1)|(i&1?lim>>1:0);
a[0]=1; ntt(a,3);
for(int i=0;i<lim;i++) a[i]=1ll*a[i]*a[i]%mod;
ntt(a,(mod+1)/3);
for(int i=1;i<=m;i++)
if(a[i] and !b[i]) return printf("NO"),0;
puts("YES"); int tot=0;
for(int i=1;i<=m;i++)
if(a[i]==2) tot++;
printf("%d\n",tot);
for(int i=1;i<=m;i++)
if(a[i]==2) printf("%d ",i);
return 0;
}