#ifndef call_from_test #include using namespace std; #define call_from_test #ifndef call_from_test #include using namespace std; #endif //BEGIN CUT HERE namespace FFT{ using dbl = double; struct num{ dbl x,y; num(){x=y=0;} num(dbl x,dbl y):x(x),y(y){} }; inline num operator+(num a,num b){ return num(a.x+b.x,a.y+b.y); } inline num operator-(num a,num b){ return num(a.x-b.x,a.y-b.y); } inline num operator*(num a,num b){ return num(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x); } inline num conj(num a){ return num(a.x,-a.y); } int base=1; vector rts={{0,0},{1,0}}; vector rev={0,1}; const dbl PI=asinl(1)*2; void ensure_base(int nbase){ if(nbase<=base) return; rev.resize(1<>1]>>1)+((i&1)<<(nbase-1)); rts.resize(1< &as){ int n=as.size(); assert((n&(n-1))==0); int zeros=__builtin_ctz(n); ensure_base(zeros); int shift=base-zeros; for(int i=0;i>shift)) swap(as[i],as[rev[i]>>shift]); for(int k=1;k multiply(vector &as,vector &bs){ int need=as.size()+bs.size()-1; int nbase=0; while((1< fa(sz); for(int i=0;i>1);i++){ int j=(sz-i)&(sz-1); num z=(fa[j]*fa[j]-conj(fa[i]*fa[i]))*r; if(i!=j) fa[j]=(fa[i]*fa[i]-conj(fa[j]*fa[j]))*r; fa[i]=z; } fft(fa); vector res(need); for(int i=0;i as(n+1,0),bs(n+1,0); for(int i=1;i<=n;i++) scanf("%d %d",&as[i],&bs[i]); auto cs=FFT::multiply(as,bs); for(int i=1;i<=n*2;i++) printf("%lld\n",cs[i]); return 0; } /* verified on 2019/11/26 https://atcoder.jp/contests/atc001/tasks/fft_c */ #endif #undef call_from_test #endif //BEGIN CUT HERE struct bigint { using ll = long long; using vll = vector; constexpr static ll base = 1000000000; constexpr static ll base_digits = 9; vll a; ll sign; bigint():sign(1){} bigint(ll v){*this=v;} bigint(const string &s){read(s);} static bigint add_identity(){return bigint(0);} static bigint mul_identity(){return bigint(1);} void operator=(ll v){ sign=1; if(v<0) sign=-1,v=-v; for(;v>0;v=v/base) a.push_back(v%base); } bigint operator+(const bigint &v) const{ if(sign==v.sign){ bigint res=v; for(ll i=0,carry=0;i<(ll)max(a.size(),v.a.size())||carry;++i){ if(i==(ll)res.a.size()) res.a.push_back(0); res.a[i]+=carry+(i<(ll)a.size()?a[i]:0); carry=res.a[i]>=base; if(carry) res.a[i]-=base; } return res; } return *this -(-v); } bigint operator-(const bigint &v) const{ if(sign==v.sign){ if(abs()>=v.abs()){ bigint res=*this; for(ll i=0,carry=0;i<(ll)v.a.size()||carry;++i){ res.a[i]-=carry+(i<(ll)v.a.size()?v.a[i]:0); carry=res.a[i]<0; if(carry) res.a[i]+=base; } res.trim(); return res; } return -(v-*this); } return *this+(-v); } void operator*=(ll v){ if(v<0) sign=-sign,v=-v; for(ll i=0,carry=0;i<(ll)a.size()|| carry;++i){ if(i ==(ll)a.size()) a.push_back(0); ll cur=a[i] *(ll)v+carry; carry=(ll)(cur/base); a[i]=(ll)(cur%base); // asm("divl %%ecx" : "=a"(carry),"=d"(a[i]) : "A"(cur),"c"(base)); } trim(); } bigint operator*(ll v) const{ bigint res=*this; res *=v; return res; } friend pair divmod(const bigint &a1,const bigint &b1){ ll norm=base/(b1.a.back()+1); bigint a=a1.abs()*norm; bigint b=b1.abs()*norm; bigint q,r; q.a.resize(a.a.size()); for(ll i=a.a.size()-1;i>=0;i--){ r *=base; r+=a.a[i]; ll s1=r.a.size()<=b.a.size() ? 0 : r.a[b.a.size()]; ll s2=r.a.size()<=b.a.size()-1 ? 0 : r.a[b.a.size()-1]; ll d=((ll)base*s1+s2)/b.a.back(); r-=b*d; while(r<0) r+=b,--d; q.a[i]=d; } q.sign=a1.sign*b1.sign; r.sign=a1.sign; q.trim(); r.trim(); return make_pair(q,r/norm); } bigint operator/(const bigint &v) const{ return divmod(*this,v).first; } bigint operator%(const bigint &v) const{ return divmod(*this,v).second; } void operator/=(ll v){ if(v<0) sign=-sign,v=-v; for(ll i=(ll)a.size()-1,rem=0;i>=0;--i){ ll cur=a[i]+rem *(ll)base; a[i]=(ll)(cur/v); rem=(ll)(cur%v); } trim(); } bigint operator/(ll v) const{ bigint res=*this; res/=v; return res; } ll operator%(ll v) const{ if(v<0) v=-v; ll m=0; for(ll i=a.size()-1;i>=0;--i) m=(a[i]+m *(ll)base)%v; return m*sign; } void operator+=(const bigint &v){ *this=*this+v; } void operator-=(const bigint &v){ *this=*this-v; } void operator*=(const bigint &v){ *this=*this*v; } void operator/=(const bigint &v){ *this=*this/v; } bool operator<(const bigint &v) const{ if(sign!=v.sign) return sign=0;i--) if(a[i]!=v.a[i]) return a[i]*sign(const bigint &v) const{ return v<*this; } bool operator<=(const bigint &v) const{ return !(v<*this); } bool operator>=(const bigint &v) const{ return !(*this=0;i--) res=res*base+a[i]; return res*sign; } friend bigint gcd(const bigint &a,const bigint &b){ return b.isZero() ? a : gcd(b,a%b); } friend bigint lcm(const bigint &a,const bigint &b){ return a/gcd(a,b)*b; } void read(const string &s){ sign=1; a.clear(); ll pos=0; while(pos<(ll)s.size()&&(s[pos]=='-'|| s[pos]=='+')){ if(s[pos]=='-') sign=-sign; ++pos; } for(ll i=s.size()-1;i>=pos;i-=base_digits){ ll x=0; for(ll j=max(pos,i-base_digits+1);j<=i;j++) x=x*10+s[j]-'0'; a.push_back(x); } trim(); } friend istream &operator>>(istream &stream,bigint &v){ string s; stream>>s; v.read(s); return stream; } friend ostream &operator<<(ostream &stream,const bigint &v){ if(v.sign==-1) stream<<'-'; stream<<(v.a.empty()?0:v.a.back()); for(ll i=(ll)v.a.size()-2;i>=0;--i) stream<=new_digits){ res.push_back(signed(cur%p[new_digits])); cur/=p[new_digits]; cur_digits-=new_digits; } } res.push_back((signed)cur); while(!res.empty()&&!res.back()) res.pop_back(); return res; } static vll karatsubaMultiply(const vll &a,const vll &b){ ll n=a.size(); vll res(n+n); if(n<=32){ for(ll i=0;i>1; vll a1(a.begin(),a.begin()+k); vll a2(a.begin()+k,a.end()); vll b1(b.begin(),b.begin()+k); vll b2(b.begin()+k,b.end()); vll a1b1=karatsubaMultiply(a1,b1); vll a2b2=karatsubaMultiply(a2,b2); for(ll i=0;ia,base_digits,nbase_digits); vll b=convert_base(v.a,base_digits,nbase_digits); //*/ while(a.size()0) c.push_back(0); } res.a=convert_base(res.a,nbase_digits,base_digits); res.trim(); return res; } }; //END CUT HERE #ifndef call_from_test #define call_from_test #ifndef call_from_test #include using namespace std; #endif //BEGIN CUT HERE template struct BIT{ int n; vector bit; //1-indexed BIT(int n_):n(n_+1),bit(n,0){} T sum(int i){ T s(0); for(int x=i;x>0;x-=(x&-x)) s+=bit[x]; return s; } void add(int i,T a){ if(i==0) return; for(int x=i;x<=n;x+=(x&-x)) bit[x]+=a; } int lower_bound(int w){ if(w<=0) return 0; int x=0,r=1; while(r0;k>>=1){ if(x+k<=n&&bit[x+k]>q; BIT bit(2e5); while(q--){ int t,x; cin>>t>>x; if(t==1) bit.add(x,1); if(t==2){ int k=bit.lower_bound(x); bit.add(k,-1); cout< using namespace std; #endif //BEGIN CUT HERE template struct SquareMatrix{ typedef array arr; typedef array mat; mat dat; SquareMatrix(){ for(size_t i=0;i>=1; } return res; } }; //END CUT HERE #ifndef call_from_test #define call_from_test #include "../tools/fastio.cpp" #include "../tools/compress.cpp" #include "../mod/mint.cpp" #include "../segtree/basic/ushi.cpp" #undef call_from_test //INSERT ABOVE HERE signed YAHOO2019_FINAL_D(){ using ll = long long; using M = Mint; using MM = SquareMatrix; int n,q; cin>>n>>q; vector ts(q),ls(q),rs(q),ps(q); for(ll i=0;i>ts[i]; if(ts[i]==1) cin>>ps[i]; if(ts[i]==2) cin>>ls[i]>>rs[i]; } vector vs; for(ll i=0;i vt(vs.size()-1,A); for(ll i=0;i+1<(ll)vs.size();i++){ vt[i]=A.pow(vs[i+1]-vs[i]); } MM I=MM::mul_identity(); auto f=[&](MM a,MM b){return b*a;}; SegmentTree seg(f,I); seg.build(vt); vector used(vs.size(),0); for(ll i=0;i>s; struct M{ uint32_t v; M(){*this=add_identity();} M(uint32_t v):v(v){} M operator+(const M &a)const{return M(v+a.v);} M operator*(const M &a)const{return M(v*a.v);} static M add_identity(){return M(0);} static M mul_identity(){return M(1);} }; using SM = SquareMatrix; auto f=[](SM a,SM b){return a*b;}; SM ti=SM::mul_identity(); SegmentTree seg(f,ti); vector vt(s.size(),ti); for(int i=0;i<(int)s.size();i++){ if(s[i]=='D') vt[i][0][1]=1; if(s[i]=='I') vt[i][1][2]=1; if(s[i]=='S') vt[i][2][3]=1; if(s[i]=='C') vt[i][3][4]=1; if(s[i]=='O') vt[i][4][5]=1; } seg.build(vt); int q; cin>>q; for(int i=0;i>l>>r; l--; cout<>n>>m; n/=1000; n%=m; bigint res(1); for(int i=0;i>l; if(l==2){ cout<<3<; M A; A[0][0]=1;A[0][1]=1; A[1][0]=1;A[1][1]=0; if(l&1) cout< dfs(int l,int r){ if(l+1>=r) return make_pair(bigint(dp[l]),bigint(l+1)); int m=(l+r)>>1; auto v=dfs(l,m); auto t=dfs(m,r); t.first*=v.second; t.first+=v.first; t.second*=v.second; return t; } signed YUKI_696(){ int n; cin>>n; vector p(n); for(int i=0;i>p[i]; BIT bit(n+1); for(int i=0;i