結果
| 問題 |
No.917 Make One With GCD
|
| ユーザー |
|
| 提出日時 | 2019-09-09 05:19:55 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 931 bytes |
| コンパイル時間 | 1,636 ms |
| コンパイル使用メモリ | 177,412 KB |
| 実行使用メモリ | 6,948 KB |
| 最終ジャッジ日時 | 2024-06-28 00:33:18 |
| 合計ジャッジ時間 | 6,493 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 RE * 1 |
| other | AC * 2 WA * 1 RE * 29 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
//divisor O(sqrt(N))
set<long long> Divisor(long long N) {
set<long long> ret;
for (long long i = 1; i*i <= N; ++i) {
if (N%i == 0) {
ret.insert(i);
ret.insert(N / i);
}
}
return ret;
}
//Greatest Common Divisor
long long GCD(long long a, long long b) {
return ((b == 0) ? a : GCD(b, a % b));
}
//GCD随時計算ver. logで落ちそう
int main() {
int N,MAX_A = 2000; cin >> N;
vector<int> A(N);
for(int i = 0; i < N; ++i) cin >> A[i];
// dp[i][j] := (最後にA[i]を使ってgcdがjになるような場合の数)
vector<vector<long long>> dp(N, vector<long long>(MAX_A+1,0));
int ans = 0;
for(int i = 0; i < N; ++i) {
dp[i][A[i]]++;
for(int j = i+1; j < N; ++j) {
for(int k = 1; k <= MAX_A; ++k){
dp[j][GCD(k,A[j])] += dp[i][k];
}
}
ans += dp[i][1];
}
cout << ans << endl;
return 0;
}