結果

問題 No.300 平方数
コンテスト
ユーザー kimiyuki
提出日時 2016-07-05 18:12:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 9 ms / 1,000 ms
コード長 1,220 bytes
コンパイル時間 755 ms
コンパイル使用メモリ 68,424 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-12 20:17:20
合計ジャッジ時間 2,067 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cstdio>
#include <vector>
#include <set>
#include <cmath>
typedef long long ll;
using namespace std;

vector<int> sieve_of_eratosthenes(int n) { // enumerate primes in [2,n] with O(n log log n)
    vector<bool> is_prime(n+1, true);
    is_prime[0] = is_prime[1] = false;
    for (int i = 2; i*i <= n; ++i)
        if (is_prime[i])
            for (int k = i+i; k <= n; k += i)
                is_prime[k] = false;
    vector<int> primes;
    for (int i = 2; i <= n; ++i)
        if (is_prime[i])
            primes.push_back(i);
    return primes;
}
vector<ll> factors(ll n, vector<int> const & primes) {
    vector<ll> result;
    for (int p : primes) {
        if (n < p *(ll) p) break;
        while (n % p == 0) {
            result.push_back(p);
            n /= p;
        }
    }
    if (n != 1) result.push_back(n);
    return result;
}

int main() {
    ll n; scanf("%lld", &n);
    auto primes = sieve_of_eratosthenes(sqrt(n));
    vector<ll> ps = factors(n, primes);
    set<ll> qs;
    for (ll p : ps) {
        if (qs.count(p)) {
            qs.erase(p);
        } else {
            qs.insert(p);
        }
    }
    ll r = 1;
    for (ll q : qs) r *= q;
    printf("%lld\n", r);
    return 0;
}
0