hdu 2197 【本原串】

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
  
  

 

 

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<math.h>
#include<map>
using namespace std;
map<int,int>mp;
const int MOD=2008;
int fun(int a,int n){  //a的n次方
   if(n==1)return a;
   int num=fun(a,n/2);
   num=num*num%MOD;
   if(n%2==1)num=num*a%MOD;
   return (num+MOD)%MOD;    
}
int cal(int n){    计算长度为n的本原串的个数
    if(mp[n]!=0)return mp[n];
    mp[n]=fun(2,n)-2;
    for(int i=2;i*i<=n;i++){
        if(n%i==0){
              mp[n]=(mp[n]-cal(i)+MOD)%MOD;
              if(i*i!=n)mp[n]=(mp[n]-cal(n/i)+MOD)%MOD;     
           }           
             
    }
   return mp[n];   //用map来存储
}



int main()
{
    int n;
    while(scanf("%d",&n)!=EOF){
       mp[0]=0;
       mp[1]=2;
       mp[2]=2;
       if(n<=2){cout<<mp[n]<<endl;continue;}
       int ans=cal(n);
       cout<<ans<<endl;                                   
    }    
    return 0;
}