結果
問題 | No.2151 3 on Torus-Lohkous |
ユーザー | hotman78 |
提出日時 | 2022-12-15 09:30:31 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 877 bytes |
コンパイル時間 | 3,269 ms |
コンパイル使用メモリ | 229,476 KB |
実行使用メモリ | 11,188 KB |
最終ジャッジ日時 | 2024-04-26 11:56:13 |
合計ジャッジ時間 | 4,531 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | AC | 140 ms
10,980 KB |
testcase_04 | AC | 120 ms
11,008 KB |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | AC | 138 ms
11,184 KB |
ソースコード
#include<bits/stdc++.h> #include<atcoder/modint> #include<atcoder/convolution> using namespace std; using lint = long long; #define rep(i,n) for(lint i=0;i<(n);++i) using mint=atcoder::static_modint<998244353>; int main(){ vector<mint>table(1000000,1),inv(1000000); rep(i,1000000)table[i+1]=table[i]*(i+1); inv[999999]=table[999999].inv(); for(lint i=999998;i>=0;--i)inv[i]=inv[i+1]*(i+1); auto comb=[&](auto a,auto b){ return table[a]*inv[b]*inv[a-b]; }; lint t; cin>>t; while(t--){ lint h,w; cin>>h>>w; lint n=gcd(h,w); mint ans=mint(h)*w; ans+=2*n*(h/n)*(w/n); if(h*w/n%3==0){ if(n>4)ans+=n*12; } if(h%4==0&&w%4==0){ ans+=16; } if(n%2==0){ rep(p,2){ mint x=0; for(lint i=(n+3*p+5)%6+1;i<=n;i+=6){ if(0<=n/2-5*i){ x+=comb((n/2-5*i)/3+i-1,i-1)/i; } } ans+=x*x*n*n; } } cout<<ans.val()<<endl; } }