【清华夏令营2016模拟5.31】图森破

题目

Description
小火车励志成为一名辣鸡出题人，但是要成为一名辣鸡出题人，代码必须跑得比谁都快，这样就能把他们都卡常数了！为了锻炼自己，他找到了一位长者——罗长者，罗长者说：“你啊，toosimple！不要想弄一个大新闻，说现在已经‘钦定’了，然后把我批评一番。”小火车坐在高高的骨灰旁边，听长者讲那钦定的事情。
如果一个n位数X（可以有前导0）对所有的k(1<=k<n)都满足X10k mod 10n>X，长者就会钦定这个X，然后长者会产生n2的愉悦度。长者会指定一个n，然后钦定满足条件的1位数，2位数，…，n位数。求长者的愉悦度之和，答案模258280327（2317+1，一个质数）。

Input
一行一个正整数表示n。

Output
一行一个整数表示答案。

Sample Input
2

Sample Output
190

Data Constraint
所有的一位数都满足，共10个。
满足条件的两位数的有01到09,12到19,23到29,34到39,45到49,56到59,67到69,78,79,89总共45个。

Hint
对于10%的数据，n=7。
对于30%的数据，n<=10000。
对于60%的数据，n<=10000000。
对于100%的数据，n<=10000000000。

思路

考虑一个串，如果它没有循环同构串和它相同，那么它的最小循环同构串就是一个合法的答案。

记f[n]为长度为n的没有循环同构串和自己相同的串个数，那么f[n]/n就是长度n时的答案。

f[n]=10ⁿ−∑d|n and d≠n f[d]

直接统计10…的贡献，不难得到：
f[n]=∑d|n10^d∗μ(n/d)

杜教筛处理μ(i)∗i的前缀和，后面是高中数列题。

代码

#include<bits/stdc++.h>
using namespace std;

const int mod = 258280327;

ll power(ll x,ll y) {
	ll s = 1;
	for(; y; y /= 2,x = x * x % mod)
		if(y & 1) s = s * x % mod;
	return s;
}

const ll ni2 = power(2,mod - 2);

const int N = 1e6 + 5;

int bz[N],p[N],p0,mu[N];
ll s[N];

void init(int n) {
	mu[1] = 1;
	fo(i,2,n) {
		if(!bz[i]) p[++ p0] = i,mu[i] = -1;
		for(int j = 1; i * p[j] <= n; j ++) {
			int k = i * p[j] % mod;
			bz[k] = 1;
			if(i % p[j] == 0) {
				mu[k] = 0; break;
			}
			mu[k] = -mu[i];
		}
	}
	fo(i,1,n) s[i] = (s[i - 1] + mu[i] * i) % mod;
}

ll n;

ll cn1(ll x) {
	x %= mod;
	return x * (x + 1) % mod * ni2 % mod;
}

const ll ni9 = power(9,mod - 2);

ll calc(ll x) {
	ll s = (power(10,x + 1) - 10) * ni9 - power(10,x + 1) * (x % mod);
	s = (-s % mod + mod) * ni9 % mod;
	return s;
}

const int M = 1960817;

struct hash {
	ll h[M],f[M];
	int fi[M],nt[M],tot;
	ll &operator [] (ll n) {
		int y = n % M;
		for(int p = fi[y]; p; p = nt[p])
			if(h[p] == n) return f[p];
		nt[++ tot] = fi[y],fi[y] = tot,h[tot] = n;
		return f[tot];
	}
	int find(ll n) {
		int y = n % M;
		for(int p = fi[y]; p; p = nt[p])
			if(h[p] == n) return 1;
		return 0;
	}
} h;

ll dg(ll n) {
	if(n <= 1e6) return s[n];
	if(h.find(n)) return h[n];
	ll s = 1;
	for(ll i = 2,j; i <= n; i = j + 1) {
		j = n / (n / i);
		s = (s - (cn1(j) - cn1(i - 1)) * dg(n / i)) % mod;
	}
	s = (s % mod + mod) % mod;
	return h[n] = s;
}

int main() 
{
	init(1e6);
	scanf("%lld",&n);
	ll ans = 0;
	for(ll i = 1,j; i <= n; i = j + 1) {
		j = n / (n / i);
		ans = (ans + (dg(j) - dg(i - 1)) * calc(n / i))	% mod;
	}
	ans = (ans % mod + mod) % mod;
	pp("%lld\n",ans);
}