【贪心】Calling Extraterrestrial Intelligence Again POJ 1411

题目链接：http://poj.org/problem?id=1411

题目大意：找两个素数p，q满足a/b<=p/q<=1 且p*q<=m，求p*q最大的一组素数对。

第一次想的是筛素数，然后枚举，后来发现没有超时（汗，这让我们写正解的很GG）。

正解是先筛素数，然后用mul数组对于两个素数的积存下较小的那个素数。最后从m向下枚举。

枚举代码：

#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>using namespace std;const int maxm = 50000;int prime[6000], cnt, m, a, b;bool flag[maxm+5];void pre(){for(int i=2;i<=50000;i++)if(!flag[i]){prime[++cnt]=i;for(int j=i*i;j<=maxm&&j>0;j+=i)//j可能会溢出flag[j]=1;}return ;}int main(){pre();while(scanf("%d%d%d",&m,&a,&b),m|a|b){int ans1=0,ans2=0;for(int i=1;i<=cnt;i++)for(int j=i;j<=cnt;j++){int p=prime[i],q=prime[j];if(a*q<=b*p&&p*q<=m){if(p*q>ans1*ans2)ans1=p,ans2=q;}else break;}printf("%d %d\n",ans1,ans2);}}

正解代码：

#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>using namespace std;const int maxm = 100000;int prime[10000], cnt, m, a, b, mul[maxm+5];bool flag[maxm+5];void pre(){for(int i=2;i<=maxm;i++)if(!flag[i]){prime[++cnt]=i;for(int j=i*i;j<=maxm&&j>0;j+=i)//j可能会溢出flag[j]=1;for(int j=1;j<=cnt;j++){if(prime[j]*prime[cnt]>maxm) break;mul[prime[j]*prime[cnt]]=prime[j];}}}int main(){pre();while(scanf("%d%d%d",&m,&a,&b),m|a|b){for(int i=m;i>=1;i--){if(!mul[i]) continue;int p=mul[i],q=i/p;if(a*q<=b*p){printf("%d %d\n",p,q);break;}}}}