#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define I_MAX 2147483647 #define LL_MAX 9223372036854775807 #define ll long long #define ld long double struct XX{ int l; int r; int m; }; class xxGreater { public: bool operator()(const XX& riLeft, const XX& riRight) const { //第2条件 if((riLeft.r) == (riRight.r)){ return riLeft.l < riRight.l;//<:昇順(小さいものから順番)、>:降順(大きいものから順番) //プライオリティキューの場合は > で、top()すると値の小さいものがとれる } //第1条件 return (riLeft.r) < (riRight.r); } }; //map prime_f(long long n){ // mapres; // for(int i=2;i*i<=n;i++){ // while(n%i==0){ // ++res[i]; // n/=i; // } // } // if(n!=1)res[n]=1; // return res; //} //int n; ////int dat[2*10000000]; ////int dat2[2*10000000]; //int dat[10]; //int dat2[10]; // //void init(int n_){ // n=1; // while(n0){ // k=(k-1)/2; // dat[k]=dat[k*2+1]+dat[k*2+2]; // } //} // ////[a,b)の間を[l,r]区間で比較しアップデート ////引数のindexに注意 ////nは固定。initで計算すみ ////update2(L[i],R[i]+1,0,0,n,D[i]); //void update2(int a,int b,int k,int l,int r,int v){//v更新値、区間は0-index // if(r<=a || b<=l)return; // if(a<=l && r<=b){ // dat[k]+=dat2[k]; // if(r-l>1){ // dat2[k*2+1]+=dat2[k]/2; // dat2[k*2+1]+=dat2[k]/2; // } // dat2[k]=v*(r-l); // return; // }else{ // update2(a,b,k*2+1,l,(l+r)/2,v); // update2(a,b,k*2+2,(l+r)/2,r,v); // return; // } //} // //int query(int a,int b,int k,int l,int r){ // if(r<=a || b<=l)return 0; // if(a<=l && r<=b){ // dat[k]+=dat2[k]; // if(r-l>1){ // dat2[k*2+1]+=dat2[k]/2; // dat2[k*2+1]+=dat2[k]/2; // } // dat2[k]=0; // return dat[k]; // } // else{ // int vl=query(a,b,k*2+1,l,(l+r)/2); // int vr=query(a,b,k*2+2,(l+r)/2,r); // return vl+vr; // } //} //void printb(unsigned int v) { // unsigned int mask = (int)1 << (sizeof(v) * CHAR_BIT - 1); // do putchar(mask & v ? '1' : '0'); // while (mask >>= 1); //} #ifdef DEBUG #else #endif void dft(vector>& p, int inverse) { int sz = (int)p.size(); if (sz == 1)return; vector> x0, x1; for(int i=0;i now = 1; complex zeta = polar(1.0, inverse * (2.0*acos(-1))/sz);//複素数の極形式 角度(2*pi)/n for(int i=0;i multiply(vector> f, vector> g) { vector> p, q; int jisu = 1; while (jisu < f.size()+g.size()){ jisu *= 2; } p.resize(jisu); q.resize(jisu); for(int i=0;i res; for(int i=0;i> a; //ここから //入力高速化 ios::sync_with_stdio(false); cin.tie(0); int n,q; cin>>n>>q; vector>f,g; for(int i=0;i>t1; f.push_back(complex(t1,0.0)); } for(int i=0;i>t1; keisu[n-t1]+=1.0; } for(int i=0;i<=n;i++){ g.push_back(complex(keisu[i],0.0)); } vector ret=multiply(f,g); for(int i=n;i<2*n;i++){ printf("%.0f",ret[i]);//誤差が出るのでこの方法で出力 cout << " "; } cout << endl; //ここまで //cout << "ans" << endl; //cout << " " << "ans" << endl; //printf("%.0f\n",ans);//小数点以下表示なし //printf("%.7f\n",p); return 0; }