Entries

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

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

コメント

コメントの投稿

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

トラックバック

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

Beta Round #26 A問題 - Almost Prime (この記事を編集する[管理者用])

Source

Codeforces Beta Round #26 A問題
Problem description
Beta Round #26の自分の参加記録

問題概要

Almost Primeとは,素数の約数の数が厳密に2個のもの.
1~NまでにAlmost Primeが何個あるか求める問題.
Nは3000以下.

解法

実際に1~NまでのそれぞれがAlmost Primeかどうかを愚直に判定する.

C言語のスパゲッティなコード

このコードは実際に本番でサブミットしたコードです.見にくいと思われます.

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

#define N 4000

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[N];

int is_ok(int n){
  int i,j,k=0;

  rep(i,ps) if(n%p[i]==0){
    while(n%p[i]==0) n/=p[i];
    k++;
  }
  if(k==2 && n==1) return 1;
  return 0;
}

int main(){
  int i,j,k,l,m,n;
  int res=0;

  ps=getPrime(N,p);

  scanf("%d",&n);
  REP(i,5,n+1) res += is_ok(i);
  printf("%d\n",res);

  return 0;
}

コメント

コメントの投稿

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

トラックバック

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

Appendix

Recent Articles

ブログ内検索

Ads


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