BZOJ 1008 越狱

Posted on 2016-03-08 23:18  ziliuziliu  阅读(64)  评论(0编辑  收藏

这TM就过了?

公式m^n-m*(m-1)^(n-1)加一个快速幂

#include<iostream>
#include<cstdio>
#include<cstring>
#define mod 100003
using namespace std;
long long n,m;
long long pow(long long x,long long y)
{
long long base=x,ans=1;
while (y)
{
if (y&1) ans=(ans*base)%mod;
base=(base*base)%mod;
y=y>>1;
}
return ans%mod;
}
int main()
{
scanf("%lld%lld",&m,&n);
long long a,b;
a=pow(m,n)%mod;
b=((m%mod)*(pow(m-1,n-1)%mod))%mod;
a=(a-b+mod)%mod;
printf("%lld\n",a);
return 0;
}