結果
問題 |
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; }