NumberTheory/Aizu/aizu_alds1_1_c_check_prime.test.cpp

View this file on GitHub
Last update: 2026-07-26 19:50:18+07:00
Problem: https://onlinejudge.u-aizu.ac.jp/courses/lesson/1/ALDS1/1/ALDS1_1_C

Depends on

Code

#define PROBLEM "https://onlinejudge.u-aizu.ac.jp/courses/lesson/1/ALDS1/1/ALDS1_1_C"

#include "../../template.h"

#include "../Math/CheckPrime.h"


void solve() {
    int n; cin >> n;
    int cnt = 0;
    while (n--) {
        ll x; cin >> x;
        cnt += BruteForce(x);
    }
    cout << cnt << '\n';
}

#line 1 "NumberTheory/Aizu/aizu_alds1_1_c_check_prime.test.cpp"
#define PROBLEM "https://onlinejudge.u-aizu.ac.jp/courses/lesson/1/ALDS1/1/ALDS1_1_C"

#line 2 "template.h"

#include <bits/stdc++.h>

using namespace std;
 
#define ll long long
#define MOD (ll)(1e9+7)
#define all(x) (x).begin(),(x).end()
#define unique(x) x.erase(unique(all(x)), x.end())
#define INF32 ((1ull<<31)-1)
#define INF64 ((1ull<<63)-1)
#define inf (ll)1e18

#define vi vector<int>
#define pii pair<int, int>
#define pll pair<ll, ll>
#define fi first
#define se second

mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
ll get_rand(ll r) { return uniform_int_distribution<ll>(0, r - 1)(rng); }

const int mod = 998244353;

void solve();

int main(){
    ios_base::sync_with_stdio(false);cin.tie(NULL);
    // cin.exceptions(cin.failbit);

    // int t; cin >> t;

    // while(t--)

        solve();
    cerr << "\nTime run: " << 1000 * clock() / CLOCKS_PER_SEC << "ms" << '\n';
    return 0;
}
#line 2 "NumberTheory/Math/CheckPrime.h"


bool BruteForce(ll n) {
    if (n == 2 || n == 3) return true;
    if (n <= 1 || n % 2 == 0 || n % 3 == 0) return false;
    ll lim = sqrtl(n);
    while (lim > 0 && lim > n / lim) lim--;
    for (ll i = 5; i <= lim; i += 6) {
        if (n % i == 0 || n % (i+2) == 0) return false;
    }
    return true;
}

// https://codeforces.com/blog/entry/91632

ll Meissel(ll n) {
    vector<ll> v;
    for (ll i = 1; i*i <= n; i++) {
        v.push_back(i);
        v.push_back(n / i);
    }
    sort(all(v));
    unique(v);
    ll sq = sqrt(n);
    auto geti = [&](ll x) {
        if (x <= sq) return x-1;
        return (int)v.size() - n / x;
    };
    vector<ll> dp = v;
    ll a = 0;
    for (ll p = 2; p*p <= n; p++) {
        if (dp[geti(p)] == dp[geti(p-1)]) continue;
        a++;
        for (int i = (int)v.size()-1; i >= 0; i--) {
            if (v[i] < p * p) break;
            dp[i] -= dp[geti(v[i] / p)] - a;
        }
    }
    return dp[geti(n)] - 1;
}
#line 5 "NumberTheory/Aizu/aizu_alds1_1_c_check_prime.test.cpp"

void solve() {
    int n; cin >> n;
    int cnt = 0;
    while (n--) {
        ll x; cin >> x;
        cnt += BruteForce(x);
    }
    cout << cnt << '\n';
}