求root(N,K)

时间限制：1秒 空间限制：65536K 热度指数：2826

算法知识视频讲解

题目描述

    N<k时，root(N,k) = N，否则，root(N,k) = root(N',k)。N'为N的k进制表示的各位数字之和。输入x,y,k，输出root(x^y,k)的值 (这里^为乘方，不是异或)，2=<k<=16，0<x,y<2000000000，有一半的测试点里 x^y 会溢出int的范围(>=2000000000) 

输入描述:

    每组测试数据包括一行，x(0<x<2000000000), y(0<y<2000000000), k(2<=k<=16)

输出描述:

    输入可能有多组数据，对于每一组数据，root(x^y, k)的值

示例1

输入

复制

4 4 10

输出

复制

4

 

//这道题本来用的一般方法，本来想怎么也不会溢出吧，毕竟long long 都用上了，但是事实证明还是溢出了，所以在网上搜了一下，才发现还是有算法的，快速模取算法和二分模取算法。

主要公式为

N=a0+a1*k+a2*k^2+……an*k^n

      N'=a0+a1+a2+……+an

      N-N'=a1*(k-1)+a2*(k-1)^2+a3*(k-1)^3+......+an*(k-1)^n

      (N-N')%(k-1)=0

      (N'-N'')%(k-1)=0

     .....

      (N(r-1)-N(r))%(k-1)=0

      相加得(N-N(r))%(k-1)=0

      N(r)=N%(k-1)

      故(x^y)%(k-1)就是我们要求的。

      当(x^y)%(k-1)=0时，注意结果为k-1

还看不懂的话可以看这个，解释的比较清楚 https://blog.csdn.net/yangyafeiac/article/details/8707079

//超时算法

/*#include<bits/stdc++.h>

using namespace std;

long long jinzhi(long long temp,long long k)

{

long long count=0;

while(temp!=0)

{

count+=temp%k;

temp=temp/k;

}

temp=count;

return temp;

}

int compare(long long temp,long long k)

{

if(temp<k)

return 0;

return 1;

}

int main()

{

long long x,y,k,count=0,temp;

while(cin>>x>>y>>k)

{

temp=pow(x,y);

while(compare(temp,k))

{

temp=jinzhi(temp,k);

}

cout<<temp<<endl;

}

system("pause");

return 0;

}*/

#include<bits/stdc++.h>

using namespace std;

long long quickpower(long long a,long long b,long long c)//快速模取算法

{

long long count=1;

while(b)

{

if(b%2==1)//奇数要单独乘一次

count=count*a%c;

b>>=1;

a=a*a%c;//偶数直接翻倍

}

return count;

}

int main()

{

long long a,b,c;

while(cin>>a>>b>>c)

{

long long ans=quickpower(a,b,c-1);//为什么要c-1呢，自己看公式

if(ans==0)//等于0说明最终的数和c-1是同一个数

ans=c-1;

cout<<ans<<endl;

}

system("pause");

return 0;

}