結果

問題 No.917 Make One With GCD
ユーザー ningenMeningenMe
提出日時 2019-09-09 05:02:55
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,618 bytes
コンパイル時間 1,802 ms
コンパイル使用メモリ 183,404 KB
実行使用メモリ 19,200 KB
最終ジャッジ日時 2024-06-28 00:32:05
合計ジャッジ時間 7,423 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 AC 238 ms
19,200 KB
testcase_04 AC 238 ms
19,200 KB
testcase_05 AC 237 ms
19,200 KB
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 AC 241 ms
19,072 KB
testcase_28 AC 233 ms
19,200 KB
testcase_29 AC 235 ms
19,200 KB
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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));
}

int main() {
	int N,MAX_A = 2000; cin >> N;
	assert(1<=N && N <= 40);
	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] <= MAX_A);
	
	// gcdを前計算 O(MAX_A*MAX_A*log(MAX_A))
	vector<vector<int>> gcd(MAX_A+1,vector<int>(MAX_A+1,0));
	for(int i = 1; i <= MAX_A; ++i) for(int j = i; j <= MAX_A; ++j) gcd[i][j] = gcd[j][i] = GCD(i,j);
	
	//約数のindexをメモするテーブル
	vector<int> table(MAX_A+1,0);

	//左端決め打ち全探索 O(N*N*N*M+N*sqrt(A)+N*M*logM)
	long long ans = 0;
	for(int i = 0; i < N; ++i) {

		//約数列挙 O(sqrt(A))
		auto d = Divisor(A[i]);
		int M = d.size(),idx = 0;
		vector<int> val;

		//添え字前計算 O(M*logM)
		for(auto& e:d) {
			table[e] = idx++;
			val.push_back(e);
		}

		// dp[i][j] := (最後にA[i]を使ってgcdがval[j]になるような場合の数) 
		vector<vector<long long>> dp(N, vector<long long>(M,0));
		dp[i][table[A[i]]] = 1;

		// 遷移 O(N*N*M)
		for(int j = i; j < N; ++j) {
			for(int k = j + 1; k < N; ++k){
				for(int l = 0; l < M; ++l){
					dp[k][table[gcd[val[l]][A[k]]]] += dp[j][l];
				}
			}
			ans += dp[j][table[1]];
		}
	}
	cout << ans << endl;
	return 0;
}
0