結果
問題 | No.917 Make One With GCD |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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; }