結果

問題 No.843 Triple Primes
ユーザー rokahikou1rokahikou1
提出日時 2021-06-30 14:21:14
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 17 ms / 2,000 ms
コード長 11,300 bytes
コンパイル時間 1,894 ms
コンパイル使用メモリ 186,712 KB
実行使用メモリ 5,492 KB
最終ジャッジ日時 2023-09-09 07:40:29
合計ジャッジ時間 4,524 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 16 ms
5,308 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,384 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 13 ms
5,032 KB
testcase_08 AC 14 ms
5,192 KB
testcase_09 AC 16 ms
5,364 KB
testcase_10 AC 14 ms
5,012 KB
testcase_11 AC 15 ms
5,272 KB
testcase_12 AC 17 ms
5,268 KB
testcase_13 AC 16 ms
5,360 KB
testcase_14 AC 16 ms
5,320 KB
testcase_15 AC 13 ms
5,008 KB
testcase_16 AC 13 ms
4,988 KB
testcase_17 AC 2 ms
4,380 KB
testcase_18 AC 1 ms
4,380 KB
testcase_19 AC 2 ms
4,380 KB
testcase_20 AC 4 ms
4,380 KB
testcase_21 AC 3 ms
4,376 KB
testcase_22 AC 9 ms
4,384 KB
testcase_23 AC 9 ms
4,432 KB
testcase_24 AC 5 ms
4,380 KB
testcase_25 AC 4 ms
4,376 KB
testcase_26 AC 17 ms
5,356 KB
testcase_27 AC 2 ms
4,376 KB
testcase_28 AC 16 ms
5,360 KB
testcase_29 AC 5 ms
4,376 KB
testcase_30 AC 16 ms
5,344 KB
testcase_31 AC 3 ms
4,376 KB
testcase_32 AC 2 ms
4,380 KB
testcase_33 AC 5 ms
4,376 KB
testcase_34 AC 8 ms
4,380 KB
testcase_35 AC 16 ms
5,420 KB
testcase_36 AC 3 ms
4,376 KB
testcase_37 AC 12 ms
4,812 KB
testcase_38 AC 10 ms
4,652 KB
testcase_39 AC 16 ms
5,492 KB
testcase_40 AC 2 ms
4,376 KB
testcase_41 AC 1 ms
4,376 KB
testcase_42 AC 14 ms
5,080 KB
testcase_43 AC 15 ms
5,140 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: 関数 ‘std::pair<long long int, long long int> math::crt2(const std::pair<long long int, long long int>&, const std::pair<long long int, long long int>&)’ 内:
main.cpp:256:10: 警告: structured bindings only available with ‘-std=c++17’ or ‘-std=gnu++17’ [-Wc++17-extensions]
  256 |     auto [r1, m1] = rm1;
      |          ^
main.cpp:257:10: 警告: structured bindings only available with ‘-std=c++17’ or ‘-std=gnu++17’ [-Wc++17-extensions]
  257 |     auto [r2, m2] = rm2;
      |          ^
main.cpp: 関数 ‘long long int math::garner(const std::vector<std::pair<long long int, long long int> >&)’ 内:
main.cpp:279:14: 警告: structured bindings only available with ‘-std=c++17’ or ‘-std=gnu++17’ [-Wc++17-extensions]
  279 |         auto [r, m] = rm[i];
      |              ^
main.cpp: 関数 ‘long long int math::garner_mod(std::vector<std::pair<long long int, long long int> >, long long int)’ 内:
main.cpp:300:14: 警告: structured bindings only available with ‘-std=c++17’ or ‘-std=gnu++17’ [-Wc++17-extensions]
  300 |         auto [r, m] = rm[i];
      |              ^

ソースコード

diff # 

#pragma region Macros
#include <bits/stdc++.h>
#define rep(i, n) for(int(i) = 0; (i) < (n); (i)++)
#define rrep(i, n) for(int(i) = (n)-1; (i) >= 0; (i)--)
#define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++)
#define ROF(i, m, n) for(int(i) = (n)-1; (i) >= (m); (i)--)
#define ALL(v) (v).begin(), (v).end()
#define LLA(v) (v).rbegin(), (v).rend()
#define SZ(v) (int)(v).size()
#define INT(...)     \
    int __VA_ARGS__; \
    read(__VA_ARGS__)
#define LL(...)     \
    ll __VA_ARGS__; \
    read(__VA_ARGS__)
#define DOUBLE(...)     \
    double __VA_ARGS__; \
    read(__VA_ARGS__)
#define CHAR(...)     \
    char __VA_ARGS__; \
    read(__VA_ARGS__)
#define STRING(...)     \
    string __VA_ARGS__; \
    read(__VA_ARGS__)
#define VEC(type, name, size) \
    vector<type> name(size);  \
    read(name)
#define VEC2(type, name, height, width)                     \
    vector<vector<type>> name(height, vector<type>(width)); \
    read(name)
#define DVEC(type, name1, name2, size)     \
    vector<type> name1(size), name2(size); \
    read(name1, name2);
#define TVEC(type, name1, name2, name3, size)           \
    vector<type> name1(size), name2(size), name3(size); \
    read(name1, name2, name3);
using namespace std;
using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
const int INF = 1 << 30;
const ll LINF = 1LL << 60;
const int MOD = 1e9 + 7;
const char newl = '\n';
const int dx[] = {1, 0, -1, 0};
const int dy[] = {0, 1, 0, -1};
template <class T> inline bool between(T x, T l, T r) { return l <= x && x < r; }
template <class T> inline vector<T> make_vec(size_t a, T val) { return vector<T>(a, val); }
template <class... Ts> inline auto make_vec(size_t a, Ts... ts) { return vector<decltype(make_vec(ts...))>(a, make_vec(ts...)); }
void read() {}
template <class T> inline void read(T &a) { cin >> a; }
template <class T, class S> inline void read(pair<T, S> &p) { read(p.first), read(p.second); }
template <class T> inline void read(vector<T> &v) {
    for(auto &&a : v) read(a);
}
template <class T, class U> inline void read(vector<T> &a, vector<U> &b) {
    for(int i = 0; i < a.size(); i++) {
        read(a[i]);
        read(b[i]);
    }
}
template <class T, class U, class V> inline void read(vector<T> &a, vector<U> &b, vector<V> &c) {
    for(int i = 0; i < a.size(); i++) {
        read(a[i]);
        read(b[i]);
        read(c[i]);
    }
}
template <class Head, class... Tail> inline void read(Head &head, Tail &...tail) { read(head), read(tail...); }
template <class T> void write(const T &a) { cout << a << '\n'; }
template <class T> void write(const vector<T> &a) {
    for(int i = 0; i < a.size(); i++) cout << a[i] << (i + 1 == a.size() ? '\n' : ' ');
}
template <class Head, class... Tail> void write(const Head &head, const Tail &...tail) {
    cout << head << ' ';
    write(tail...);
}
template <class T> void writel(const T &a) { cout << a << '\n'; }
template <class T> void writel(const vector<T> &a) {
    for(int i = 0; i < a.size(); i++) cout << a[i] << '\n';
}
template <class Head, class... Tail> void writel(const Head &head, const Tail &...tail) {
    cout << head << '\n';
    write(tail...);
}
template <class T> auto sum(const vector<T> &a) { return accumulate(ALL(a), T(0)); }
template <class T> auto min(const vector<T> &a) { return *min_element(ALL(a)); }
template <class T> auto max(const vector<T> &a) { return *max_element(ALL(a)); }
template <class T, class U> void msort(vector<T> &a, vector<U> &b) {
    assert(a.size() == b.size());
    vector<pair<T, U>> ab(a.size());
    for(int i = 0; i < a.size(); i++) ab[i] = {a[i], b[i]};
    sort(ALL(ab));
    for(int i = 0; i < a.size(); i++) {
        a[i] = ab[i].first;
        b[i] = ab[i].second;
    }
}
template <class T, class U, class V> void msort(vector<T> &a, vector<U> &b, vector<V> &c) {
    assert(a.size() == b.size() && b.size() == c.size());
    vector<tuple<T, U, V>> abc(a.size());
    for(int i = 0; i < a.size(); i++) abc[i] = {a[i], b[i], c[i]};
    sort(ALL(abc));
    for(int i = 0; i < a.size(); i++) {
        a[i] = get<0>(abc[i]);
        b[i] = get<1>(abc[i]);
        c[i] = get<2>(abc[i]);
    }
}
template <class T, class U> inline bool chmax(T &a, U b) {
    if(a < b) {
        a = b;
        return 1;
    }
    return 0;
}
template <class T, class U> inline bool chmin(T &a, U b) {
    if(a > b) {
        a = b;
        return 1;
    }
    return 0;
}
inline int bsf(int v) { return __builtin_ctz(v); } // 最下位の1が下から何番目か
inline int bsf(ll v) { return __builtin_ctzll(v); }
inline int bsr(int v) { return 31 - __builtin_clz(v); } // 最上位の1が下から何番目か
inline int bsr(ll v) { return 63 - __builtin_clzll(v); }
inline int lsb(int v) { return v & -v; } // 最上位の1だけ残す
inline ll lsb(ll v) { return v & -v; }
inline int msb(int v) { return 1 << bsr(v); } // 最上位の1だけ残す
inline ll msb(ll v) { return 1LL << bsr(v); }
struct IO {
    IO() {
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
        cout << fixed << setprecision(10);
    }
} io;
#pragma endregion

namespace math {

std::set<long long> divisor(long long n) {
    std::set<long long> ret;
    for(long long i = 1; i * i <= n; i++) {
        if(n % i == 0) {
            ret.insert(i);
            if(i * i != n)
                ret.insert(n / i);
        }
    }
    return ret;
}

// 素数篩＋素因数分解
// 初期化O(N),素因数分解O(logN)
// たくさん素因数分解するときはこっち
struct Sieve {
    int N;
    std::vector<int> sieve;
    Sieve(int n) : N(n + 1), sieve(n + 1) { init(); }
    void init() {
        std::iota(sieve.begin(), sieve.end(), 0);
        for(int i = 2; i * i <= N; i++) {
            if(sieve[i] < i)
                continue;
            for(int j = i * i; j <= N; j += i) {
                if(sieve[j] == j)
                    sieve[j] = i;
            }
        }
    }

    bool is_prime(int x) {
        assert(x <= N);
        if(x == 1) return false;
        return sieve[x] == x;
    }

    std::map<long long, int> prime_factorize(long long n) {
        assert(n <= N);
        std::map<long long, int> ret;
        while(n > 1) {
            ret[sieve[n]]++;
            n = n / sieve[n];
        }
        return ret;
    }
};

// 素因数分解
// O(sqrt(N))
std::map<long long, int> prime_factor(long long n) {
    std::map<long long, int> ret;
    for(long long i = 2; i * i <= n; i++) {
        while(n % i == 0) {
            ret[i]++;
            n /= i;
        }
    }
    if(n != 1)
        ret[n] = 1;
    return ret;
}

long long mod_pow(long long x, long long n, long long mod) {
    if(n == 0)
        return 1;
    long long res = mod_pow(x * x % mod, n / 2, mod);
    if(n & 1)
        res = res * x % mod;
    return res;
}

long long euler_phi(long long n) {
    long long ret = n;
    for(long long i = 2; i * i <= n; i++) {
        if(n % i == 0) {
            ret -= ret / i;
            while(n % i == 0)
                n /= i;
        }
    }
    if(n > 1)
        ret -= ret / n;
    return ret;
}

long long extgcd(long long a, long long b, long long &x, long long &y) {
    if(b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    long long d = extgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

long long mod_inv(long long a, long long mod) {
    long long b = mod, u = 1, v = 0;
    while(b) {
        long long t = a / b;
        std::swap(a -= t * b, b);
        std::swap(u -= t * b, v);
    }
    u %= mod;
    if(u < 0)
        u += mod;
    return u;
}

std::pair<long long, long long> crt2(const std::pair<long long, long long> &rm1,
                                     const std::pair<long long, long long> &rm2) {
    long long p, q;
    auto [r1, m1] = rm1;
    auto [r2, m2] = rm2;
    long long d = extgcd(r1, m2, p, q);
    if((r2 - r1) % d != 0)
        return {0, -1};
    long long m = m1 * (m2 / d);
    long long tmp = (r2 - r1) / d * p % (m2 / d);
    long long c = (r1 + m1 * tmp);
    long long r = (c % m + m) % m;
    return {r, m};
}

/**
 * @brief  Garnerのアルゴリズム. x % m[i] == b[i] % m[i]を満たす最小のxを求める.
 *
 * @param rm {あまり,法}の配列. !!法はすべて互いに素!!
 * @return long long x % m[i] == b[i] % m[i]を満たす最小のx.
 */
long long garner(const std::vector<std::pair<long long, long long>> &rm) {
    int n = rm.size();
    long long m_prod = 1;
    long long res = rm[0].first % rm[0].second;
    for(int i = 1; i < n; i++) {
        auto [r, m] = rm[i];
        m_prod *= rm[i - 1].second;
        long long t = (((r - res) * mod_inv(m_prod, m)) % m + m) % m;
        res += t * m_prod;
    }
    return res;
}

/**
 * @brief  Garnerのアルゴリズム. x % m[i] == b[i] % m[i]を満たす最小のxを求める.
 *
 * @param rm {あまり,法}の配列. !!法はすべて互いに素!!
 * @param mod 答えの法.
 * @return long long x % m[i] == b[i] % m[i]を満たす最小のx (mod M).
 */
long long garner_mod(std::vector<std::pair<long long, long long>> rm,
                     long long mod) {
    rm.emplace_back(0, mod);
    int n = rm.size();
    std::vector<long long> coffs(n, 1), constants(n, 0);
    for(int i = 0; i < n - 1; i++) {
        auto [r, m] = rm[i];
        long long t = ((r - constants[i]) * mod_inv(coffs[i], m) % m + m) & m;
        for(int j = i + 1; j < n; j++) {
            (constants[j] += coffs[j] * t) %= rm[j].second;
            (coffs[j] *= m) %= rm[j].second;
        }
    }
    return constants.back();
}

long long floor_div(long long a, long long b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? a / b : (a - (b - 1)) / b;
}

long long ceil_div(long long a, long long b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? (a + b - 1) / b : a / b;
}

int bit_length(long long x) {
    assert(x >= 0);
    int ret = 0;
    while(x > 0) {
        ret++;
        x /= 2;
    }
    return ret;
}

long long ll_sqrt(long long x) {
    long long ret = sqrt(x) - 1;
    while((ret + 1) * (ret + 1) <= x)
        ret++;
    return ret;
}

// multiple of d in [left,right],left>0
long long cmul(long long left, long long right, long long d) {
    if(left == 0)
        return right / d;
    return right / d - (left - 1) / d;
}

// multiple of d in [left,right]
long long count_multiple(long long left, long long right, long long d) {
    if(right < 0)
        return cmul(-right, -left, d);
    if(left > 0)
        return cmul(left, right, d);
    return cmul(0, right, d) + cmul(0, -left, d) - 1;
}

}; // namespace math

using namespace math;

void solve() {
    INT(n);
    Sieve sieve(n);
    vector<int> primes;
    FOR(i, 2, n + 1) {
        if(sieve.is_prime(i)) primes.push_back(i);
    }
    int m = primes.size();
    int res = 0;
    for(int i = 0; i < m; i++) {
        ll r2 = primes[i] * primes[i];
        if(r2 > n * 2) break;
        for(int j = 0; j < m; j++) {
            if(primes[i] * primes[i] - primes[j] > n || primes[i] * primes[i] - primes[j] < 0) continue;
            if(sieve.is_prime(primes[i] * primes[i] - primes[j])) res++;
        }
    }
    write(res);
}

int main() {
    // INT(t);
    int t = 1;
    while(t--) {
        solve();
    }
}
0