結果

問題 No.407 鴨等素数間隔列の数え上げ
コンテスト
ユーザー daku9640
提出日時 2018-04-29 23:41:29
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 37 ms / 1,000 ms
コード長 2,174 bytes
コンパイル時間 1,556 ms
コンパイル使用メモリ 169,872 KB
実行使用メモリ 10,356 KB
最終ジャッジ日時 2024-06-27 23:44:31
合計ジャッジ時間 2,727 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
struct cww{cww(){ios::sync_with_stdio(false);cin.tie(0);}}star;
#define all(c) (c).begin(), (c).end()
#define rall(c) (c).rbegin(), (c).rend()
#define vv(type, c, m, n, i) vector<vector<type>> c(m, vector<type>(n, i))
#define rep(i,a,n) for(int i=(a), i##_len=(n); i<i##_len; ++i)
#define rrep(i,a,n) for(int i=(a), i##_len=(n); i>i##_len; --i)
#define each(x,y) for(auto &(x):(y))
#define len(x) ((int)(x).size())
#define mp make_pair
#define eb emplace_back
#define em emplace
#define pb pop_back
#define fi first
#define se second
#define get(a, i) get<i>(a)
typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef vector<double> vd;
typedef vector<long long> vll;
typedef vector<string> vs;
typedef vector<bool> vb;
typedef pair<int, int> pii;
#define mini min_element
#define maxi max_element
#define sum(a, b) accumulate(a, b)
#define chmax(a,b) if(a<b)a=b
#define chmin(a,b) if(a>b)a=b
#define ce(x, y) (x + y - 1) / y
template <class T> inline void print(T t){cout << t << "\n";}
template <class H, class... T> inline void print(H h, T... t){cout << h << " ";print(t...);}
bool sieve[int(1e7)];
vi AtkinSieve(int limit) {
	vi results;
	for(int i: {2, 3, 5})
		if (i < limit)
			results.eb(i);
	int factor = int(sqrt(limit)) + 1;
	int n;
	rep(i, 1, factor) {
		rep(j, 1, factor) {
			n = 4*(i * i) + (j * j);
			if ((n <= limit) && (n % 12 == 1 || n % 12 == 5))
				sieve[n] ^= true;
			n = 3 * (i * i) + (j * j);
			if ((n <= limit) && (n % 12 == 7))
				sieve[n] ^= true;
			if (i > j) {
				n = 3 * (i * i) - (j * j);
				if ((n <= limit) && (n % 12 == 11))
					sieve[n] ^= true;
			}
		}
	}
	rep(index, 5, factor)
		if (sieve[index])
			for (int jndex = index * index; jndex < limit; jndex += index * index)
				sieve[jndex] = false;
	for (int index = 7; index < limit; index += 2)
		if (sieve[index])
			results.eb(index);
	return results;
}
ll solve() {
	ll n, l; cin >> n >> l;
	if (n * 2 - 1 > l) {
		return 0;
	}
	vi ret = AtkinSieve(int(l / (n - 1) + 1));
	ll ans = 0;
	rep(i, 0, len(ret)) {
		ans += l - ret[i] * (n - 1) + 1;
	}
	return ans;
}
int main() {
	print(solve());
	return 0;
}
0