Entries

スポンサーサイト (この記事を編集する[管理者用])

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/2336-e9d9abc4
この記事にトラックバックする(FC2ブログユーザー)

Round #147 DIV2 C問題 - Primes on Interval (この記事を編集する[管理者用])

Source

Codeforces Round #147 DIV2 C問題 (1500pt)
Problem description

問題概要

正整数a, b, kが与えられる.bは10^6以下.
区間[a, b]の任意のL個の連続する整数は,少なくてもk個の素数を含む.
そのような最小のL (b-a+1以下) を求める問題.
存在しないならそれを指摘する.

解法

最初の数字を固定したら,k個先の素数の位置は単調増加なのでしゃくとりメソッドできて,それをやった.端に注意.
Lについて2分探索して,累積和を用いてある区間の素数の数をO(1)で求めれるようにしておくのがミスしなさそう.
答えは,範囲内の隣合う素数の差の最大値,区間の最初と最初の素数までの差,区間の最後と最後の素数までの差を使ってかけるので,それでも良い.

C言語のスパゲッティなコード
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

int getPrime(int n,int p[]){int i,j,n2=n/2;rep(i,n2)p[i]=1;for(i=3;i*i<n;i+=2)if(p[i>>1])for(j=(i*i)>>1;j<n2;j+=i)p[j]=0;j=1;p[0]=2;REP(i,1,n2)if(p[i])p[j++]=i*2+1;return j;}
int ps, p[1000010], isp[1000010];

int main(){
  int i,j,k,l,m,n;
  int a, b, res, sum;

  ps = getPrime(1000000, p);
  rep(i,1000001) isp[i] = 0;
  rep(i,ps) isp[p[i]] = 1;
  
  scanf("%d%d%d",&a,&b,&k);

  m = 0;
  REP(i,a,b+1) m += isp[i];
  if(m<k){ puts("-1"); return 0; }

  res = 0;
  j = a; sum = 0;
  REP(i,a,b+1){
    while(j <= b && sum < k){
      sum += isp[j];
      j++;
    }
    if(sum < k){
      if(res < j-i+1) res = j-i+1;
      break;
    }
    if(res < j-i) res = j-i;
    sum -= isp[i];
  }

  printf("%d\n",res);

  return 0;
}

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/2336-e9d9abc4
この記事にトラックバックする(FC2ブログユーザー)

Appendix

Recent Articles

ブログ内検索

Ads


(プライバシーポリシー)
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。