結果

問題 No.407 鴨等素数間隔列の数え上げ
コンテスト
ユーザー Tatsuno
提出日時 2016-08-05 23:08:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,936 bytes
コンパイル時間 1,839 ms
コンパイル使用メモリ 171,884 KB
実行使用メモリ 13,760 KB
最終ジャッジ日時 2024-11-07 01:07:06
合計ジャッジ時間 5,192 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 15 TLE * 1 -- * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _CRT_SECURE_NO_WARNINGS
#define _USE_MATH_DEFINES

#include <bits/stdc++.h>
using namespace std;
using i32 = int32_t;  using i64 = int64_t;  using str = string;
using u32 = uint32_t; using u64 = uint64_t; using usize = size_t;
using f64 = double_t; template <typename T> using vec = vector<T>;
#define times(n, i) for (i64 i = 0; i < (n); ++i)
#define range(n, m, i) for (i64 i = (n); i < (m); ++i)
#define upto(n, m, i) for (i64 i = (n); i <= (m); ++i)
#define downto(n, m, i) for (i64 i = (n); i >= (m); --i)
#define foreach(xs, x) for (auto &x : (xs))
#define all(xs) (xs).begin(), (xs).end()
#define sortall(xs) sort(all(xs))
#define reverseall(xs) reverse(all(xs))
#define uniqueall(xs) (xs).erase(unique(all(xs)), (xs).end())
#define maximum(xs) (*max_element(all(xs)))
#define minimum(xs) (*min_element(all(xs)))
#define even(x) (((x) & 1) == 0)
#define odd(x) (((x) & 1) == 1)
#define append emplace_back
const i64 MOD = 1000000007;

i32 n, l;
vec<i32> primes;

void init() {
    i32 v = l/(n-1)+1;
    if (v < 2)
        return;

    vec<i32> d;
    if (v >= 7) {
        primes.append(2); primes.append(3);
        i32 acc = 5; i32 b = 1;
        while (acc <= v) {
            d.append(acc);
            if (b > 0)
                acc += 2;
            else
                acc += 4;
            b *= -1;
        }
    } else {
        times(v-1, i)
            d.append(i+2);
    }
    
    i32 w = sqrt(v)+1;
    while (d[0] < w) {
        i32 a = d[0];
        primes.append(a);
        d.erase(d.begin());
        times(d.size(), i) {
            if (d[i] % a == 0)
                d.erase(d.begin() + i);
        }
    }
    times(d.size(), i) {
        primes.append(d[i]);
    }
}

i32 main()
{
    cin >> n >> l;

    init();
    i64 ans = 0;
    times(primes.size(), i) {
        i32 p = primes[i];
        if (p*(n-1) <= l)
            ans += l - p*(n-1) + 1;
    }

    cout << ans << endl;
    return 0;
}
0