FZU 1759-Super A^B mod C(快速幂+大整数取模+欧拉函数)

题目链接:点击打开链接

题意:计算 a^b %c 但其中b很大,可能会达到10^1000000, 故有降幂公式 a^b %c= a^(b%phi(c)+phi(c)) %c  (b>=phi(c)) 

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define maxn 10100
#define _ll __int64
#define ll long long
#define INF 0x3f3f3f3f
#define Mod 1000000007
#define pp pair
#define ull unsigned long long
using namespace std;
int a,c;char b[1000010];
ll phi(ll n)
{
	ll m=(ll)sqrt(n+0.5),ans=n;
	for(ll i=2;i<=m;i++)
	{
		if(n%i==0)
		{
			ans=ans/i*(i-1);
			while(n%i==0)n/=i;
		}
	}
	if(n>1)ans=ans/n*(n-1);
	return ans;
}
ll pow_mod(ll a,ll n,ll p)
{
	if(n==0)return 1;
	ll ans=pow_mod(a,n/2,p);
	ans=ans*ans%p;
	if(n&1)ans=ans*a%p;
	return ans;
}
void solve()
{
	int len=strlen(b);ll tem=phi(c),sb;
	if(len<=10)
	{
		sscanf(b,"%I64d",&sb);
		if(sb>=tem)
			printf("%I64d\n",pow_mod(a,sb%tem+tem,c));
		else
			printf("%I64d\n",pow_mod(a,sb,c));
		return ;
	}
	ll ans=0;
	for(int i=0;i


你可能感兴趣的:(数学)