uva 11029 Leading and Trailing

数论题:幂取模

给出a和n,求a^n的前3位和后3位

关于求后3位,就是一个二分并且模1000即可,至于怎么求前3位,需要一些数学分析

我们来看一大数n

n可以写成n=10^b,这个a是个小数
所以b=i+d,i是整数部分,d是数部分
则n=10^(i+d)=10^i * 10^d
这里很重要,10^i其实有什么用?i是整数,所以10^i一定是1000…………000,然后*10^d,d是一个小于1的小数(别忘了是b的小数部分)
所以决定n这个大数会出现什么数字的,是10^d,10^i只是给10^d起到往后移动小数点的作用
所以我们要知道10^d
另外我们知道0<=d<1,所以 10^0 <= 10^d < 10^1   即1<=10^d<10

这样只需要10^d*100然后取整数部分,就能得到一个3位整数,这3个数字正是我们要的答案

a=10^k , a^n=(10^k)^n=10^k*n=10^i * 10^d    ,  这里i是k*n的整数部分,d是k*n的小数部分

要得到k就使用库函数fmod即可    d=fmod(k*n , 1) , 返回k*n/1的余数,也就是小数部分

另外k=log10(a) , 写在一起就是   d=fmod( log10(a)*n , 1) 

然后就是求  100*10^d=10^(d+2) 

//幂取模,保存后3位和前3位
#include <cstdio>
#include <cmath>
#define M 1000

int pow_mod(int a , int n)
{
    if(n==1)   return a%M;
    int x=pow_mod(a,n/2);
    long long ans=(long long)x*x%M;
    if(n&1)   ans=ans*a%M;
    return (int)ans;
}

int main()
{
    int T,a,n;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&a,&n);
        int left,right;
        left=(int)pow( 10.0 , 2 + fmod(n*log10(1.*a) , 1) );
        right=pow_mod(a,n);
        printf("%d...%03d\n",left,right);
    }
    return 0;
}
原文地址:https://www.cnblogs.com/scau20110726/p/2891054.html