結果
問題 | No.2605 Pickup Parentheses |
ユーザー | TKTYI |
提出日時 | 2024-01-13 10:22:48 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,420 ms / 2,000 ms |
コード長 | 3,216 bytes |
コンパイル時間 | 4,640 ms |
コンパイル使用メモリ | 290,424 KB |
実行使用メモリ | 129,420 KB |
最終ジャッジ日時 | 2024-09-30 06:37:29 |
合計ジャッジ時間 | 86,895 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1,344 ms
120,496 KB |
testcase_01 | AC | 1,247 ms
120,504 KB |
testcase_02 | AC | 1,259 ms
120,508 KB |
testcase_03 | AC | 1,244 ms
120,644 KB |
testcase_04 | AC | 1,263 ms
120,596 KB |
testcase_05 | AC | 1,241 ms
120,600 KB |
testcase_06 | AC | 1,250 ms
120,544 KB |
testcase_07 | AC | 1,240 ms
120,752 KB |
testcase_08 | AC | 1,237 ms
120,672 KB |
testcase_09 | AC | 1,255 ms
120,572 KB |
testcase_10 | AC | 1,233 ms
120,644 KB |
testcase_11 | AC | 1,256 ms
120,560 KB |
testcase_12 | AC | 1,286 ms
120,628 KB |
testcase_13 | AC | 1,247 ms
120,652 KB |
testcase_14 | AC | 1,240 ms
120,560 KB |
testcase_15 | AC | 1,249 ms
120,476 KB |
testcase_16 | AC | 1,244 ms
120,472 KB |
testcase_17 | AC | 1,267 ms
120,584 KB |
testcase_18 | AC | 1,341 ms
123,320 KB |
testcase_19 | AC | 1,309 ms
122,268 KB |
testcase_20 | AC | 1,267 ms
120,868 KB |
testcase_21 | AC | 1,247 ms
120,756 KB |
testcase_22 | AC | 1,243 ms
120,716 KB |
testcase_23 | AC | 1,362 ms
123,832 KB |
testcase_24 | AC | 1,253 ms
120,860 KB |
testcase_25 | AC | 1,322 ms
124,108 KB |
testcase_26 | AC | 1,317 ms
124,496 KB |
testcase_27 | AC | 1,320 ms
123,060 KB |
testcase_28 | AC | 2 ms
6,820 KB |
testcase_29 | AC | 1,247 ms
121,052 KB |
testcase_30 | AC | 1,268 ms
121,760 KB |
testcase_31 | AC | 1 ms
6,820 KB |
testcase_32 | AC | 1 ms
6,816 KB |
testcase_33 | AC | 1,290 ms
122,692 KB |
testcase_34 | AC | 2 ms
6,820 KB |
testcase_35 | AC | 1 ms
6,820 KB |
testcase_36 | AC | 1 ms
6,816 KB |
testcase_37 | AC | 1,290 ms
122,132 KB |
testcase_38 | AC | 1,258 ms
121,740 KB |
testcase_39 | AC | 1,241 ms
120,832 KB |
testcase_40 | AC | 1,282 ms
122,548 KB |
testcase_41 | AC | 1,267 ms
122,772 KB |
testcase_42 | AC | 1,284 ms
122,656 KB |
testcase_43 | AC | 1,277 ms
121,280 KB |
testcase_44 | AC | 1,279 ms
121,276 KB |
testcase_45 | AC | 1,262 ms
120,524 KB |
testcase_46 | AC | 1,261 ms
120,880 KB |
testcase_47 | AC | 1,273 ms
120,472 KB |
testcase_48 | AC | 1,269 ms
120,568 KB |
testcase_49 | AC | 1,296 ms
121,932 KB |
testcase_50 | AC | 1,250 ms
121,120 KB |
testcase_51 | AC | 2 ms
6,816 KB |
testcase_52 | AC | 1,258 ms
121,536 KB |
testcase_53 | AC | 1,268 ms
121,932 KB |
testcase_54 | AC | 1,268 ms
121,124 KB |
testcase_55 | AC | 1,262 ms
121,628 KB |
testcase_56 | AC | 1,251 ms
120,880 KB |
testcase_57 | AC | 2 ms
6,820 KB |
testcase_58 | AC | 1,420 ms
129,420 KB |
testcase_59 | AC | 1,302 ms
123,408 KB |
testcase_60 | AC | 1,332 ms
124,360 KB |
testcase_61 | AC | 1,318 ms
124,020 KB |
testcase_62 | AC | 1,330 ms
124,184 KB |
testcase_63 | AC | 1,314 ms
123,900 KB |
testcase_64 | AC | 1,331 ms
124,308 KB |
testcase_65 | AC | 1,322 ms
123,720 KB |
testcase_66 | AC | 1,404 ms
127,460 KB |
testcase_67 | AC | 1,313 ms
124,344 KB |
testcase_68 | AC | 1,246 ms
120,596 KB |
testcase_69 | AC | 2 ms
6,816 KB |
testcase_70 | AC | 1,247 ms
120,500 KB |
ソースコード
#include<bits/stdc++.h> #include<atcoder/all> using namespace std; using namespace atcoder; typedef long long int ll; typedef long double ld; typedef vector<ll> vi; typedef vector<vi> vvi; typedef vector<vvi> vvvi; typedef vector<vvvi> vvvvi; typedef vector<bool> vb; typedef vector<vb> vvb; typedef vector<vvb> vvvb; typedef vector<vvvb> vvvvb; typedef pair<ll,ll> pi; typedef pair<ll,pi> ppi; typedef pair<ll,ppi> pppi; typedef pair<ll,pppi> ppppi; #define FOR(i,l,r) for(ll i=l;i<r;i++) #define REP(i,n) FOR(i,0,n) #define RFOR(i,l,r) for(ll i=r-1;i>=l;i--) #define RREP(i,n) RFOR(i,0,n) #define ALL(x) x.begin(),x.end() #define F first #define S second #define BS(A,x) binary_search(ALL(A),x) #define LB(A,x) (ll)(lower_bound(ALL(A),x)-A.begin()) #define UB(A,x) (ll)(upper_bound(ALL(A),x)-A.begin()) #define COU(A,x) (UB(A,x)-LB(A,x)) #define sz(c) ((ll)(c).size()) /* #include<boost/multiprecision/cpp_int.hpp> namespace mp=boost::multiprecision; using Bint=mp::cpp_int; */ template<typename T>using min_priority_queue=priority_queue<T,vector<T>,greater<T>>; template<typename T1,typename T2>ostream&operator<<(ostream&os,pair<T1,T2>&p){os<<p.F<<" "<<p.S;return os;} template<typename T1,typename T2>istream&operator>>(istream&is,pair<T1,T2>&p){is>>p.F>>p.S;return is;} template<typename T>ostream&operator<<(ostream&os,vector<T>&v){REP(i,sz(v))os<<v[i]<<(i+1!=sz(v)?" ":"");return os;} template<typename T>istream&operator>>(istream&is,vector<T>&v){for(T&in:v)is>>in;return is;} template<class T>bool chmax(T&a,T b){if(a<b){a=b;return 1;}return 0;} template<class T>bool chmin(T&a,T b){if(b<a){a=b;return 1;}return 0;} ld dist(ld x1,ld y1,ld x2,ld y2){return sqrtl((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));} vi fast_mod_convolution(vi&a,vi&b,ll mod){ const ll m1=167772161,m2=469762049,m3=1224736769; const ll m1_inv_m2=inv_mod(m1,m2); const ll m12_inv_m3=inv_mod(m1*m2,m3); const ll m12_mod=m1*m2%mod; auto x=convolution<m1>(a,b); auto y=convolution<m2>(a,b); auto z=convolution<m3>(a,b); vector<ll>ret(sz(a)+sz(b)-1); REP(i,sz(ret)){ ll v1=(y[i]-x[i])*m1_inv_m2%m2;if(v1<0)v1+=m2; ll v2=(z[i]-(x[i]+m1*v1)%m3)*m12_inv_m3%m3;if(v2<0)v2+=m3; ret[i]=(x[i]+m1*v1+m12_mod*v2)%mod; } return ret; } const ld EPS=1e-8; //* using mint=modint998244353; const ll mod=998244353; //*/ /* using mint=modint1000000007; const ll mod=1000000007; //*/ //using mint=modint; //* typedef vector<mint> vm; typedef vector<vm> vvm; typedef vector<vvm> vvvm; typedef vector<vvvm> vvvvm; ostream&operator<<(ostream&os,mint&a){os<<a.val();return os;} //*/ int main(){ ll N,M;cin>>N>>M; if(N%2){cout<<0<<endl;return 0;} vi D; REP(i,M){ ll L,R;cin>>L>>R; if((R-L)%2)D.emplace_back(R-L+1); } M=sz(D); sort(ALL(D)); vm K(1e7,1),R(1e7),C(1e7); FOR(i,1,1e7)K[i]=i*K[i-1]; REP(i,1e7)R[i]=1/K[i]; REP(i,1e6)C[i]=K[2*i]*R[i]*R[i+1]; sort(ALL(D)); vvm E(M+1); REP(i,M){ E[i].resize(D[i]/2+1); E[i][0]=1; E[i][D[i]/2]=-C[D[i]/2]; } E[M]=vm{1}; while(sz(E)>1){ vvm F; REP(i,sz(E)/2)F.emplace_back(convolution(E[2*i],E[2*i+1])); if(sz(E)%2)F.emplace_back(E.back()); E=F; } mint ans=0; REP(i,sz(E[0]))ans+=C[N/2-i]*E[0][i]; cout<<ans<<endl; return 0; }