FZU 1759-Super A^B mod C(高速べき乗+大整数型取り+Euler関数)

1394 ワード

タイトルリンク:クリックしてリンクを開く
标题: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
",pow_mod(a,sb%tem+tem,c));
		else
			printf("%I64d
",pow_mod(a,sb,c));
		return ;
	}
	ll ans=0;
	for(int i=0;i

matlabテキストファイルを読み込む----textscan

phpcmsの基本知識