結果
問題 | No.886 Direct |
ユーザー | tempura_pp |
提出日時 | 2019-08-20 17:52:15 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 533 ms / 4,000 ms |
コード長 | 787 bytes |
コンパイル時間 | 1,509 ms |
コンパイル使用メモリ | 171,340 KB |
実行使用メモリ | 87,512 KB |
最終ジャッジ日時 | 2024-10-06 12:41:34 |
合計ジャッジ時間 | 21,052 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 533 ms
87,388 KB |
testcase_01 | AC | 492 ms
87,384 KB |
testcase_02 | AC | 510 ms
87,384 KB |
testcase_03 | AC | 493 ms
87,388 KB |
testcase_04 | AC | 487 ms
87,380 KB |
testcase_05 | AC | 476 ms
87,388 KB |
testcase_06 | AC | 499 ms
87,384 KB |
testcase_07 | AC | 500 ms
87,380 KB |
testcase_08 | AC | 525 ms
87,368 KB |
testcase_09 | AC | 524 ms
87,384 KB |
testcase_10 | AC | 518 ms
87,380 KB |
testcase_11 | AC | 508 ms
87,316 KB |
testcase_12 | AC | 528 ms
87,388 KB |
testcase_13 | AC | 495 ms
87,384 KB |
testcase_14 | AC | 485 ms
87,388 KB |
testcase_15 | AC | 517 ms
87,384 KB |
testcase_16 | AC | 505 ms
87,380 KB |
testcase_17 | AC | 498 ms
87,380 KB |
testcase_18 | AC | 482 ms
87,380 KB |
testcase_19 | AC | 513 ms
87,384 KB |
testcase_20 | AC | 500 ms
87,380 KB |
testcase_21 | AC | 491 ms
87,384 KB |
testcase_22 | AC | 481 ms
87,384 KB |
testcase_23 | AC | 484 ms
87,384 KB |
testcase_24 | AC | 492 ms
87,380 KB |
testcase_25 | AC | 518 ms
87,384 KB |
testcase_26 | AC | 486 ms
87,384 KB |
testcase_27 | AC | 529 ms
87,384 KB |
testcase_28 | AC | 487 ms
87,384 KB |
testcase_29 | AC | 487 ms
87,380 KB |
testcase_30 | AC | 474 ms
87,512 KB |
testcase_31 | AC | 527 ms
87,384 KB |
testcase_32 | AC | 499 ms
87,384 KB |
testcase_33 | AC | 500 ms
87,512 KB |
testcase_34 | AC | 505 ms
87,384 KB |
testcase_35 | AC | 475 ms
87,384 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll=long long; constexpr ll mod = 1e9+7; int main(){ ll h, w; cin>>h>>w; int n=1<<22; vector<ll> cnt1(n), cnt2(n); for(int i=1;i<h;++i)cnt1[i]=h-i; for(int i=1;i<w;++i)cnt2[i]=w-i; vector<int> prime, isprime(n,1); for(int i=2;i<n;++i){ if(!isprime[i])continue; prime.push_back(i); for(int j=2*i;j<n;j+=i){ isprime[j]=false; } } for(auto p:prime){ int m=(n-1)/p; for(int i=m;i>0;--i){ (cnt1[i]+=cnt1[i*p])%=mod; (cnt2[i]+=cnt2[i*p])%=mod; } } for(int i=1;i<n;++i)(cnt1[i]*=cnt2[i])%=mod; for(auto p:prime){ int m=(n-1)/p; for(int i=1;i<=m;++i){ (cnt1[i]+=mod-cnt1[i*p])%=mod; } } cout<<(2*cnt1[1]+h*(w-1)+w*(h-1))%mod<<endl; }