結果
| 問題 |
No.917 Make One With GCD
|
| ユーザー |
|
| 提出日時 | 2019-09-26 06:27:04 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 3 ms / 2,000 ms |
| コード長 | 758 bytes |
| コンパイル時間 | 1,641 ms |
| コンパイル使用メモリ | 178,156 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-06-24 11:03:11 |
| 合計ジャッジ時間 | 2,598 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 32 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
//Greatest Common Divisor
long long GCD(long long a, long long b) {
return ((b == 0) ? a : GCD(b, a % b));
}
int main() {
int N; cin >> N;
assert(1<= N && N <= 50);
vector<int> A(N);
for(int i = 0; i < N; ++i) cin >> A[i];
for(int i = 0; i < N; ++i) assert(1 <= A[i] && A[i] <= 100000000);
// dp[i][j] := (最後にA[i]を使ってgcdがjになるような場合の数)
vector<map<int,long long>> dp(N);
long long ans = 0;
//dp O(N*N*d(A)*(log(A)+log(d(A))) d(A) Aの約数の個数
for(int i = 0; i < N; ++i) {
dp[i][A[i]]++;
for(int j = i+1; j < N; ++j) {
for(auto e:dp[i]){
dp[j][GCD(e.first,A[j])] += e.second;
}
}
ans += dp[i][1];
}
cout << ans << endl;
return 0;
}