【HDU - 2197】本原串 （数学）

kuangbin大神的博客

由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

思路：

一道拖了好久的题目，今天总算是看明白了。

一个串如果是非本原串，那么他一定有循环节，这个循环节的长度就是它因子的长度。

所以求本原串的个数可以用长度为n的串的总个数减去非本原串的个数。

f[n]=2^n-sum-2。（f[i]表示长度为i的串有多少个本原串）

sum是n的所有因子的本原串的个数的加和。

因为一个非本原串可以找到一个本原串组成他。本原串找不到个，所以只需要减去每个因子大小的本原串的个数就可以。因为非本原串也会被包含在内了，如果在减去非本原串，那就多算了。

比如10101010

可以由1010组成，也可以有10组成，但是算10组成时，已经包含了1010组成的情况了。

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=2008;
map<int,int> mp;
int qsm(ll a,ll b)
{
	ll t=1;
	while(b)
	{
		if(b&1)
		t=(t*a)%mod;
		a=(a*a)%mod;
		b>>=1;
	}
	return (int)(t%mod);
}
int solve(int x)
{
	if(mp[x]!=0) return mp[x];
	int ans=qsm(2,x);
	ans=(ans-2)%mod;
	for(int i=2;i*i<=x;i++)
	{
		if(x%i!=0) continue;
		if(i*i==x)
		{
			ans-=solve(i);
			ans%=mod;
		}
		else
		{
			ans-=solve(i);
			ans-=solve(x/i);
			ans%=mod;
		}
	}
	return mp[x]=(ans+mod)%mod;
}
int main()
{
	mp[1]=2;
	int n;
	while(~scanf("%d",&n))
	{
		printf("%d\n",solve(n));
		
	} 
}