結果

問題 No.917 Make One With GCD
コンテスト
ユーザー tkmst201
提出日時 2019-10-26 16:35:15
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 11 ms / 2,000 ms
コード長 1,402 bytes
コンパイル時間 1,870 ms
コンパイル使用メモリ 185,516 KB
実行使用メモリ 11,880 KB
最終ジャッジ日時 2024-06-24 11:09:38
合計ジャッジ時間 3,021 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define FOR(i,a,b) for(int i=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
#define fi first
#define se second
template<typename A, typename B> inline bool chmax(A &a, B b) { if (a<b) { a=b; return 1; } return 0; }
template<typename A, typename B> inline bool chmin(A &a, B b) { if (a>b) { a=b; return 1; } return 0; }
using ll = long long;
using pii = pair<int, int>;
constexpr ll INF = 1ll<<30;
constexpr ll longINF = 1ll<<60;
constexpr ll MOD = 1000000007;
constexpr bool debug = 0;
//---------------------------------//

vector<ll> findDivisor(ll n) {
	vector<ll> divisor;
	for (ll i = 1; i * i <= n; i++) {
		if (n % i == 0) {
			divisor.push_back(i);
			if (i * i != n) divisor.push_back(n / i);
		}
	}
	sort(ALL(divisor));
	return divisor;
}

ll gcd(ll a, ll b) {
	if (b == 0) return a;
	return gcd(b, a % b);
}


int N;
int A[50];
ll dp[51][25601];

int main() {
	cin >> N;
	REP(i, N) cin >> A[i];
	map<int, int> mm;
	vector<int> rev;
	REP(i, N) {
		vector<ll> v = findDivisor(A[i]);
		for (ll cur : v) if (mm.find(cur) == mm.end()) {
			int sz = mm.size();
			mm[cur] = sz;
			rev.push_back(cur);
		}
	}
	
	REP(i, N) {
		REP(j, rev.size()) if (dp[i][j]) {
			dp[i + 1][mm[gcd(A[i], rev[j])]] += dp[i][j];
			dp[i + 1][j] += dp[i][j];
		}
		dp[i + 1][mm[A[i]]]++;
	}
	
	cout << dp[N][mm[1]] << endl;
	
	return 0;
}
0