【题解】bzoj-3884 上帝与集合的正确做法

模指数塔问题解析

最新推荐文章于 2019-04-02 18:11:29 发布

原创最新推荐文章于 2019-04-02 18:11:29 发布 · 270 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#欧拉定理 #欧拉函数

杂题集同时被 2 个专栏收录

48 篇文章

订阅专栏

--------------数学--------------

24 篇文章

订阅专栏

本文探讨了一种特殊的模指数塔问题，即求解形如2^(2^(2^(...))) mod p的问题。通过引入欧拉函数及其性质，文章提供了一种有效的递归求解策略，并给出了具体的实现代码。

Problem

to bzoj

求值：

2 2 2 2 . . . mod p

$2^{2^{2^{2^{...}}}} \bmod p$

Solution

看到无限的 $2$ 就很 $2$ 啊，最初想法是感觉是找循环节，觉得应该会稳定在一个数

看了skywalkert dalao的博客蒟蒻才渐渐顿悟

$2^{2^{2^{2^{...}}}} \pmod p=\frac{2^{2^{2^{2^{...}}}}}{{2^{\phi (p)*t}}}\pmod p=2^{2^{2^{2^{...}}}\pmod {\phi (p)}} \pmod p$

这样在指数的位置 $2^{2^{2^{...}}}\bmod \phi (p)$ 又和原式一样了

只有 $p\rightarrow \phi p$ ，当 $p=0$ 时，答案为 $0$

这样递归求解即可

但

这样要求 $2$ 与 $p$ 互质，如果 $p$ 为偶数，算法就不成立

那就把 $p$ 中的 $2$ 因子拆分：令 $p=2^k·d$

则 $2^{2^{2^{2^{...}}}} \pmod {2^k·d}=2^k·2^{2^{2^{2^{...}}}-k} \pmod {\phi(d)}$

$2^k$ 放在外面，再递归求解

2 2 2 2 . . . - k (mod ϕ (d)) = (2 2 2 2 . . . mod ϕ (d) - k mod ϕ (d)) mod ϕ (d)

$2^{2^{2^{2^{...}}}}-k \pmod {\phi (d)}=(2^{2^{2^{2^{...}}}} \bmod {\phi(d)}-k \bmod {\phi(d)})\bmod {\phi(d)}$

其中 $2^{2^{2^{2^{...}}}} \bmod {\phi(d)}$ 可以递归求解，直到 $p=1\Leftrightarrow ans=0$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rg register

template <typename _Tp> inline void read(_Tp&x){
    rg char c11=getchar(),ob=0;x=0;
    while(c11^'-'&&!isdigit(c11))c11=getchar();if(c11=='-')c11=getchar(),ob=1;
    while(isdigit(c11))x=x*10+c11-'0',c11=getchar();if(ob)x=-x;return ;
}

inline int phi(int x){
    int res=x;
    for(rg int i=2;i*i<=x;++i)
        if(x%i==0){
            res=res/i*(i-1);
            while(x%i==0)x/=i;
        }
    if(x^1)res=res/x*(x-1);
    return res;
}

int qpow(int A,int B,int p){
    int ans=1;
    while(B){
        if(B&1)ans=(ll)ans*A%p;
        A=(ll)A*A%p;
        B>>=1;
    }
    return ans;
}

int work(int p){
    if(p==1)return 0;
    rg int tmp=0,ans=0;
    while(!(p&1))++tmp,p>>=1;
    rg int phip=phi(p);
    ans=(ll)(work(phip)+phip-tmp%phip)%phip;
    ans=qpow(2,ans,p);
    return ans<<tmp;
}

int main(){
    rg int T;rg int p;
    read(T);
    while(T--)read(p),printf("%d\n",work(p)%p);
    return 0;
}