UVa1521 GCD Guessing Game（唯一分解定理+素筛）

题目链接

---

题目大意

给出一个数$n$，现在让我们猜$1-n$内的一个数$y$，唯一的提示是我们每猜一个数$x$，就会提示$gcd(x,y)$的值，现在让我们求给定$n$的最小猜的次数。

解题思路

首先根据质因数分解定理，除了$1$以外任何一个正整数都可以表示成若干个质因数的乘积，那么显然我们一旦猜了一个质数，那么含该质因子的合数都会因为给出的$gcd$被判断出来是不是在下一个猜的范围内。但是这还不够，不难发现最坏的情况是猜完了所有质数才确定了答案（因为$n$范围内的质数起码构成了包括$n$在内的合数，也就是$n$内的合数分解质因子的范围就在$n$以内，具体说是$sqrt(n+0.5)$）。那怎么才能更快呢，一次猜多个质数！比如$6=2\ast 3$，那么我们猜$6$的时候，因为$6$含有$2$和$3$的质因子，因此答案是否含有$2$和$3$就可以确定，如果不含有那么就直接排除了所有含$2$和$3$质因子的数。那么如果想更快，就必须排除更多的数，那么问题就变成了$n$以内的多个质数相乘不超过$n$就可以被分为一组，也就是判断一次即可，能分成多少组就能判断多少次（当所有质数都判断了$n$以内的所有数也都判断了）

代码：

#include <iostream>
#include <cstring>
#include <math.h>
using namespace std;
const int maxn=1e4+10;
int prime[maxn];
bool isprime[maxn];
int num,n;

void init(){  //埃氏筛法
    memset(isprime,true,sizeof(isprime));
    num=0;
    for(int i=2;i<maxn;i++){
        if(isprime[i]){
            prime[num++]=i;
            for(int j=i*i;j<maxn;j+=i)
                isprime[j]=false;

        }
    }
}

int main()
{
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    init();
    while(cin>>n){
        int l,r,k=1,ans=0;
        while(prime[k]<=n){  //找最大的小于n的质数的下标
            if(k>=num) break;
            k++;
        }
        l=0,r=k-1;
        while(l<=r){
            int res=prime[r];
            while(prime[l]*res<=n) res*=prime[l++];
            ans++;
            r--;
        }
        cout<<ans<<endl;
    }
    return 0;
}