题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=2197
本原串
Problem Description
由0和1组成的串中,不能表示为由几个相同的较小的串连接成的串,称为本原串,有多少个长为n(n<=100000000)的本原串?
答案mod2008.
例如,100100不是本原串,因为他是由两个100组成,而1101是本原串。
答案mod2008.
例如,100100不是本原串,因为他是由两个100组成,而1101是本原串。
Input
输入包括多个数据,每个数据一行,包括一个整数n,代表串的长度。
Output
对于每个测试数据,输出一行,代表有多少个符合要求本原串,答案mod2008.
Sample Input
1 2 3 4
Sample Output
2 2 6 12
题目分析:
长度为n的串有2^n个,减去非本原串就是本原串的数目,非本原串肯定是由长度为k字串不断重复u次得到的,那么k必然是n的约数。 F[n]=2^n-F[k]...
f[n]=2^n-∑f[i](i为n的除1和本身的因子)f[i]代表长度为i的非本源串个数 .
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
const int mod=2008;
int f[100000005];
int pow_mod(int a,int mi)
{
int ans=1;
while(mi)
{
if(mi&1) ans=ans*a%mod;
a=a*a%mod;
mi>>=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;
if(i*i!=n)
ans=(ans-get(n/i)+mod)%mod;
}
}
return f[n]=(ans-2+mod)%mod;
}
int main()
{
int n;
while(~scanf("%d",&n))
printf("%d\n",get(n));
return 0;
}