题目链接:http://lightoj.com/login_main.php?url=volume_showproblem.php?problem=1282
题目大意:给出n^k的前三位和后三位数字
题目思路:后三位简单,直接快速幂%100即可,记得如果不够3位要补0,前三位两种做法
做法一:好像大部分都是采用方法一,方法是设m满足等式,x和y分别表示m的整数部分和小数部分,取一个log,那么
,由于
,所以有了这玩意儿直接*100取个整就行,y就是
减去其取整的结果
做法二:很巧妙的做法,魔改快速幂。由于double精度有限,但是保住最后几位没问题,所以直接每次乘完都一直除,从而保证整数部分一直都小于1000,来保证最后十几位是对的,从而最后三位是对的
以下是代码:
做法一:
#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(ll i=a;i<=b;i++)
#define per(i,a,b) for(ll i=a;i>=b;i--)
#define ll long long
const ll MAXN = 1e6+5;
const ll MOD =1000;
ll powmod(ll x,ll y){
ll rst=1;
for(;y;y>>=1){
if(y&1)rst=rst*x%MOD;
x=x*x%MOD;
}
return rst;
}
int main(){
ll t;
ll n,k;
scanf("%lld",&t);
rep(_,1,t){
scanf("%lld%lld",&n,&k);
double temp=k*log10(n);
double y=temp-(ll)temp;
ll ans1=pow(10,y)*100;
ll ans2=powmod(n,k);
printf("Case %lld: %lld %03lld\n",_,ans1,ans2);
}
return 0;
}
做法二:
#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(ll i=a;i<=b;i++)
#define per(i,a,b) for(ll i=a;i>=b;i--)
#define ll long long
const ll MAXN = 1e6+5;
const ll MOD =1000;
ll powmod(ll x,ll y){
ll rst=1;
for(;y;y>>=1){
if(y&1)rst=rst*x%MOD;
x=x*x%MOD;
}
return rst;
}
ll powmod2(double x,ll y){
double rst=1;
for(;y;y>>=1){
if(y&1)rst=rst*x;
while(rst>=1000)rst/=10;
x=x*x;
while(x>=1000)x/=10;
}
while(rst>=1000)rst/=10;
return rst;
}
int main(){
ll t;
ll n,k;
scanf("%lld",&t);
rep(_,1,t){
scanf("%lld%lld",&n,&k);
ll ans1=powmod2(n,k);
ll ans2=powmod(n,k);
printf("Case %lld: %lld %03lld\n",_,ans1,ans2);
}
return 0;
}