[ZJOJ6271]整除(中国剩余定理+原根)

题意

[题目描述]

整除符号为$|$, $d|n$在计算机语言中可被描述为$n\%d==0$。
现有一算式$n|x^m-x$,给定$n$,$m$,求$[1,n]$以内$x$解的个数。
解可能很大，输出取模$998244353$。

[输入格式]

从文件$division.in$中读入数据。

其中$n$的给定方式是由$c$个不超过$t$的质数的乘积给出的，$c$和$t$的范围会在数据范围中给出。
第一行一个$id$表示这个数据点的标号。
多组数据，其中第二行一个整数$T$表示数据组数。
对于每一组数据：
第一行两个整数$c$和$m$。
第二行$c$个整数，这些整数都是质数，且两两不同，他们的乘积即为$n$。
由于你可以通过输入求出$t$,输入不再给出。

[输出格式]

输出到文件$division.out$中。

对于每组数据输出一行，表示解的个数。

[样例1输入]

0
1
2 3
2 3

[样例1输出]

6

[数据范围与提示]

测试点编号	$c\le$	$t\le$	$m\le$	$T\le$
1	2	$10^3$	$2$	50
2	2	$10^3$	$10^9$	50
3	2	$10^2$	$2$	10000
4	1	$10^4$	$2$	50
5	2	$10^4$	$2$	50
6-8	10	$10^4$	$10^9$	50
9-10	50	$10^4$	$10^9$	50

题解

$x^{m-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}n)$的解的个数,在[1,n-1]范围内就是$gcd(m-1,n-1)$。
原问题求$x^m\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}{p}_1{p}_2...p_c)$
可以将其化为
$$\{\begin{array}{l}{x}^{m-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_1)\\ x^{m-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_2)\\ ...\\ x^{m-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_c)\end{array}$$
注意到从每个方程的一个合法解中任取一个乘起来，就能够组成一个原方程合法的解

因此我们就只用算每个方程的解就好了…

而单个方程的解的个数就为$\gcd (m-1,p-1)+1$。

乘起来即可。

证明:

我们要求满足$x\in [1,p]$时 $x^m\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的 $x$ 的个数。

首先 $x=p$ 一定满足，故不特殊考虑，最后答案 $+1$ 即可。接下来只考虑 $x\in [1,p-1]$ 的情况。由于 $p$ 为质数，因此$p$存在一个原根 $g$，$[1,p-1]$内的任意数在 $\mathrm{m}\mathrm{o}\mathrm{d}p$ 意义下都可以表示为 $g^y$ 的形式。这样，原方程就转化为：

$$g^{my}\equiv g^y(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

根据费马小定理：

$$my\equiv y(\mathrm{m}\mathrm{o}\mathrm{d}p-1)\Rightarrow (m-1)y\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p-1)$$

令 $k=\gcd (m-1,p-1)$，两边同时除以 kk，得：

$$\frac{m-1}{k}y\equiv 0(mod\frac{p-1}{k})$$

由于此时 $gcd(\frac{m-1}{k},\frac{p-1}{k})=1$，因此 $y$ 一定有 $\frac{p-1}{k}|y$。由于 $y\in [0,p-2]$，显然，$\frac{p-1}{k}$ 的任意小于 $k$ 的非负整数倍（$0\sim k-1$ 倍）均满足条件，因此 $y$ 共有 $k$ 种合法取值。

因此满足 ${\sum }_{x=1}^p[x^m\equiv x\mathrm{m}\mathrm{o}\mathrm{d}p]$ 的 $i$ 共有 $k+1$个。

[来自ImagineC大佬]

实现

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define MAXN 10000000
#define LL long long
#define MOD 998244353
int T,n,m;
int gcd(int a,int b){
	return (!b)?a:gcd(b,a%b);
}
int main()
{
	freopen("division.in","r",stdin);
	freopen("division.out","w",stdout);
	scanf("%*d%d",&T);
	while(T--){
		scanf("%d%d",&n,&m);
		LL ans=1;
		for(int i=1;i<=n;i++){
			int x;scanf("%d",&x);
			ans=ans*(gcd(x-1,m-1)+1)%MOD;
		}printf("%lld\n",ans);
	}
}