結果

問題 No.983 Convolution
ユーザー koi_kotya
提出日時 2020-02-11 14:47:39
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,850 bytes
コンパイル時間 1,639 ms
コンパイル使用メモリ 183,344 KB
実行使用メモリ 18,944 KB
最終ジャッジ日時 2024-10-01 07:45:20
合計ジャッジ時間 7,355 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18 WA * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int,int> P;

#define p_ary(ary,a,b) do { cout << "["; for (int count = (a);count < (b);++count) cout << ary[count] << ((b)-1 == count ? "" : ", "); cout << "]\n"; } while(0)
#define p_map(map,it) do {cout << "{";for (auto (it) = map.begin();;++(it)) {if ((it) == map.end()) {cout << "}\n";break;}else cout << "" << (it)->first << "=>" << (it)->second << ", ";}}while(0)

ll gcd(ll a,ll b) {
    while (b) {
        a %= b;
        swap(a,b);
    }
    return a;
}


int const MAX_N = 1000;
vector<ll> prime;
vector<bool> is_prime(MAX_N,true);
void Eratosthenes() {
    is_prime[0] = is_prime[1] = false;
    for (ll i = 2;i*i < MAX_N;++i) if (is_prime[i]) for (ll j = 2*i;j < MAX_N;j += i) is_prime[j] = false;
    for (ll i = 0;i < MAX_N;++i) if (is_prime[i]) prime.push_back(i);
}

int main() {
    Eratosthenes();
    int n;cin >> n;
    vector<ll> a(1000000),c;
    for (int i = 0;i < n;++i) cin >> a[i];
    c = a;
    for (int i = 0;i < 19;++i) for (int j = 0;j < n;++j) if (!(j&(1<<i)) && c[j] > -1) {
        if (c[i|(1<<i)] > -1) c[j] += c[i|(1<<i)];
        else c[j] = -1;
    }
    map<ll,int> mp;
    for (int i = 0;i < n;++i) if (a[i] > -1) {
        map<ll,int> mp2;
        mp2[1]++;
        for (ll j : prime) {
            if (j*j > a[i]) break;
            while (a[i]%j == 0) mp2[j]++,a[i] /= j;
        }
        if (a[i] > 1) mp2[a[i]]++;
        for (ll j : prime) {
            if (j*j > c[i]) break;
            while (c[i]%j == 0) mp2[j]++,c[i] /= j;
        }
        if (c[i] > 1) mp2[c[i]]++;
        if (mp.size() == 0) mp = mp2;
        else for (auto& p : mp) p.second = min(p.second,mp2[p.first]);
    }
    ll ans = 1;
    for (auto p : mp) for (int i = 0;i < p.second;++i) ans *= p.first;
    cout << (mp.size() ? ans : -1) << endl;
}
0