結果
問題 | No.886 Direct |
ユーザー |
![]() |
提出日時 | 2019-08-20 17:52:15 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.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 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 32 |
ソースコード
#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; }