数字游戏
Time Limit:1000MS Memory Limit:1024K
Description:
最近Catcher对数字游戏很感兴趣,他发明了一种新的游戏,给出一个数N,问N是否能表示成某个正整数X的K次幂(K>1),N可能有多种表示方法,请找出最大的X并输出相应的K。例如 16=2^4=4^2,64=4^3=2^6=8^2则16应表示为4^2,64应表示为8^2。
Input:
每行一个正整数N,输入文件以0为结束标志。(0<N<10^8)
Output:
每行有两个整数,如果能表示,则输出X K,(中间用一个空格隔开)如果不能,则输出0 0;
Sample Input:
5
4
16
27
0
Sample Output:
0 0
2 2
4 2
3 3
解法一:最初想法,也是最容易想到,算法最简单的,就是从i=N/2出发,找到一个可以让N整除的,然后继续循环整除下去,只要最后出现1,那么计算它的次数,答案就出来了。
#include<iostream>
using namespace std;
int main()
{
int n;
cin >> n;
while(n)
{
int temp = 0 ;
int time = 0;
for(int i = n/2;i >1;i --)//从n/2开始,寻找能被n整除的数
{
if(n % i == 0)//当i能被n整除,便开始进行循环整除
{
temp = n / i;
time ++;
while(temp %i == 0)
{
temp = temp /i;
time ++ ;
}
if(temp == 1)//最后结果是为1,即n是i的k次方
{
cout << i <<" " << time << endl;
break;
}
else //结果不是1,即n不是i的k次方,则全部归0,重新计算
{
temp = 0;
time = 0;
}
}
}
if(temp == 0)//如果找不到相应的i与k,则输出0 0
cout << temp << " "<<temp << endl;
cin >> n;
}
return 0;
}
解法二:把N=i^k,则logi(N)=k,而logi(N)=logN/logi,可以用一个for循环,从根号i=N(sqrt(N*10))开始找到一个整数的temp,即是k=temp,答案即出来。
原本想用temp-(int)temp==0来判断temp是否为整数,结果发现WA。
在计算方法中学过,temp-(int)temp==0不一定成立,因为计算log(N)/log(i)是会有误差的,应该改为fabs(temp-(int)temp)<e,而这个e要设置为一个充分小的数,这个可以根据题目设的数字范围找到。比如,设log(N+1)与log(N)之差的最小区别,如e=log(10000)/log(100)-log(9999)/log(100);
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
int n;
while(cin>>n&&n)
{
double temp;
bool flag=1;
for(int i=sqrt(n*1.0);i>=2;i--)//从i=sqrt(n*1.0)开始寻找
{
temp=log(n)/log(i);
if(fabs(temp-int(temp+0.5))<(log(10000)/log(100)-log(9999)/log(100)))//判断temp是否为整数,如果是,即是答案,则输出
{
flag=0; //判断是否找到i
cout << i << " "<<temp <<endl;
break;
}
}
if(flag) //找不到i,即没有答案,输出0 0
cout<<"0 0"<<endl;
}
return 0;
}
解法三:先构建结构体m,包括x,y,则题目要求的x与k,然后用一个数组把1<n<10000中的所有数的x与y装起来,如果没有相应的x与y,则为0。
#include<iostream>
using namespace std;
struct m
{
int x;
int y;
};
int main()
{
m a[10000];
for(int k=1;k<10000;k++)//将数组a全附上值为0
{
a[k].x=0;
a[k].y=0;
}
for(int i=2;i<100;i++)//从2开始,计算i的j次方,并把相应的i,j放进数组内
{
int temp = i*i;
for(int j=2;temp<10000;j++,temp*=i)//temp为i的K次方,并根据题目要求,小于10000
{
a[temp].x=i; //temp可能不断被赋予新的值,因为题目要求输出最大的X
a[temp].y=j;
}
}
int n;
while(cin>>n&&n)
{
cout << a[n].x<<" "<<a[n].y<<endl;
}
return 0;
}