結果
| 問題 | No.886 Direct |
| ユーザー |
 pockyny
|
| 提出日時 | 2019-09-13 23:06:40 |
| 言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 915 bytes |
| コンパイル時間 | 730 ms |
| コンパイル使用メモリ | 74,108 KB |
| 実行使用メモリ | 14,976 KB |
| 最終ジャッジ日時 | 2024-07-04 10:21:44 |
| 合計ジャッジ時間 | 11,246 ms |
|
ジャッジサーバーID （参考情報） |
judge1 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | TLE | - |
| testcase_01 | -- | - |
| testcase_02 | -- | - |
| testcase_03 | -- | - |
| testcase_04 | -- | - |
| testcase_05 | -- | - |
| testcase_06 | -- | - |
| testcase_07 | -- | - |
| testcase_08 | -- | - |
| testcase_09 | -- | - |
| testcase_10 | -- | - |
| testcase_11 | -- | - |
| testcase_12 | -- | - |
| testcase_13 | -- | - |
| testcase_14 | -- | - |
| testcase_15 | -- | - |
| testcase_16 | -- | - |
| testcase_17 | -- | - |
| testcase_18 | -- | - |
| testcase_19 | -- | - |
| testcase_20 | -- | - |
| testcase_21 | -- | - |
| testcase_22 | -- | - |
| testcase_23 | -- | - |
| testcase_24 | -- | - |
| testcase_25 | -- | - |
| testcase_26 | -- | - |
| testcase_27 | -- | - |
| testcase_28 | -- | - |
| testcase_29 | -- | - |
| testcase_30 | -- | - |
| testcase_31 | -- | - |
| testcase_32 | -- | - |
| testcase_33 | -- | - |
| testcase_34 | -- | - |
| testcase_35 | -- | - |
コンパイルメッセージ
main.cpp: In function 'll solve(int)':
main.cpp:34:1: warning: no return statement in function returning non-void [-Wreturn-type]
34 | }
| ^
ソースコード
#include <iostream>
#include <vector>
using namespace std;
typedef long long ll;
int used[300010];
ll ans = 0,mod = 1000000007,h,w;
ll solve(int x){
int y = x;
vector<int> fac;
while(y!=1){
int k = used[y];
fac.push_back(k);
while(y%k==0 && y!=1){
y /= k;
}
}
int i,j,s = fac.size();
ll num = 0,sum = 0;
for(i=0;i<(1<<s);i++){
ll z = 1,sig = 1;
for(j=0;j<s;j++){
if((1<<j) & i){
z *= fac[j];
sig *= -1;
}
}
ll y = w/z;
num += sig*y;
sum += sig*z*y*(y + 1)/2;
}
/*cout << x << endl;
cout << num << " " << sum << endl;*/
(ans += (h - x)*(w*num - sum + mod)%mod) %= mod;
}
int main(){
int i,j;
cin >> h >> w;
for(i=1;i<=300000;i++){
used[i] = -1;
}
used[1] = 1;
for(i=2;i<=h;i++){
for(j=i;j<=h;j+=i){
if(used[j]!=-1) break;;
used[j] = i;
}
}
for(i=1;i<h;i++){
solve(i);
}
ans *= 2;
(ans += h*(w - 1) + w*(h - 1)) %= mod;
cout << ans << endl;
}