結果
問題 |
No.407 鴨等素数間隔列の数え上げ
|
ユーザー |
|
提出日時 | 2016-08-08 20:09:29 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 209 ms / 1,000 ms |
コード長 | 1,132 bytes |
コンパイル時間 | 1,438 ms |
コンパイル使用メモリ | 166,548 KB |
実行使用メモリ | 45,312 KB |
最終ジャッジ日時 | 2024-12-15 23:03:13 |
合計ジャッジ時間 | 7,584 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 31 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS #define _USE_MATH_DEFINES #include "bits/stdc++.h" #define REP(i,a,b) for(i=a;i<b;++i) #define rep(i,n) REP(i,0,n) #define ll long long #define ull unsigned ll typedef long double ld; #define ALL(a) (a).begin(),(a).end() #define ifnot(a) if(not a) #define dump(x) cerr << #x << " = " << (x) << endl using namespace std; // #define int ll bool test = 0; int dx[] = { 0,1,0,-1 }; int dy[] = { 1,0,-1,0 }; #define INF (1 << 28) ull mod = (int)1e9 + 7; //..................... #define MAX (int)1e7 + 5 static int prime[MAX + 5]; static int isprime[MAX + 5]; void Eratosthenes(int N) { for (int i = 0; i < N; i++) { isprime[i] = 1; } for (int i = 2; i < sqrt((double)N); i++) { if (isprime[i]) { for (int j = 0; i * (j + 2) < N; j++) { isprime[i *(j + 2)] = 0; } } } int a = 0; for (int i = 2; i < N; i++) { if (isprime[i]) { prime[a++] = i; } } } signed main() { ll i, j, k, l; ll N,L; cin >> N >> L; Eratosthenes(MAX - 3); ll cnt = 0; rep(i, L) { if (L - prime[i] * (N - 1) < 0) break; cnt += L - prime[i] * (N - 1) + 1; } cout << cnt << endl; return 0; }