codeforces286E Ladies' Shop

题目大意：给你n个互不相同的数，保证其不大于m。要求你再确定一串数p，满足：1、p个数可以选出一些集合（相同的数可重复选），其和可以组成所有给定的n个数。2、从p中选择任意数集（相同的数可重复选），若其和不大于m，那么这个和必在给定的n个数中。3、在满足1和2的情况下p最小。

n,m<=10^6 , 8s

（若不是tag上友情提示是fft，蒟蒻怎么也不会想到正解）

不难发现p数一定是n中的数（废话），那么n中的任意一个数的任意倍（不大于m）一定也在这n个数中——否则no solution。

这样，我们构造一个多项式：k0*x0+k1*x1+k2*x2+...+km*xm,如果i在这n个数中，那么ki=1,否则ki=0,特别的k0=1

求这个多项式的平方，我们得到的多项式有什么意义呢？——在p个数中选出两个数求和，ki表示和为i的方案数*2。

选出两个数是这样，我们是不是要求三次方、四次方乃至更多的呢？

不要！其实第一步我们保证了n中任意一个数t，若kt<=m,那么kt在这n个数中。那么也就是说对于n中的两个数a,b,若pa+qb<=m，那么pa+qb前的系数必定大于1*2

最后我们只要看平方得到的系数ki，若i未在原来的n个数中出现且ki>1*2，no solution。

若i在原来的n个数中出现且ki>1*2，最后的p个数中不需要此数，否则一定需要。

求平方就用fft吧。。。。

（好题点赞！）

#include 
#include 
#include 
#include 

using namespace std;
const int Maxn=4000005;
int n,m,N,M,i,j,k,t,tmp,dig[25],rev[Maxn],ans;
bool vis[Maxn],v[Maxn];
struct CP
{
  double x,y;
  CP operator +(const CP &a)const{ return (CP){x+a.x,y+a.y}; }
  CP operator -(const CP &a)const{ return (CP){x-a.x,y-a.y}; }
  CP operator *(const CP &a)const{ return (CP){x*a.x-y*a.y,x*a.y+y*a.x}; }
} a[Maxn];

void FFT(CP a[],int flag){
  for (i=0;i0;j>>=1) dig[len++]=(j&1);
  	for (j=0;j0?' ':'\n');
  }
  return 0;
}