题解 [校内测试]图森破

题意

设 $f(n)$ 为满足以下条件的字符串个数：

1. 串长为 $S$，字符集为 $[0,9]$。
2. 设 $su{f}_i$ 为第 $i$ 个字符对应的后缀在最后补 $0$ 补成长度为 $n$ 得到的串，则 $\forall i\in [2,n],su{f}_i>su{f}_1$

试求 ${\sum }_{i=1}^n{i}^{2}f(i)$ 模 $258280327$ 得到的值。

组合

试求 $f(n)$。

首先考虑长度为 $n$ 的答案是什么。

首先所有循环节小于 $n$ 的串是不合法的，这十分显然。

然后对于一个不具有循环节的串，我们可以转它（${\text{newS}}_{i\%n+1}\leftarrow S_i$），可以转出 $n$ 个互不相同的串。下面我们证明这 $n$ 个串中有且仅有一个是合法的。

1. 这 $n$ 个串中至少有一个是合法的：

考虑设 $S$ 为 $n$ 个串中字典序最小的串，若它不合法，则存在一个后缀 $B<S$。由于 $S$ 是字典序最小的，因此 $B$ 的长度不超过 $\lfloor \frac{n}{2}\rfloor$。

由于 $S$ 是字典序最小的，因此如果设 $S$ 去掉 $B$ 为 $A$，则 $B$ 是 $A$ 的前缀。

然后考虑转到 $B$ 的字符串。

由于 $S$ 是字典序最小的，因此对应字符存在如下关系：

此时红色后缀循环移位后得到

2. 这 $n$ 个串中至多有一个是合法的：

假设有两个串合法，设为 $S_1,S_2$，考虑 $S_2$ 设 $S_2$ 在 $S_1$ 中的位置的后缀是 $F_2$，另一个是 $F_1$，知 $S_2\ge F_2,S_1\ge F_1$（补 0），而 $S_2<F_1,S_1<F_2$ （定义），故 $S_2\ge F_2>S_1\ge F_1$，与 $S_2<F_1$ 矛盾！

因此我们只需要算循环节等于 $n$ 的串的个数，设为 $g(n)$。

事实上这非常好算啊！易知 ${\sum }_{d|n}g(d)=10^n$，所以 $g(n)={\sum }_{d|n}10^{\frac{n}{d}}\times \mu (d)$ （莫比乌斯反演）！

$$f(n)=\frac{1}{n}\sum _{d|n}10^{\frac{n}{d}}\times \mu (d)$$

数论

$$\begin{array}{rl}& \sum _{n=1}^{N}n\sum _{d|n}\mu (\frac{n}{d})10^d\\ & =\sum _{d=1}^{N}d10^d\sum _{n=1}^{N/d}n\mu (n)\end{array}$$

注意到这是杜教筛简单题，因此 $\Theta (n^{2/3}\log n)$ 求解。（可以 $\Theta (n^{2/3})$，只需要把 map 和算 sum(i¹⁰i) 换掉，前者可以只记忆化 i 而不是 N/i，后者可以分块预处理来 O(1) 算 10^n）

代码

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

typedef long long ll;

const int NN = 5000000;
const int MOD = 258280327;
int mu[NN + 5], P[NN + 5], pl;
ll sum[NN + 5], ans, N;
bool np[NN + 5];
map<ll, ll> M;

ll power(ll a, ll b, ll p = MOD) {
	ll ret = 1; a %= p;
	for (; b; a = a * a % p, b >>= 1) if (b & 1) ret = ret * a % p;
	return ret;
}

const ll inv2 = power(2, MOD - 2), inv9 = power(9, MOD - 2);

ll getS(ll n) {
	if (n <= NN) return sum[n];
	if (M.find(n) != M.end()) return M[n];
	ll s = 1;
	for (ll i = 2, j; i <= n; i = j + 1) {
		j = n / (n / i);
		s = (s - getS(n / i) * (i + j) % MOD * (j - i + 1) % MOD * inv2) % MOD;
	} M[n] = s; return s;
}

inline ll S(ll n) {
	ll pw = power(10, n + 1); n %= MOD;
	return (pw * n - (pw - 10) * inv9) % MOD * inv9 % MOD;
}
int main() {
	scanf("%lld", &N), mu[1] = 1;
	for (int i = 2; i <= NN; ++i) {
		if (!np[i]) P[++pl] = i, mu[i] = -1;
		for (int j = 1; j <= pl && i * P[j] <= NN; ++j) {
			np[i * P[j]] = 1;
			if (i % P[j]) mu[i * P[j]] = mu[i] * mu[P[j]];
			else { mu[i * P[j]] = 0; break; }
		}
	}
	for (int i = 1; i <= NN; ++i) sum[i] = (sum[i - 1] + i * mu[i]) % MOD;
	getS(N);
	for (ll i = 1, j; i <= N; i = j + 1) {
		 j = N / (N / i);
		 ans = (ans + getS(N / i) * (S(j) - S(i - 1))) % MOD;
	}
	printf("%lld\n", (ans + MOD) % MOD);
	return 0;
}