結果
問題 | No.1549 [Cherry 2nd Tune] BANning Tuple |
ユーザー | 👑 Nachia |
提出日時 | 2021-05-04 00:54:34 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,776 bytes |
コンパイル時間 | 1,334 ms |
コンパイル使用メモリ | 99,172 KB |
実行使用メモリ | 10,112 KB |
最終ジャッジ日時 | 2024-12-14 21:33:03 |
合計ジャッジ時間 | 70,088 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 12 ms
5,248 KB |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | WA | - |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | TLE | - |
testcase_08 | TLE | - |
testcase_09 | TLE | - |
testcase_10 | TLE | - |
testcase_11 | TLE | - |
testcase_12 | TLE | - |
testcase_13 | TLE | - |
testcase_14 | TLE | - |
testcase_15 | TLE | - |
testcase_16 | TLE | - |
testcase_17 | TLE | - |
testcase_18 | TLE | - |
testcase_19 | WA | - |
ソースコード
#include <iostream> #include <vector> #include <map> #include <algorithm> using namespace std; using ll = long long; using ull = unsigned long long; const ull M = 57384961; ull powm(ull a, ull i) { if(i==0) return 1; ull r=powm(a*a%M,i/2); if(i&1) r=r*a%M; return r; } ull invm(ull a){ return powm(a,M-2); } ull primitiveRoot(){ vector<ull> D; ull m = M-1; for(ull p=2; p*p<=m; p++) if(m%p==0){ D.push_back(p); while(m%p==0) m /= p; } if(m != 1) D.push_back(m); for(ull g=2; g<=M-1; g++){ bool ok = true; for(ull d : D) if(powm(g,(M-1)/d) == 1) ok = false; if(ok) return g; } exit(1); } void NTT(vector<ull>& A, ull g){ int N=A.size(); for(int i=0,j=0; j<N; j++){ if(i<j) swap(A[i],A[j]); for(int k=N>>1; k>(i^=k); k>>=1); } for(int i=1; i<N; i<<=1){ ull q=powm(g,(M-1)/i/2), qj=1; for(int j=0; j<i; j++){ for(int k=j; k<N; k+=i*2){ ull l=A[k],r=A[k+i]*qj%M; A[k]=(l+r); if(A[k]>=M) A[k]-=M; A[k+i]=(l+M-r); if(A[k+i]>=M) A[k+i]-=M; } qj=qj*q%M; } } } const int maxT = 3000; const int NTTZ = 1<<13; const ull NTTg = primitiveRoot(); const ull NTTinvZ = invm(NTTZ); const ull NTTinvg = invm(NTTg); vector<ull> convolute(vector<ull> l, vector<ull> r){ NTT(l,NTTg); NTT(r,NTTg); vector<ull> ans(NTTZ,0); for(int i=0; i<NTTZ; i++) ans[i] = l[i] * r[i] % M * NTTinvZ % M; NTT(ans,NTTinvg); for(int i=maxT+1; i<NTTZ; i++) ans[i] = 0; return move(ans); } ull N; int Q; map<ull,int> Idx; vector<vector<ull>> VecList; vector<ull> AnsArr; int main(){ cin>>N>>Q; for(int q=0; q<Q; q++){ ull K; cin>>K; int A,B,S,T; cin>>A>>B>>S>>T; if(Idx.count(K) == 0){ Idx[K] = VecList.size(); VecList.push_back(vector<ull>(NTTZ,0)); for(int i=0; i<=3000; i++) VecList.back()[i] = 1; } int idx = Idx[K]; for(int x=A; x<=B; x++) VecList[idx][x] = 0; AnsArr = vector<ull>(NTTZ,0); AnsArr[0] = 1; for(auto& vec : VecList) AnsArr = convolute(AnsArr, vec); ull leftN = N - VecList.size(); vector<ull> ansCoeff(NTTZ+1); ansCoeff[0] = 1; for(int i=1; i<=NTTZ; i++) ansCoeff[i] = ansCoeff[i-1] * ((leftN+i-1) % M) % M * invm(i) % M; reverse(ansCoeff.begin(),ansCoeff.end()); ansCoeff.push_back(0); reverse(ansCoeff.begin(),ansCoeff.end()); for(int i=0; i<=NTTZ; i++) ansCoeff[i+1] = (ansCoeff[i] + ansCoeff[i+1]) % M; ull ans = 0; for(int s=0; s<NTTZ; s++){ ull coeff = ansCoeff[max(0,T-s+1)] + M - ansCoeff[max(0,S-s)]; ans = (ans + AnsArr[s] * coeff) % M; } cout << (ans % M) << "\n"; } return 0; } struct ios_do_not_sync{ ios_do_not_sync(){ ios::sync_with_stdio(false); cin.tie(nullptr); } } ios_do_not_sync_instance;