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

 

思路：题目要求找出本原串个数，那我们也可以找出非本原串的个数，再用2 ^ n减去就好了， 容易得出，组成非本原串的字符串的长度一定是n的因子，所以我们枚举所有因子，然后再处理一下重复情况就好了

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
LL INF = 1000000000;
const LL MAX = 1e8 + 50;

LL mod = 2008;
LL qpow(LL b, LL n){
    LL re = 1;
    b %= mod;

    while(n){
        if(n & 1){
            re = re * b % mod;
        }
        n >>= 1;
        b = b * b % mod;
    }

    return re % mod;
}
LL a[MAX];
LL solve(LL n){
    //prLLf("n = %d\n", n);
    if(a[n] != 0){
        return a[n];
    }
    a[n] = qpow(2, n) - 2;

    for(LL i = 2; i * i <= n; i++){
        if(n % i == 0){
            a[n] = (a[n] - solve(i) + mod) % mod;
            if(i * i != n){
                a[n] = (a[n] - solve(n / i) + mod) % mod;
            }
        }
        
    }
    return a[n];
}
int main()
{
    LL n;
    a[0] = 0;
    a[1] = 2;
    a[2] = 2;
    while(~scanf("%I64d", &n)){
        if(n > 2){
            a[n] = solve(n);
        }
        printf("%I64d\n", a[n]);
    }
    return 0;
}