SPOJ 9117

题目链接：http://www.spoj.pl/problems/GCPC11A/

——————————————————————————————————————————————

题目描述：

n = k^i ,求一最大i，满足左式。

——————————————————————————————————————————————

题目思路：

1、筛素数。把k的因子筛出来。

2、通过计算( [n/p ] + [n/p^2] + [n/p^3] ....)/(k中p因子的个数), 并对所有的k的质因子计算一遍，求一个最小的值，即为本题解。

——————————————————————————————————————————————

题目细节：

1、对k 筛质数的时候要注意，筛到最后，所有的质数都筛除来以后，退出循环的时候，temp仍有可能不为1，可能为一大素数，所以要把剩下的继续存到数组里。

2、long long int的最大值为 k*10^18，所以 n可能很接近这个值。所以如果将

  for (int i=0;i<num;i++)
            kans[i]=get(n,ks[i]);

改成：

for(i = 0;i<num;i++)
{
    for(j = 1;ks[i]^j <=n;j++)
      kans[i] += ks[i]^j;
}

在n很大的时候就会成过了，即由于ks[i]^j 很大，再多乘以一个ks[i]，会超过longlongint 的范围，这样又会满足小于n的循环，答案会多算很多遍。

今后当n很大的时候一定要注意这个问题。

3、

long long get(long long n,long long int  p)

对于本函数，参数传递的时候要注意，由于k最大的为10^12，是大于int型的，k本身可能是素数，且k的因子里可能有大于10^6的因子（这点在第一条里就应该注意到了），所以p很可能是超过int的范围的，需要定义成long long int

4、开数组时，对于knum，ks，kans，不必将数组开的很大。将前几个素数连乘，乘到37的时候就已经可以超过10^12次方了，此时只用了前11个素数。所以数组开到20就满够了。

——————————————————————————————————————————————

源代码：

#include <iostream>
#include<stdio.h>
#include<stdlib.h>
#include<cstring>

using namespace std;

#define maxn 1000010
#define minfun(a,b) ((a)>(b)?(b):(a))

int prime[maxn];
int num_prime = 0;
int a[maxn];
long long int knum[100],ks[100],num = 0;
long long int kans[100];

long long get(long long n,long long int  p)
{
    if (!n) return 0;
    return n/p+get(n/p,p);
}

int main()
{
    long long int n = 0,k = 0;
    int t = 0;
    int i = 0,j = 0;
    long long int min = 0;

    memset(a,0,sizeof(a));
    for(i = 2;i<=maxn;i++)
    {
        if(!a[i])
          prime[num_prime++] = i;

        for(j = 0;j<num_prime &&i*prime[j]<=maxn;j++)
        {
             a[i*prime[j]] = 1;
             if(!(i%prime[j]))
               break;
        }
    }

    scanf("%d",&t);
    while(t--)
    {
        num = 0;
        scanf("%lld%lld",&n,&k);
        long long int temp = k;
        for(i = 0;i<num_prime;i++)
        {
            if(temp%prime[i] == 0)
            {
                ks[num] = prime[i];
                knum[num] = 1;
                temp = temp/prime[i];
                num++;
            }

            while(temp%prime[i] == 0)
            {
               knum[num-1]++;
               temp = temp/prime[i];
            }

            if(temp == 1)
              break;
        }

        if(temp!=1)                                          //本处加星星！！
        {
            ks[num] = temp;
            knum[num] = 1;
            num++;
        }
        for (int i=0;i<num;i++)
            kans[i]=get(n,ks[i]);

        min = kans[0]/knum[0];
        for(i = 1;i<num;i++)
           min = minfun(min,kans[i]/knum[i]);

        printf("%lld\n",min);
    }

    return 0;
}