同余
模运算
一般的,对于任意整数 \(a,b\),\(b>0\)。则有:
性质:
- \(\frac{a}{k} \bmod n=\frac{a\bmod kn}{k}\)。
简证:原式等价于 \(\frac{a}{k}-\lfloor \frac{a}{kn} \rfloor \times n=(a-\lfloor \frac{a}{kn} \rfloor \times kn)\div k\).
等号右边乘法分配律拆分,命题显然成立。
证毕。
- \((a+b) \bmod n=(a\bmod n+b\bmod n)\bmod n\)
证明:
原式可写为 \(a+b-\lfloor \frac{a+b}{n} \rfloor \times n=a-\lfloor \frac{a}{n} \rfloor \times n+b-\lfloor \frac{b}{n} \rfloor \times n-\lfloor \frac{a-\lfloor \frac{a}{n} \rfloor \times n+b-\lfloor \frac{b}{n} \rfloor \times n}{n} \rfloor \times n\).把下取整号内整数提出:\(a+b-\lfloor \frac{a+b}{n} \rfloor \times n=a-\lfloor \frac{a}{n} \rfloor \times n+b-\lfloor \frac{b}{n} \rfloor \times n-\lfloor \frac{a+b}{n} \rfloor \times n+\lfloor \frac{a}{n} \rfloor \times n+\lfloor \frac{b}{n} \rfloor \times n\).
消掉一部分:\(a+b-\lfloor \frac{a+b}{n} \rfloor \times n=a+b-\lfloor \frac{a+b}{n} \rfloor \times n\).
显然,等式成立。
证毕。
- \(a\cdot b \bmod n=(a \bmod n)\times (b\bmod n)\bmod n\)。
证明与上面一条差不多。
- 若 \(a\equiv b\pmod p,c\equiv d\pmod p\),则有 \(a\cdot c\equiv b\cdot d\pmod p\)。
证明:一定有 \(a+k_1\times p=b,c+k_2\times p=d\)。
所以 \(b\times d=(a+k_1\times p)\times (c+k_2\times p)=a\times c+k_1pc+k_2pd+k_1k_2p^2\)。
而在模 \(p\) 意义下后面带 \(p\) 的项为 \(0\),故得证。
同余
若 \(a\bmod n=b\bmod n\),则称\(a\) 与 \(b\) 在模 \(n\) 意义下同余,写作 \(a \equiv b \pmod n\)。
定理:
-
\(a\equiv b\pmod n\iff (a-b)\equiv 0\pmod n\iff n\mid (a-b)\)。
-
若 \(a\equiv b\pmod n\),那么 \(a+c\equiv b+c\pmod n\),\(a\cdot c \equiv b\cdot c\pmod n\)。
同余类与剩余系
对于 \(a\in [0,m-1]\),集合 \(\{a+km\}(k\in \mathbb{Z})\) 中的所有数在模 \(m\) 意义下同余,余数均为 \(a\),则称该集合为一个模 \(m\) 的同余类,简记 \(\bar{a}\)。
\(m\) 的同余类一共有 \(\bar{0},\bar{1},\bar{2}\dots \bar{m-1}\),共 \(m\) 个,它们组成 \(m\) 的完全剩余系。
其中余数与 \(m\) 互质的一共有 \(\varphi(m)\) 个,其组成 \(m\) 的简化剩余系。
简化剩余系关于模 \(m\) 乘法封闭。
定理:
费马小定理:若 \(p\) 为质数,则 \(a^p\equiv a\pmod p\)。
欧拉定理:若 \(\gcd(a,n)=1\),则 \(a^{\varphi(n)}\equiv 1\pmod n\)。
证明:设 \(n\) 的简化剩余系为 \(\{\bar{s_1},\bar{s_2}\dots \bar{s_{\varphi(n)}}\}\)。
对于任意 \(x,y\),若 \(a\cdot x\equiv a\cdot y\pmod n\),故 \(a\cdot(x-y)\equiv 0\pmod n\)。
故 \((x-y)\equiv 0\pmod n\),即 \(x\equiv y\pmod n\)。
所以若 \(a\cdot x\) 与 \(a\cdot y\) 出自相同同余类,则 \(x,y\) 也出自相同同余类中,反之亦然。
由于简化剩余系运算关于模 \(m\) 乘法封闭,故 \(\bar{a\cdot s_i}\) 也属于简化剩余系。
所以集合 \(A=\{\bar{s_1\cdot a},\bar{s_2\cdot a}\dots \bar{s_{\varphi(n)}\cdot a}\}\) 与 \(B=\{\bar{s_1},\bar{s_2}\dots \bar{s_{\varphi(n)}}\}\) 在模 \(n\) 意义下相等。
故 \(\prod _{s\in A} s\equiv a^{\varphi(n)}\times \prod_{s\in B} \equiv \prod_{s\in B}\pmod n\)。
即 \(a^{\varphi(n)}\equiv 1\pmod n\)。
欧拉定理证毕。
若 \(n\) 为质数,则 \(\varphi(n)=n-1\),代入原式费马小定理成立。
若 \(\gcd(a,n)>1\),那么 \(a\) 只能是 \(n\) 的倍数,显然成立。
费马小定理证毕。
欧拉定理推论
推论1(降幂公式):
\(\forall \gcd(a,n)=1,a^b\equiv a^{b\bmod \varphi(n)}\pmod n\).
证明:若 \(b<\varphi(n)\) 显然成立。
否则,设 \(b=k\times \varphi(n)+r\),其中 \(r=b\bmod \varphi(n)\)。
所以:
\[\begin{aligned} a^b&\equiv a^{k\times \varphi(n)+r}\\ &\equiv (a^{\varphi(n)})^k\times a^r\\&\equiv 1^k\times a^r\\&\equiv a^{b\bmod \varphi(n)} \end{aligned} \]证毕。
推论2:
若 \(\gcd(a,n)=1\),则使得 \(a^{x}\equiv 1\pmod n\) 成立的最小正整数 \(x\) 满足 \(x\mid \varphi(n)\)。
证明:
反证法,假设最小正整数 \(x\) 满足 \(x\nmid \varphi(n)\)。
设 \(\varphi(n)=qx+r(0<r<x)\),因为 \(a^x\equiv 1\pmod n\),所以 \(a^{qx}\equiv 1\pmod n\),根据欧拉定理,\(a^{qx+r}\equiv 1\pmod n\),故 \(a^r\equiv 1\pmod n\)。
因为 \(r<x\),所以 \(r\) 显然是满足条件的更小的正整数,与假设矛盾,故假设不成立,原命题成立。
证毕。
扩展欧拉定理
若 \(a,m\in \mathbb{z}\),则有:
证明:
前者显然成立,主要证明 \(b\ge \varphi(m)\) 的情况。
若 \(\gcd(a,m)=1\),则根据推论 \(1\) 可轻松得证。
否则,设 \(m\) 的质因数分解为 \(p_1^{c_1}p_2^{c_2}p_3^{c_3}\dots p_k^{c_k}\)。
引理:\(\forall i,x \equiv y\pmod {p_i^{c_i}}\) 成立是 \(x\equiv y\pmod m\) 成立的充要条件。
证明:显然 \(\forall i\ne j,\gcd(p_i^{c_i},p_j^{c_j})=1\) 恒成立。
考虑对于两个同余式 \(x\equiv y\pmod n,x\equiv y\pmod m\),若 \(\gcd(n,m)=1\),则有 \(x\equiv y \pmod {nm}\)。
证明很简单,移项有 \(n\mid (x-y),m\mid (x-y)\),所以有 \(\operatorname{lcm}(n,m)\mid (x-y)\) 即 \((nm)\mid (x-y)\),再次移项即可证明,必要性也可以这么证。
通过上述定理,显然可以将若干模数两两互质且两个同余号左右内容一致的同余方程合并成一个。故引理成立。
证毕。
若 \(\gcd(a,p_i^{c_i})=1\),则根据降幂公式,有 \(a^b\equiv a^{b\bmod \varphi(m)+\varphi(m)}\pmod m\),所以 \(a^b\equiv a^{b\bmod \varphi(m)+\varphi(m)}\pmod {p_i^{c_i}}\)
若 \(\gcd(a,p_i^{c_i})\ne 1\),则一定有 \(p_i\mid a\),记 \(a=z\times p_i\)。
引理:\(b\ge \varphi(m)\ge \varphi(p_i^{c_i})\ge c_i\).
证明:显然前两个不等号成立,这里主要证明第三个。
因为 \(p_i\ge 2\),所以 \(p_i\) 越大时,\(\varphi(p_i^{c_i})=p_i^{c_i-1}\times (p_i-1)\) 越大,而 \(c_i\) 不变,故只考虑 \(p_i=2\) 时即可。
考虑分讨,\(c_i=1\) 时 \(2^0=1\) 成立,\(c_i=2\) 时 \(2^1=2\) 成立。
\(c_i>2\) 时,有 \(2^{c_i-1}>c_i\) 恒成立。综上引理成立,证毕。
所以 \(a^b=z^b\times p_i^b\),综合引理有 \((p_i^{c_i})\mid a^b\)。
同理,\((p_i^{c_i})\mid a^{b\bmod \varphi(m)+\varphi(m)}\)。
故 \(a^b\equiv a^{b\bmod \varphi(m)+\varphi(m)}\equiv 0\pmod {p_i^{c_i}}\)。
综上所述,扩展欧拉定理成立。
证毕。
例1.上帝与集合的正确用法
其实就是求 :
根据扩展欧拉定理,有:
所以只需求:
发现可以递归求。
引理:\(\varphi(\varphi(\varphi(\dots p \dots)))\) 收敛到 \(1\),至多在外面套上 \(O(\log p)\) 个 \(\varphi\)。
证明:首先根据 \(\varphi(n)\) 的计算式,\(\varphi(n)=\prod _{i=1}^m p_i^{c_i-1}\times (p_i-1)\),如果 \(n=1\) 或者 \(n=2\) 显然是奇数,
如果 \(\varphi(n)\) 是偶数,则一定有以下条件成立其中一个:
\(\exist i,p_i\) 是奇数。
\(n\) 的质因数分解中,质因子 \(2\) 的次幂 \(>1\)。
很显然,上述两个条件成立其中一个必然会让 \(\varphi(n)\) 乘上一个偶数,从而使得整个值一定为偶数。
显然,除了 \(1,2\) 其余所有正整数均满足上述两个条件。
所以 \(\forall x\in[2,+\infty)\cap\mathbb{Z}\),\(\varphi(x)\equiv 0\pmod 2\)。
而又因为如果 \(n\) 是偶数,则比 \(n\) 小的所有偶数都不与 \(n\) 互质,故 \(\varphi(n)\le \frac{n}{2}\)。
综合上述两点,命题成立,证毕。
const int N=1e7+10;
int t,Mo;
signed phi[N],v[N];
vector<signed>prime;
int power(int a,int b,int p) {
int res=1;
for(;b;b>>=1) {
if(b&1) res=res*a%p;
a=a*a%p;
}
return res%p;
}
void Euler(int n) {
FOR(i,1,n) phi[i]=i;
for(int i=2;i<=n;++i) {
if(!v[i]) v[i]=i,prime.pb(i),phi[i]=i-1;
for(int p:prime) {
if(p*i>n||v[i]<p) break;
v[i*p]=p;
phi[p*i]=phi[i]*(i%p?p-1:p);
}
}
return ;
}
int dfs(int x) {
if(x==1) return 0ll;
return power(2,dfs(phi[x])+phi[x],x);
}
void solve() {
cin>>Mo;
cout<<dfs(Mo)<<"\n";
}
main() {
cin>>t;
Euler(N-10);
while(t--) solve();
return 0;
}
模意义下的除法
由于除法在模意义下不一定满足 若 \(a\equiv b\pmod n\) 等价于 \(\frac{a}{c}\equiv \frac{b}{c} \pmod n\),所以要找到一种办法去求除法取模意义下的值。
逆元:一般的,对于一个数 \(x\),若存在 \(x'\) 使得同余方程 \(x'\cdot x\equiv 1\pmod n\) 成立,则称 \(x'\) 为 \(x\) 在模 \(n\) 意义下的逆元,记为 \(x^{-1}\)。
引理:\(x\) 存在模 \(n\) 意义下的逆元,当且仅当 \(\gcd(x,n)=1\)。
证明:反证法,假设 \(\gcd(x,n)\ne 1\)。设 \(d=\gcd(x,n),x\times x'=k_1\times d,n=k_2\times d\)
\(x\cdot x' \bmod n=x\cdot x'-\lfloor \frac{x\cdot x'}{n} \rfloor\times n=k_1\times d-\lfloor \frac{k_1}{k_2} \rfloor\times k_2\times d=(k_1-\lfloor \frac{k_1}{k_2}\rfloor\times k_2)\times d\)。
又因为 \(d>1\),且 \(x\times x'\bmod n\)是 \(d\) 的倍数,故 \(x\times x' \bmod n\) 不可能为 \(1\)。故假设不成立,原命题成立。
证毕。
费马小定理求逆元
若 \(n\) 为质数,则求 \(a\) 在模 \(n\) 意义下的逆元:
根据费马小定理,\(a^n\equiv a\pmod n\),所以 \(a^{n-1}\equiv 1\pmod n\)。
将同余号左写成 \(a\times k\) 的形式。显然,\(k=a^{n-2}\)。
所以 \(k=a^{n-2}\) 是方程 \(a\times l\equiv 1\pmod n\) 的一组解。
所以快速幂直接求 \(a^{n-2}\bmod n\) 即可,时间复杂度 \(O(\log n)\)。
欧拉定理求逆元
与费马小定理几乎一致。
因为要求解 \(n\) 的欧拉函数值,所以时间复杂度为 \(O(\sqrt n)\),显然不是很优秀。
但此方法可以求任意一对 \(a,n\) 的逆元(当然要满足 \(\gcd(a,n)=1\))。
线性递推 \(1\sim n\) 的所有逆元
第一种方法:
显然 \(1^{-1}=1\)。
考虑当前要求 \(n^{-1}\)。
显然有 \(\lfloor \frac{p}{n} \rfloor\times n+p\bmod n=p\)。
即 \(\lfloor \frac{p}{n} \rfloor\times n+p\bmod n \equiv 0 \pmod p\)。
令 \(\lfloor \frac{p}{n} \rfloor=a,p\bmod n=b\)。
同时乘上 \(n^{-1}\cdot b^{-1}\)。
而由于 \(p\bmod n\le n\),所以只要把 \(1\sim n-1\) 的逆元全部求出就一定可以求出 \(n\) 的逆元。
但是,当 \(p\) 为合数时或者 \(n=p\) 时,有可能 \(n\mid p\),此时 \(p\bmod n=0\),\(0^{-1}\) 一定不存在小学老师告诉过你 \(0\) 能做除数吗。
所以线性递推逆元的方法(包括下面的阶乘递推逆元)都只适用于 \(n\le p\) 且 \(p\) 为质数的情况。
第二种方法:
预处理出 \(1\sim n\) 所有数的阶乘在 \(\bmod p\) 意义下的值。
显然 \(n!^{-1}=(n+1)!^{-1}\times (n+1)\)。
这样就可以预处理出所有 \(1\sim n\) 阶乘的逆元。
由于 \(i^{-1}=i!^{-1}\times (i-1)!\)。
所以也可以线性求。
一般用第二种方法,因为第二种方法会顺便求出阶乘以及阶乘逆元,这在求组合数中十分有用处。
裴蜀定理(Bézout定理)
对于任意整数 \(a,b\),存在整数 \(x,y\),使得 \(ax+by=\gcd(a,b)\)。
证明:
在欧几里得求 \(\gcd\) 的最后一步,即 \(b=0\) 时,一定存在一组解 \(x=1,y\in \mathbb{Z}\),使得 \(ax+by=a+0=a=\gcd(a,b)\)。
假设存在一组解 \(x,y\) 使得 \(bx+(a\bmod b)y=\gcd(b,a\bmod b)\)。
根据欧几里得算法,\(\gcd(a,b)=\gcd(b,a\bmod b)\)。
又因为:
\[\begin{aligned} &bx+(a\bmod b)y\\ =&bx+(a-\lfloor \frac{a}{b} \rfloor\times b)y\\ =&bx+ay-\lfloor \frac{a}{b} \rfloor\times by\\ =&ay+b(x-\lfloor \frac{a}{b} \rfloor\times y) \end{aligned} \]则此时令 \(x'=y,y'=x-\lfloor \frac{a}{b} \rfloor\times y\),就得到了使得 \(ax'+by'=\gcd(a,b)\) 成立的一组解。
对此使用数学归纳法,命题得证。
该证明同时还给出了整数 \(x,y\) 的计算方法,我们称其为扩展欧几里得算法。
void Exgcd(int a,int b,int &x,int &y) {
if(!b) return x=1,y=0,void();
Exgcd(b,a%b,x,y);
int t=x;x=y;y=t-(a/b)*y;
}
对于更为一般的情况,即 \(ax+by=c\),该方程有解当且仅当 \(\gcd(a,b)\mid c\)。
证明:
设 \(d=\gcd(a,b),a=k_1\cdot d,b=k_2\cdot d\)。
则:
\[\begin{aligned} &ax+by=c\\ \Rightarrow &dk_1x+dk_2y=c\\ \Rightarrow &d(k_1x+k_2y)=c\\ \Rightarrow &k_1x+k_2y=\frac{c}{d} \end{aligned} \]显然左侧为整数,故右侧也应为整数,即 \(d\mid c\)。
证毕。
先用扩展欧几里得算法求出一组方程 \(ax+by=\gcd(a,b)\) 的特解,将求出的特解 \(x_0,y_0\) 同时乘上 \(\frac{c}{\gcd(a,b)}\) 即可求出原方程的一组特解。
原方程通解即:
证明:
- 当 \(a,b\) 互质。
则原式可写为 \(y=\frac{c-a\times x}{b}\)。
若要满足整数的条件,则必须满足 \(c\equiv a\times x \pmod b\)。
显然,\(c\) 不变,变的数字只有同余号右侧,假设改变后的数字是 \(a\times x+m\)。
而 \(x\) 改变只会让 \(a\times x\) 变成 \(a\) 的倍数,所以有 \(a\mid m\)。
又因为只有当 \(b\mid m\) 时,\(a\times x+m\) 在模 \(b\) 意义下的值等于 \(a\times x\)。所以有 \(b\mid m\)。
所以改变后的数字可以被写成 \(a\times x+k\times \operatorname{lcm}(a,b)=a\times x+k\times ab\)。
此时 \(x\) 增长了 \(k\times b\)。
同时解出对应的 \(y\),增长了 \(-k\times a\)。
而 \(k\) 可以取任意正整数。
- 当 \(a,b\) 不互质。
不互质则可以通过将原方程 \(a,b,c\) 同时除以 \(\gcd(a,b)\) 得到,解集显然不变。
证毕。
例2. [NOIP2012 提高组] 同余方程
考虑对原式进行变形:
换元,令 \(y=-k\):
这是什么?不定方程!
它有解当且仅当 \(\gcd(a,b)=1\)。
直接用扩展欧几里得算法求出一组特解,再通过取模运算求得最小正整数解。
时间复杂度 \(O(\log (a+b))\)。
LL a,b;
/*
ax=1(mod b)
b|ax-1
ax+by=1
*/
int exgcd(LL a,LL b,LL &x,LL &y) {
if(!b) return x=1,y=6+1,a;
int d=exgcd(b,a%b,x,y);
LL t=x;x=y;y=t-(a/b)*y;
return d;
}
main() {
cin>>a>>b;
LL x,y;int d=exgcd(a,b,x,y);
cout<<(x%b+b)%b<<"\n";
return 0;
}
中国剩余定理(CRT)与线性同余方程组
设 \(m_1,m_2,m_3\dots m_n\) 是两两互质的整数,设 \(m=\prod _{i=1}^n m_i,M_i=\frac{m}{m_i},t_i\) 为 \(M_i\) 在模 \(m_i\) 意义下的逆元(可取任意解)。
对于任意 \(n\) 个整数 \(a_1,a_2,a_3\dots a_n\) 方程组:
有整数解,解为 \(x=\left(\sum_{i=1}^n a_iM_it_i\right)+km,k\in \mathbb{Z}\)。
证明:
考虑 \(k=0\) 的特解 \(x_0\)。
由于 \(\forall k\ne i,a_iM_it_i\equiv 0\pmod{m_k}\),且 \(\forall i,a_iM_it_i\equiv a_i\pmod {m_i}\)。
故 \(\forall i,x_0\equiv a_i\pmod{m_i}\)。
又因为模数两两互质可以合并,故得证。
中国剩余定理给出了模数两两互质的方程组的一组通解。
一般的,使用中国剩余定理求解线性方程组的时间复杂度为 \(O(n \log V)\),其中 \(V\) 为值域。
(其实,即使模数不满足两两互质,我们也有方法求解方程组,并有着与上述同样优秀的时间复杂度)
例3.曹冲养猪
自信点,把题目中“假定”去掉。
直接中国剩余定理即可。
const int N=11;
int n,a[N],c[N];
int exgcd(int a,int b,int &x,int &y) {
if(!b) return x=1,y=6+1,a;
int d=exgcd(b,a%b,x,y);
int z=x;x=y;y=z-(a/b)*y;
return d;
}
int mul(int x,int y,int p) {
int ans=0;
for(;y;y>>=1) {
if(y&1) ans=(ans+x)%p;
x=(x+x)%p;
}
return ans%p;
}
int CRT() {
int M=1ll,ans=0;
FOR(i,1,n) M*=a[i];
FOR(i,1,n) {
int MI=M/a[i],b=a[i],x,y;
int d=exgcd(MI,b,x,y);
x=(x%a[i]+a[i])%a[i];
ans=(ans+mul(mul(x,c[i],M),MI,M))%M;
}
return (ans+M)%M;
}
main() {
cin>>n;
FOR(i,1,n) cin>>a[i]>>c[i];
cout<<CRT()<<"\n";
return 0;
}
更加一般的情况
现在我们考虑更加一般的情况。
例4.扩展中国剩余定理(EXCRT)
问题:求解上述线性方程组,但不保证 \(m_i\) 两两互质。
考虑使用数学归纳法,不断合并方程组,最后归为一个方程并求出答案。
先考虑将两个方程合并成一个:
转化:
故:
换元,令 \(c\gets a_2-a_1,x_0\gets k_1,y_0\gets -k_2\)。
显然可以扩欧,若没有解则原方程组无解。
求出 \(x_0,y_0\) 的一组特解,将其中一个代入原式即可求出一个同时满足两个方程的特解 \(x'\):
考虑求出通解,发现要满足原来的式子加减的数字只能同时是 \(m_1\) 和 \(m_2\) 的倍数。
覆盖的最大范围显然是最小公倍数。
综上,我们求得了一个合并后的方程组:
合并 \(n\) 次后即可求解。
const int N=1e5+10;
int mul(int a,int b,int p) {
int ans=0;
for(;b;b>>=1) {
if(b&1) ans=(ans+a)%p;
a=(a+a)%p;
}
return ans%p;
}
int exgcd(int a,int b,int &x,int &y) {
if(!b) return x=1,y=6+1,a;
int d=exgcd(b,a%b,x,y);
int z=x;x=y;y=z-(a/b)*y;
return d;
}
int n,r[N],m[N];
int EXCRT() {
int M=m[1],R=r[1],d,x,y;
FOR(i,2,n) {
int c=((r[i]-R)%m[i]+m[i])%m[i];
d=exgcd(M,m[i],x,y);
if(c%d) return -1;
x=mul(x,c/d,m[i]/d);
R+=x*M;
M/=d;M*=m[i];
R=(R%M+M)%M;
}
return (R%M+M)%M;
}
main() {
cin>>n;
FOR(i,1,n) scanf("%lld%lld",&m[i],&r[i]);
cout<<EXCRT()<<"\n";
return 0;
}
综合应用
例5.abc193_e Oversleeping
设公交循环了 \(a\) 次 “从 A 出发到 B,在 B 停 y 秒后从 B 出发到 A,在 A 停 y 秒”,此人重复睡去醒来的次数为 \(b\),公交车已经停在 B 地 \(m\) 秒,此人醒来 \(n\) 秒下车。
则有:
显然这是一个不定方程的形式。
枚举 \(n,m\),使用扩展欧几里得求出一组 \(a,b\) 的最小正整数解并带入原式更新答案即可。
复杂度 \(O(T\times nm \log(x+y+p+q))\)。
例6.[六省联考 2017] 相逢是问候
暴力显然不可做。
根据例1,一个数字最多更新 \(O(\log p)\) 次在模 \(p\) 意义下的值就不会变化了。
所以可以用一个势能分析线段树维护。
const int N=5e4+10,M=55;
int a[N],n,m,P,c,tot,use[N],pos[M],cp[M][N],cp1[M][N];
bool b[M][N],b1[M][N],flg;
/*
a i表示原数列。
use i表示 a i 做了多少次操作。
pos i表示p做 i 次欧拉函数的值。
cp ij表示 c^j 对 pos i 取Mod的值。 (b)
cp1 ij表示 c^(10000j) 对 pos i 取Mod的值。(b1)
*/
int Phi(int x) {
int res=x,cnt=sqrt(x);
FOR(i,2,cnt) {
if(x%i) continue;
while(x%i==0) x/=i;
res/=i;res*=(i-1);
}
if(x>1) res/=x,res*=(x-1);
return res;
}
struct SegmentTree {
int l,r,sum;
bool all;
}t[N<<2];
void Pre() {
FOR(i,0,tot) {
cp[i][0]=1;
FOR(j,1,10000) {
cp[i][j]=cp[i][j-1]*c;
if(cp[i][j]>=pos[i]) b[i][j]=1,cp[i][j]%=pos[i];
b[i][j]|=b[i][j-1];
}
}
FOR(i,0,tot) {
cp1[i][0]=1;
b1[i][1]=b[i][10000];
FOR(j,1,10000) {
cp1[i][j]=cp1[i][j-1]*cp[i][10000];
if(cp1[i][j]>=pos[i]) b1[i][j]=1,cp1[i][j]%=pos[i];
b1[i][j]|=b1[i][j-1];
}
}
}
int power(int p,int ind) {
int p1=p/10000,p2=p%10000;
int res=cp[ind][p2]*cp1[ind][p1];
if(res>=pos[ind]) flg=1,res%=pos[ind];
flg|=b[ind][p2]|b1[ind][p1];
return res;
}
int dfs(int x,int dep,int ind) {
flg=0;
if(dep==ind) {
if(x>=pos[ind]) flg=1,x%=pos[ind];
return x;
}
int y=dfs(x,dep+1,ind);
return power(y+flg*pos[dep+1],dep);
}
void push_up(int p) {
t[p].all=t[p<<1].all;t[p].all&=t[p<<1|1].all;
t[p].sum=(t[p<<1].sum+t[p<<1|1].sum)%P;
}
void Build(int p,int l,int r) {
t[p].l=l,t[p].r=r;
if(l==r) return (void)(t[p].sum=a[l],t[p].all=use[l]=0);
int mid=l+r>>1;
Build(p<<1,l,mid);
Build(p<<1|1,mid+1,r);
push_up(p);
}
void change(int p,int l,int r) {
if(t[p].all) return ;
if(t[p].l==t[p].r) {
use[t[p].l]++;
if(use[t[p].l]==tot) t[p].all=1;
t[p].sum=dfs(a[t[p].l],0,use[t[p].l]);
return ;
}
int mid=t[p].l+t[p].r>>1;
if(l<=mid) change(p<<1,l,r);
if(r>mid) change(p<<1|1,l,r);
push_up(p);
}
int ask(int p,int l,int r) {
if(l<=t[p].l&&t[p].r<=r) return t[p].sum;
int mid=t[p].l+t[p].r>>1,val=0;
if(l<=mid) val+=ask(p<<1,l,r);
if(r>mid) val+=ask(p<<1|1,l,r);
return val%P;
}
main() {
cin>>n>>m>>P>>c;
FOR(i,1,n) a[i]=read();
pos[0]=P;
while(pos[tot]!=1) ++tot,pos[tot]=Phi(pos[tot-1]);
pos[++tot]=1;
Pre();
Build(1,1,n);
while(m--) {
int op=read(),l=read(),r=read();
if(!op) change(1,l,r);
else cout<<ask(1,l,r)%P<<"\n";
}
return 0;
}