【bzoj3884】上帝与集合的正确用法 【欧拉函数】

本文介绍了一种利用欧拉函数的性质解决特定数学问题的方法。通过递归计算,解决了形如222...模p的问题,并给出了完整的C++实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目传送门
题解
这道题其实就是一个欧拉函数的性质的应用:
(a,m)=1(a,m)=1a,m2a,m≥2,则aφ(m) mod m1aφ(m) mod m≡1
我们设x=222...x=222...p=2kqp=2kq(q,2)=1(q,2)=1
x mod px mod p
=2k(2xk mod q) mod p=2k(2x−k mod q) mod p
=2k(2(xk) mod φ(q)+φ(q) mod q) mod p=2k(2(x−k) mod φ(q)+φ(q) mod q) mod p
我们令f(i,j)=(xi) mod jf(i,j)=(x−i) mod j
f(0,p)=2k(2f(k,φ(q))+φ(q) mod q) mod pf(0,p)=2k(2f(k,φ(q))+φ(q) mod q) mod p
于是我们只需要递归计算即可。当q=1q=1就返回00
答案就是f(0,p)f(0,p)
代码

#include<cstdio>
#define int long long
int t,p;
int calcphi(int n){
    if(n==1||n==2){
        return 1;
    }
    int res=n;
    for(int i=2;i*i<=n;i++){
        if(n%i==0){
            res=res/i*(i-1);
            while(n%i==0){
                n/=i;
            }
        }
    }
    if(n!=1){
        res=res/n*(n-1);
    }
    return res;
}
int fastpow(int a,int x){
    int res=1;
    while(x){
        if(x&1){
            res*=a;
        }
        x>>=1;
        a*=a;
    }
    return res;
}
int fastpow2(int a,int x,int mod){
    int res=1;
    while(x){
        if(x&1){
            res=res*a%mod;
        }
        x>>=1;
        a=a*a%mod;
    }
    return res;
}
int solve(int i,int j){
    if(j==1){
        return 0;
    }
    int u=0,tmp=j;
    while(!(tmp&1)){
        u++;
        tmp>>=1;
    }
    int phi=calcphi(tmp);
    return (fastpow(2,u)*(fastpow2(2,solve(u,phi)+phi,tmp)%tmp)%j-i%j+j)%j;
}
signed main(){
    scanf("%lld",&t);
    while(t--){
        scanf("%lld",&p);
        printf("%lld\n",solve(0,p));
    }
    return 0;
}

转载于:https://www.cnblogs.com/2016gdgzoi471/p/9476855.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值