Again Prime? No Time.（UVA 10780）

关注

原创最新推荐文章于 2019-02-23 00:11:00 发布 485 阅读 AI 写同款· GEO 优化 ›

Again Prime? No Time. http://acm.hust.edu.cn/vjudge/contest/view.action?cid=31675#problem/D 涨姿势了。。。 // File Name: uva10780.cpp // Author: bo_jwolf // Created Time: 2013年09月16日星期一 21:54:27 #include #include #include # 阅读详情

题目链接
预处理答案ans为+∞，定义函数ok（x，y）表示y！中含有的x的最高次幂的指数，定义函数go（y，x）表示y中含有的x的最高次幂的指数，将m分解质因数，对于每一个质因数i，当前答案now=ok（i，n）/go（m，i），假如now更小则更新ans。
求y！中x的最高次幂的方法为求出并加上1到y中能被x整除的数的个数即y/x，然后将每个数除以x，不能整除的数删去，剩下的数就变成了1到y/x，重复该操作直到y==0。
代码实现：

int ok(int x,int y)
{
    int ans=0;
    while(y)
    {
        ans+=y/x;
        y/=x;
    }
    return ans;
}

附上AC代码：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<queue>
#include<set>
#include<vector>
#include<map>
#include<string>
#include<cmath>
#define pq priority_queue
#define Pi acos(-1.0)
#define MAXX 1000000007
using namespace std;
bool life[6666];
int v[6666],l=0;
int ok(int x,int y)
{
    int ans=0;
    while(y)
    {
        ans+=y/x;
        y/=x;
    }
    return ans;
}
int go(int y,int x)
{
    int ans=0;
    while(y%x==0)
    {
        ans++;
        y/=x;
    }
    return ans;
}
int main()
{
    for(int j=2;j<=5000;j++)
    {
        if(!life[j])
        {
            v[l++]=j;
            for(int k=2;k*j<=5000;k++)
                life[k*j]=1;
        }
    }
    int t,n,m,i=0,now;
    int ans;
    cin>>t;
    while(i<t)
    {
        i++;
        ans=100000000;
        scanf("%d%d",&m,&n);
        cout<<"Case "<<i<<":"<<endl;
        for(int j=0;j<l;j++)
        {
            if(m%v[j]==0)
            {
                ans=min(ans,ok(v[j],n)/go(m,v[j]));
            }
        }
        if(ans)
            cout<<ans<<endl;
        else
            cout<<"Impossible to divide"<<endl;
    }
    return 0;
}