NOI 2016 循环之美 题解

题目传送门

题目大意： 问有多少个值不相同的分数 $\frac{x}{y}(x\in [1,n],y\in [1,m])$，满足在 $k$ 进制下是个纯循环小数。纯循环：即小数存在循环节（可以是 $0$），且从小数点后第一位就开始循环。

题解

假设 $k$ 为 $10$，循环节为 $x_1{x}_2...x_c$，那么小数点部分就可以表示为这样的分数：$\frac{x_1{x}_2...x_c}{999...99}$，下面是 $c$ 个 $9$。

也就是说，对于分数 $\frac{x}{y}$，只有当 $y$ 的某个倍数是若干个 $9$ 时，$\frac{x}{y}$ 才是纯循环小数。

$c$ 个 $9$ 组成的数可以表示成 $10^c-1$，也就是要存在一个 $c$ 满足 $y|(10^c-1)$，这个东西可以转化成 $10^c\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}y)$，根据扩欧，只有当 $10$ 和 $y$ 互质时，$10$ 才存在逆元，也就是存在满足条件的 $c$。

类似的，我们可以得到，只有当 $y\perp k$ 时，$\frac{x}{y}$ 才是纯循环小数。

那么此时已经可以写出答案式子了：
$$\sum _{i=1}^n\sum _{j=1}^m[i\perp j][j\perp k]$$

对 $[i\perp j]$ 使用莫反，得到：
$$\begin{array}{rl}& =\sum _{i=1}^n\sum _{j=1}^m[j\perp k]\sum _{d|\gcd (i,j)}\mu (d)\\ & =\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{i=1}^n[d|i]\sum _{j=1}^m[d|j][j\perp k]\\ & =\sum _{d=1}^{\min (n,m)}\mu (d)[d\perp k]\lfloor \frac{n}{d}\rfloor \sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }[j\perp k]\end{array}$$

一开始尝试往另外两个互质上套莫反，然而发现会越推越丑……所以直接考虑从这个式子入手求解。

注意到里面有一个 $\lfloor \frac{n}{d}\rfloor$，这提示我们进行整除分块，假如能对另外两个部分进行快速的前缀求和，那么就可以套整除分块了。

设 $g(n,k)={\sum }_{i=1}^n\mu (i)[i\perp k],f(n,k)={\sum }_{i=1}^n[i\perp k]$，注意到如果 $k$ 有 $>1$ 个质因子 $p_i$，那么只保留一个其他全部除掉是不会影响这两个式子的，分析一下互质的性质就能发现，于是我们预处理一下这个 $k$，现在 $k$ 满足不包含平方因子。

$f,g$ 没有什么很好的直接求的方法，所以只能尝试递推。对于 $k$ 的一个质因子 $p$，考虑除掉 $p$ 之后答案的变化，先看 $f$：
$$\begin{array}{rl}f(n,k)& =\sum _{i=1}^n[i\perp k]\\ & =\sum _{i=1}^n[i\perp \frac{k}{p}]-\sum _{i=1}^n[p|i][i\perp \frac{k}{p}]\\ & =f(n,\frac{k}{p})-\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }[ip\perp \frac{k}{p}]\\ & =f(n,\frac{k}{p})-\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }[i\perp \frac{k}{p}]\\ & =f(n,\frac{k}{p})-f(\lfloor \frac{n}{p}\rfloor ,\frac{k}{p})\end{array}$$

中间 $[ip\perp \frac{k}{p}]=[i\perp \frac{k}{p}]$ 这一步就用到了k没有平方因子这个性质。

$g$ 的递推是类似的：
$$\begin{array}{rl}g(n,k)& =\sum _{i=1}^n\mu (i)[i\perp k]\\ & =\sum _{i=1}^n\mu (i)[i\perp \frac{k}{p}]-\sum _{i=1}^n\mu (i)[p|i][i\perp \frac{k}{p}]\\ & =g(n,\frac{k}{p})-\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\mu (ip)[ip\perp \frac{k}{p}]\end{array}$$

$\mu$ 有一个性质，$\mu (ij)=\mu (i)\mu (j)[i\perp j]$，证明的话，若不满足 $i\perp j$，那么 $ij$ 就含有平方因子，则 $\mu (ij)=0$；如果满足那么根据积性函数定义就显然可以拆开。

代入得到：
$$\begin{array}{rl}& =g(n,\frac{k}{p})-\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\mu (i)\mu (p)[i\perp p][i\perp \frac{k}{p}]\\ & =g(n,\frac{k}{p})+\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\mu (i)[i\perp k]\\ & =g(n,\frac{k}{p})+g(\lfloor \frac{n}{p}\rfloor ,k)\end{array}$$

$f,g$ 的个数是 $\sqrt{n}\times \frac{\log k}{\log \log k}$ 的，记忆化一下递推就好，当 $k=1$ 时，$f(n,k)$ 显然为 $n$，$g(n,k)$ 就是 $\mu$ 的前 $n$ 项之和，弄一个杜教筛就好。

然而存 $f,g$ 时我写的哈希冲突很多所以跑得飞慢，最后直接换了map反而来的快……代码如下：

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

int n,m,k;
int prime[maxn],t=0,mu[maxn];
bool v[maxn];
void SieveInit(){
	mu[1]=1;
	for(int i=2;i<=maxn-10;i++){
		if(!v[i])prime[++t]=i,mu[i]=-1;
		for(int j=1;j<=t&&i*prime[j]<=maxn-10;j++){
			v[i*prime[j]]=true;
			if(i%prime[j]==0)break;
			mu[i*prime[j]]=-mu[i];
		}
		mu[i]+=mu[i-1];
	}
}
map<int,int> Smu;
struct par{
	int x,y;
	bool operator <(const par &B)const{return x==B.x?y<B.y:x<B.x;}
};
map<par,int> f,g;
int calc_mu(int x){
	if(x<=maxn-10)return mu[x];
	if(Smu.count(x))return Smu[x];
	long long tot=1;
	for(int l=2,r;l<=x;l=r+1){
		r=x/(x/l);
		tot-=calc_mu(x/l)*(r-l+1);
	}
	Smu[x]=tot;
	return tot;
}
int fac[maxn];
int calc_f(int x,int y){
	if(x==0)return 0;
	if(y==1)return x;
	par p=(par){x,y};int re;
	if(f.count(p))return f[p];
	re=calc_f(x,y/fac[y])-calc_f(x/fac[y],y/fac[y]);
	f[p]=re;
	return re;
}
int calc_g(int x,int y){
	if(x==0)return 0;
	if(y==1)return calc_mu(x);
	par p=(par){x,y};int re;
	if(g.count(p))return g[p];
	re=calc_g(x,y/fac[y])+calc_g(x/fac[y],y);
	g[p]=re;
	return re;
}

int main()
{
	scanf("%d %d %d",&n,&m,&k);SieveInit();
	for(int i=2;i<=k;i++)//去掉多余的质因子
		while(k%(i*i)==0)k/=i;
	int k_=k;for(int i=2;i<=k;i++)
		if(k_%i==0)fac[k_]=i,k_/=i;
	long long ans=0;
	for(int l=1,r;l<=min(n,m);l=r+1){
		r=min(n/(n/l),m/(m/l));
		ans+=1ll*(calc_g(r,k)-calc_g(l-1,k))*(n/l)*calc_f(m/l,k);
	}
	printf("%lld",ans);
}