hdu 2197 本原串（递推，快速幂取模）

本原串

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 1135    Accepted Submission(s): 364

Problem Description

由0和1组成的串中，不能表示为由几个相同的较小的串连接成的串，称为本原串，有多少个长为n（n<=100000000)的本原串？
答案mod2008.
例如，100100不是本原串，因为他是由两个100组成，而1101是本原串。

 

Input

输入包括多个数据，每个数据一行，包括一个整数n，代表串的长度。

 

Output

对于每个测试数据，输出一行，代表有多少个符合要求本原串，答案mod2008.

 

Sample Input

  
  

   
   1
2
3
4
  
  

 

Sample Output

  
  

   
   2
2
6
12
  
  

思路：先找出递推方程  假设f[n]表示长度为n的本原串的个数  那么有f[n]=2^n-2-f[i]  其中i是n的约数（不包括1和本身） 减2是因为要减去2^n中的全0和全1

然后就是优化的问题了，首先求2^n肯定是要快速幂的。注意此题不能打表，1E打表实在太慢。

可能是因为此题询问没有很多次，所以只要每次询问时做一次查找就可以，数组记录上上次找过的结果，相同的数字就不用重复查找了。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define LL long long
#define mod 2008
#define N 100000005
int f[N];
int pow_mod(int a,int n)
{
    LL ans=1;
    while(n)
    {
        if(n&1) ans=ans*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return ans;
}
int get(int n)
{
    if(f[n]) return f[n];
    if(n==1) return 2;
    int ans=pow_mod(2,n);
    for(int i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans=(ans-get(i))%mod;
            int t=n/i;
            if(i*i!=n) ans=(ans-get(t)+mod)%mod;
        }
    }
    return f[n]=(ans-2+mod)%mod;
}
int main()
{
    int n;
    while(~scanf("%d",&n))
        printf("%d\n",get(n));
    return 0;
}