题目
Description
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]=10n−∑d|n and d≠n f[d]
直接统计10…的贡献,不难得到:
f[n]=∑d|n10d∗μ(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);
}