1. 2019南京网络赛super_log: 求(fuck[0]^fuck[1]^fuck[2]^...^fuck[b-1] mod m )>=b 其中fuck数组=a; 数据范围:1<=a<=1e6: 0<=b<=1e6,1<=m<=1e6
上代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAX = 1e6+100;
ll phi[MAX],fuck[MAX];
void getphi(){
phi[1]=1;
for (ll i = 2; i < MAX;i++){
if(!phi[i]){
for (ll j = i; j < MAX;j+=i){
if(!phi[j]) phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}
}
}
}
ll quick(ll a,ll b,ll c){
ll ans=1;
if(a>=c) a=a%c+c;
while(b){
if(b&1){
ans*=a;
if(ans>=c) ans=ans%c+c;
}
a*=a;
if(a>=c) a=a%c+c;
b>>=1;
}
return ans;
}
ll a,b,m;
ll zsxhj(ll id,ll mod){// 求(fuck[0]^fuck[1]^fuck[2]^...^fuck[b-1] mod m )>= b
if(b==0) return 1%m;
if(mod==1) return 1;
if (id==b-1){
if(fuck[id]>=mod) return fuck[id]%mod+mod;
else return fuck[id];
}
else return quick(fuck[id],zsxhj(id+1,phi[mod]),mod);
}
int main(){
getphi();
int t;
scanf("%d",&t);
while(t--){
scanf("%lld%lld%lld",&a,&b,&m);
for (int i = 0; i < b;i++){
fuck[i]=a;
}
ll w=zsxhj(0,m);
printf("%lld\n",w%m);
}
return 0;
}
2. FZU-1759 求a^b%mod 数据范围:(1<=a,c<=1000000000,1<=b<=10^1000000).
上代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll get_euler(ll n){
ll res=n,a=n;
for (ll i = 2; i*i<=a;i++){
if(a%i==0){
res=res/i*(i-1);
while(a%i==0) a=a/i;
}
}
if(a>1) res=res/a*(a-1);
return res;
}
ll quick(ll a,ll b,ll c){
ll ans=1;
if(a>=c) a=a%c+c;
while(b){
if(b&1){
ans*=a;
if(ans>=c) ans=ans%c+c;
}
a*=a;
if(a>=c) a=a%c+c;
b>>=1;
}
return ans;
}
int main(){
ll a,c;
string b;
while(~scanf("%lld",&a)){
cin >> b;
scanf("%lld",&c);
ll p=get_euler(c);
int len=b.size();
ll sum=0;
for (int i = 0; i < len;i++){
sum=sum*10+(b[i]-'0');
sum%=p;
}
printf("%lld\n",quick(a,sum,c)%c);
}
return 0;
}
3. 计蒜客 Exponial 求(n^(n-1)^(n-2)^...^1 )mod m 数据范围:1<=n<=1e9,1<=m<=1e9
上代码:
#include<bits/stdc++.h>
typedef long long ll;
int get_euler(int x){
int ans=x;
for(int i=2;i*i<=x;i++){
if(x%i==0){
ans-=ans/i;
while(x%i==0) x/=i;
}
}
if(x>1) ans-=ans/x;
return ans;
}
ll quick(ll a,ll b,ll c){
ll ans=1;
if(a>=c) a=a%c+c;
while(b){
if(b&1){
ans*=a;
if(ans>=c) ans=ans%c+c;
}
a*=a;
if(a>=c) a=a%c+c;
b>>=1;
}
return ans;
}
ll exponial(int x,int y){//求(x^(x-1)^(x-2)^...^1 )mod y
if(x==1) return 1%y;
if(x%y==0) return quick(x,x-1,y);
return quick(x,exponial(x-1,get_euler(y)),y);
}
int main(){
int n,m;
while(~scanf("%d%d",&n,&m)){
printf("%lld\n",exponial(n,m)%m);
}
return 0;
}