結果

問題 No.12 限定された素数
コンテスト
ユーザー pekempey
提出日時 2015-08-18 17:12:17
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 61 ms / 5,000 ms
コード長 1,544 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,261 ms
コンパイル使用メモリ 163,416 KB
実行使用メモリ 11,196 KB
最終ジャッジ日時 2024-11-24 08:28:02
合計ジャッジ時間 3,379 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define rep(i, a) for (int i = 0; i < (a); i++)
#define rep2(i, a, b) for (int i = (a); i < (b); i++)
#define repr(i, a) for (int i = (a) - 1; i >= 0; i--)
#define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--)
using namespace std;
typedef long long ll;
const ll inf = 1e9;
const ll mod = 1e9 + 7;

const int maxn = 5000001;
bool isprime[maxn];
vector<int> prime;
vector<int> digit;

int decomp(int n) {
	int res = 0;
	while (n) {
		res |= 1 << n % 10;
		n /= 10;
	}
	return res;
}

void sieve() {
	fill(isprime, isprime + maxn, true);
	for (ll i = 2; i < maxn; i++) {
		if (isprime[i]) {
			prime.push_back(i);
			digit.push_back(decomp(i));
			for (ll j = i * i; j < maxn; j += i) {
				isprime[j] = false;
			}
		}
	}
}

bool issubset(int a, int b) {
	return (a | b) == b;	
}

int main() {
	sieve();
	int N;
	cin >> N;
	int A = 0;
	rep (i, N) {
		int temp;
		cin >> temp;
		A |= 1 << temp;
	}

	int r = 0;
	int curr = 0;
	vector<int> num(10);
	int ans = -1;
	rep (i, prime.size()) {
		r = max(r, i);
		while (r < prime.size() && issubset(curr | digit[r], A)) {
			curr |= digit[r];
			rep (j, 10) {
				if (digit[r] >> j & 1) num[j]++;
			}
			r++;
		}

		if (curr == A) {
			int ub, lb;
			if (r == prime.size()) {
				ub = 5000000;
			} else {
				ub = prime[r] - 1;
			}
			if (i == 0) {
				lb = 1;
			} else {
				lb = prime[i - 1] + 1;
			}
			ans = max(ans, ub - lb);
		}

		curr = 0;
		rep (j, 10) {
			if (digit[i] >> j & 1) num[j]--;
			if (num[j] > 0) curr |= 1 << j;	
		}
	}

	cout << ans << endl;
}
0