HDU 2197 本原串 （规律+快速幂）

本原串

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

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
  
  

 

Author

scnu

 

此题用反面情况求解，长度为n的串有2^n中，减去非本原串，非本原串肯定是由长度为v字串不断重复u次得到的，那么v必然是n的约数。

应该可以打表加快速度，但我懒得写，可以ac。

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
#define ll long long
const ll p = 2008;
ll qkm(ll base, ll mi)
{
    ll ans=1;
    while(mi)
    {
        if(mi&1) (ans=ans*base)%=p;
        mi>>=1;
        (base*=base)%=p;
    }
    return ans;
}
ll d(ll x)
{
    if(x==1) return 2;
    ll ans=qkm(2,x)-2;
    for(ll i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            (ans=ans-d(i)+p)%=p;
            if(i*i!=x) (ans=ans-d(x/i)+p)%=p;
        }
    }
    return ans;
}

int main()
{
    ll n;
    while(~scanf("%lld",&n))
    {
        printf("%lld\n",d(n));
    }
    
}