結果
問題 |
No.407 鴨等素数間隔列の数え上げ
|
ユーザー |
![]() |
提出日時 | 2019-06-13 18:12:23 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 62 ms / 1,000 ms |
コード長 | 1,389 bytes |
コンパイル時間 | 1,369 ms |
コンパイル使用メモリ | 166,092 KB |
実行使用メモリ | 13,312 KB |
最終ジャッジ日時 | 2024-10-15 00:00:51 |
合計ジャッジ時間 | 4,552 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 31 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } #define FOR(i,a,b) for(ll i=(a);i<(b);++i) #define ALL(v) (v).begin(), (v).end() #define p(s) cout<<(s)<<endl #define p2(s, t) cout << (s) << " " << (t) << endl #define br() p("") #define pn(s) cout << (#s) << " " << (s) << endl const int N_MAX = 1e7 + 10; bool is_prime[N_MAX]; // エラトステネスのふるい void Eratosthenes(){ FOR(i, 0, N_MAX){ is_prime[i] = true; } is_prime[0] = false; is_prime[1] = false; for(ll i=2; i*i<=N_MAX; i++){ if(is_prime[i]){ int p = i; int step_p = p; // それの倍数を消していく while(p < N_MAX){ p += step_p; if(p<N_MAX){ is_prime[p] = 0; } } } } } int main(){ cin.tie(0); ios::sync_with_stdio(false); // input ll N, L; cin >> N >> L; Eratosthenes(); ll count = 0; ll max_interval = L / (N-1); FOR(p, 2, max_interval+1){ if(!is_prime[p]) continue; ll len = p*(N-1); ll rest = L - len; count += rest+1; } p(count); return 0; }