nyoj 520 最大素因子(打表)
原创
©著作权归作者所有:来自51CTO博客作者leibniz_zhang的原创作品,请联系作者获取转载授权,否则将追究法律责任
最大素因子
时间限制:1000 ms | 内存限制:65535 KB
难度:2
描述
GreyAnts最近正在学习数论中的素数,但是现在他遇到了一个难题:给定一个整数n,要求我们求出n的最大素因子的序数,例如:2的序数是1,3的序数是2,5的序数是3,以此类推. 研究数论是需要很大的耐心的,为了惩罚那些没有耐心读完题目的童鞋,我们规定:1的最大素因子序数是0.
输入
有多组测试数据,每一行输入一个数字n.(0<n<=1000000)
输出
在接下来的一行,输出结果.
样例输入
2
3
4
5
样例输出
1
2
1
3
来源
hdu改编
上传者
ACM_李如兵
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MA 1000100
int fafe[1000001];
int shu[MA];
int zhi[MA];
int main()
{
memset(fafe,0,sizeof(fafe));
memset(zhi,0,sizeof(zhi));
fafe[1]=1;long long i,j,kp;
for (i=2;i<=1000000;i++)
{
if (!fafe[i])
{
/* for (j=2;j*i<MA;j++)
fafe[i*j]=1;*/
/* for(j=i;i*j<MA;j++)
fafe[i*j]=1;*/
for (j=i*i;j<MA;j=j+i)
fafe[j]=1;
}
}
kp=1;
for (i=2;i<MA;i++)
{
if (!fafe[i])
shu[i]=kp++;
}
kp=0;shu[0]=0;
for (i=1;i<MA;i++)
{
if (!fafe[i])
{
for (j=1;j*i<MA;j++)
zhi[i*j]=shu[i];
}
}
while(~scanf("%d",&kp))
printf("%d\n",zhi[kp]);
return 0;
}