HDU Problem1452 Happy 2004（费马小定理+逆元）

题目：

这题，内容与题目严重不符，说好的happy呢？这道题对初学者一点都不友好，尤其是像我这种没什么知识储备的人，但对于学过数论那就随便过了。

分析：

题目要先求2004^X，再求所得数的约数和，再对29取余。input X非常大，肯定会爆，还要求约数和？怎么办？
从后往前走：约数和有公式
S=(p1⁰+p1¹+…+p1^k1)(p2⁰+…+p2^k2)…(pn⁰+…+pn^kn)
(其中p1、p2、……、pn为所求数的质因数，k1、k2…kn是其次数)
根据唯一质数分解定理可以把 p 求出来，剩下的就好求了，求质因数可以用试除法：

    for(long long i=2;i<=n/i;i++){
      if(n%i == 0){
        p[cnt++]=i; //记录质因数
        while(i%n == 0){
          nub[i]++;
          n/=i;
        }
      }
   }
   if(n > 1) {  p[cnt++]=n; nub[n]++;//别忘记最后没除完的数哦

好了，算到这边题目应该变简单了吧，再用快速幂和费马小定理吧 n=2004^X 算出来，题目就算完成了。可是有没有发现，题目似乎没给我们数字对2004^X去取模。
但这种思路是错了，我在这条路上走了一天。。。。
重新来看，2004^X的质因数不就是2004的质因数做X次幂吗，2004=2²*3¹*167¹；所以2004^X的质因数就只有2、3、167，次数分别为2x，x，x；再用等比数列通项公式就可以将其求出来，这里要用到费马小定理和逆元，因为通项公式是（aⁿ⁺¹-1)/(a-1),为了不爆long long 必须要取除数的逆元，如果a=3，要除以2，2 mod（29）逆元是15，同理 167-1=166的逆元是18，不多说，上代码

#include<stdio.h> 

long long slow(long long n,int x){//快速幂
	long long ans=1;
	n%=28;     //费马小定理
    for(int i=1;i<=n;i++)
    {
	 ans=(ans*x)%29;
    }
	return ans-1;
}

int main(){
	long long n;
	while(scanf("%lld",&n) && n){
		long long p1,p2,p3,ans;
	      p1=slow(2*n+1, 2);
	      p2=slow(n+1, 3)*15;   //15是2的逆元
	      p3=slow(n+1, 167)*18; //18是166的逆元
	    ans=(p1*p2*p3)%29;
		printf("%lld\n",ans);
	}
	return 0;
}