LightOJ - 1282 Leading and Trailing(快速幂+数学)

题目链接：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;
}