結果

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

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
using ll     = long long;
using pii    = pair<int, int>;
using pll    = pair<ll, ll>;
using vi     = vector<int>;
using vl     = vector<ll>;
using vvi    = vector<vi>;
using vvl    = vector<vl>;
const ll INF = 1LL << 60;
const ll MOD = 1000000007;
template <class T>
bool chmax(T &a, const T &b) {
    return (a < b) ? (a = b, 1) : 0;
}
template <class T>
bool chmin(T &a, const T &b) {
    return (b < a) ? (a = b, 1) : 0;
}
template <class C>
void print(const C &c, std::ostream &os = std::cout) {
    std::copy(std::begin(c), std::end(c), std::ostream_iterator<typename C::value_type>(os, " "));
    os << std::endl;
}

template <typename T = int>
void primeFactors(T a, map<T, int> &facs) {
    double sqrtA = sqrt(a);
    for (int i = 2; i <= sqrtA + 1e-10; ++i) {
        while (a % i == 0) {
            facs[i]++;
            a /= i;
        }
    }
    if (a > sqrtA)
        facs[a]++;
    return;
}

// list up all factors
template <typename T>
set<T> factors(T a) {
    set<T> facs;
    for (T i = 1; i * i <= a; ++i) {
        if (a % i == 0) {
            facs.insert(i);
            facs.insert(a / i);
        }
    }
    return facs;
}

int main() {
    int n;
    cin >> n;
    vl a(n);
    for (int i = 0; i < n; ++i) {
        cin >> a[i];
    }
    map<ll, int> facs;
    for (int i = 0; i < n; ++i) {
        auto fac = factors(a[i]);
        for (auto &f : fac) {
            facs[f]++;
        }
    }
    ll ret = 1;
    for (int i = 0; i < n; ++i) {
        ret *= 2;
    }
    ret -= 1; // empty set
    for (auto &p : facs) {
        if (p.first == 1)
            continue;
        ll temp = 1;
        for (int i = 0; i < p.second; ++i) {
            temp *= 2;
        }
        temp--;
        map<ll, int> primefac;
        primeFactors(p.first, primefac);
        for (auto &pp : primefac) {
            if (pp.second > 1) {
                temp = 0;
            }
        }
        temp *= int(primefac.size()) % 2 == 0 ? -1 : 1;
        ret -= temp;
    }
    cout << ret << "\n";
    return 0;
}
0