結果
問題 | No.2506 Sum of Weighted Powers |
ユーザー |
|
提出日時 | 2023-10-14 00:48:52 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 191 ms / 2,000 ms |
コード長 | 994 bytes |
コンパイル時間 | 2,213 ms |
コンパイル使用メモリ | 112,484 KB |
実行使用メモリ | 14,868 KB |
最終ジャッジ日時 | 2024-09-15 20:35:47 |
合計ジャッジ時間 | 5,853 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 42 |
ソースコード
#include<iostream> #include<vector> #include<algorithm> #include<cassert> #include<atcoder/modint> #include<atcoder/convolution> using namespace std; using mint=atcoder::static_modint<943718401>; int N; mint X,A[2<<17],B[2<<17],C[2<<17]; istream&operator>>(istream&is,mint&v) { int val; is>>val; v=mint::raw(val); return is; } mint TX[2<<17],TXinv[2<<17]; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cin>>N>>X; for(int i=0;i<=N;i++)cin>>A[i]; for(int i=0;i<=N;i++)cin>>B[i]; for(int i=0;i<=N;i++)cin>>C[i]; if(X.val()==0) { mint ans=0; for(int j=0;j<=N;j++)ans+=A[j]*B[j]*C[0]; for(int k=1;k<=N;k++)ans+=A[k]*B[0]*C[k]; cout<<ans.val()<<endl; return 0; } for(int i=0;i<=N;i++) { long t=(long)i*(i+1)*(i-1)/3; TX[i]=X.pow(t); TXinv[i]=TX[i].inv(); B[i]*=TXinv[i]; C[i]*=TXinv[i]; } vector<mint>b(B,B+N+1),c(C,C+N+1); vector<mint>a=atcoder::convolution(b,c); mint ans=0; for(int i=0;i<=N;i++)ans+=A[i]*TX[i]*a[i]; cout<<ans.val()<<"\n"; }