记录编号 141110 评测结果 AAAAAAAAAAAAAAAAAAAA
题目名称 [COCI 2011] BROJ 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.166 s
提交时间 2014-11-29 10:54:11 内存使用 1.87 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<bitset>
using namespace std;
typedef long long LL;
const int C=71;
const LL SUP=1000000000;
int checkprime(int p){
	for(int i=2;i<p;i++) if(p%i==0) return false;
	return true;
}
int N,P;
LL mx;
int cnt;
vector<int> prime;
void DFS(LL now,int k,int sgn){
	if(k==prime.size()){
		cnt+=(mx/now)*sgn;
		return;
	}
	if(now*prime[k]<=mx) DFS(now*prime[k],k+1,-sgn);
	DFS(now,k+1,sgn);
}
void work_small(void){
	for(int i=2;i<P;i++){
		if(checkprime(i)) prime.push_back(i);
	}
	LL l=1,r=SUP/P+1;
	while(l<r){
		LL mid=(l+r)>>1;
		cnt=0,mx=mid;
		DFS(1,0,1);
		//cnt=最小素因子为P的个数
		if(cnt<N) l=mid+1;
		else r=mid;
	}
	LL ans=l*P;
	if(ans>SUP) ans=0;
	printf("%lld\n",ans);
}
bitset<13000000> flag;
void work_large(void){
	LL mx=SUP/P;
	cnt=1;
	if(cnt==N){
		printf("%d\n",P);
		return;
	}
	for(LL i=2;i<=mx;i++){
		if(!flag[i]){
			if(i<P){
				for(LL j=i*i;j<=mx;j+=i) flag[j]=true;
			}
			else{
				cnt++;
				if(cnt==N){
					printf("%lld\n",i*P);
					return;
				}
			}
		}
	}
	printf("0\n");
}
int main(){
	freopen("broj.in","r",stdin);
	freopen("broj.out","w",stdout);
	scanf("%d%d",&N,&P);
	if(P<=73) work_small();
	else work_large();
	return 0;
}