結果
問題 | 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";}