結果

問題 No.12 限定された素数
コンテスト
ユーザー k
提出日時 2014-10-30 13:40:34
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,880 bytes
コンパイル時間 810 ms
コンパイル使用メモリ 75,008 KB
実行使用メモリ 23,144 KB
最終ジャッジ日時 2024-12-30 14:00:05
合計ジャッジ時間 4,189 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21 WA * 5
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <vector>
#include <queue>
#include <deque>
#include <stack>
#include <algorithm>
#include <functional>
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <cctype>
#include <string>
#include <cstring>

using namespace std;

typedef long long ll;

bool not_prime[5000001];
int memo[10][348515];

int main() {
  int n;
  cin >> n;
  int a[n];
  bool ba[10];
  memset(ba, false, sizeof(ba));
  for (int i = 0; i < n; i++) cin >> a[i];
  for (int i = 0; i < n; i++) ba[a[i]] = true;


  ll cnt = 0;
  for (int i = 2; i <= 5000000; i++) {
    if (not_prime[i]) continue;
    for (int j = i+i; j <= 5000000; j+=i) {
      not_prime[j] = true;
    }
  }

  vector <int> prime;
  for (int i = 2; i <= 5000000; i++) if (!not_prime[i]) prime.push_back(i);

  memset(memo, 0, sizeof(memo));
  for (int i = 1; i <= prime.size(); i++) {
    for (int j = prime[i-1]; j != 0; j/=10) {
      memo[j%10][i]++;
    }
    for (int j = 0; j < 10; j++) memo[j][i+1] += memo[j][i];
  }

  int ans = -1;
  for (int l = 1, r = 1; r <= prime.size(); r++) {
    bool ok = true;
    for (int i = 0; i < n; i++) {
      if (memo[a[i]][r]-memo[a[i]][l-1] == 0) {
	ok = false;
	break;
      }
    }
    if (ok) {
      bool yes = true;
      for (; l <= r; l++) {
	yes = true;
	for (int j = 0; j < 10; j++) {
	  if ((memo[j][r]-memo[j][l-1] == 0 && ba[j]) ||
	      (memo[j][r]-memo[j][l-1] > 0 && !ba[j])) {
	    yes = false;
	    break;
	  }
	}
	if (yes) break;
      }
      if (yes) {
	int tl = prime[l-1], tr = prime[r-1];
	if (l-1 == 0) {
	  tl = 1;
	}else {
	  tl = prime[l-1-1]+1;
	}
	if (r == prime.size()) {
	  tr = 5000000;
	}else {
	  tr = prime[r-1+1]-1;
	}
	ans = max(ans, tr-tl);
      }else {
	l = r;
      }
    }
  }
  cout << ans << endl;
}
0