指数循环节练习

最新推荐文章于 2024-06-14 14:40:57 发布
最新推荐文章于 2024-06-14 14:40:57 发布
阅读量133 收藏
点赞数
分类专栏： 数论 # 基础数论 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/CK1513710764/article/details/124790890
版权

1.HDU4335

我们先将n^{n!}\equiv b(mod\ p)转化为n^{n!%\varphi (p)+\varphi (p)}\equiv b(mod\ p).

然后分情况讨论:

1.n! <= phi(p)时,我们直接枚举有多少n满足条件即可

2.n! > phi(p)时,注意到,当出现一个 n 使得n! % phi(p) == 0,那么对于i > n,一定都有 i! % phi(p) == 0,因为 i! 包括 n!.所以,当n! % phi(p) != 0时,我们依旧枚举有多少个 n 满足条件.

3.当 n! % phi(p) == 0时,式子变成n^{phi(p)} \equiv b(mod\ p),这部分其实是存在循环节的,并且长度为p,具体证明可以查阅网上资料,这里不再赘述.我们求出每个循环节有多少个n是满足条件的即可.

最后补充一点,这道题的题给数据很大,尤其是对于ull_max 的情况,我们要进行特判,其他的数据也要开 ull 而不只是 ll

#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <map>
#include <unordered_map>
#include <vector>
#include <stack>
#include <time.h>
#include <set>
#define IOS ios::sync_with_stdio(false), cin.tie(0)
#define PII pair<int, int>
#define PDI pair<double, int>
#define int long long
#define ll long long
#define ull unsigned long long
#define eps 1e-6
#define MOD 1000000007
#define debug(a) cout << #a << "=" << a << endl;
#define all(x) (x).begin(), (x).end()
#define mem(x, y) memset((x), (y), sizeof(x))
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
using namespace std;
ull b, p, m;
ull get_phi(int n) {
    ull ans = n;
    for (int i = 2; i <= n / i; ++i) {
        if (n % i == 0) {
            ans = ans / i * (i - 1);
            while (n % i == 0)
                n /= i;
        }
    }
    if (n > 1)
        ans = ans / n * (n - 1);
    return ans;
}
ull qpow(ull a, ull k, ull p) {
    ull res = 1;
    while (k) {
        if (k & 1)
            res = res * a % p;
        k >>= 1;
        a = a * a % p;
    }
    return res;
}
signed main() {
#ifdef LOCAL
    freopen("input.in", "r", stdin);
    freopen("output.out", "w", stdout);
#endif
    IOS;
    int T, Case = 0;
    cin >> T;
    while (T--) {
        cout << "Case #" << ++Case << ": ";
        cin >> b >> p >> m;
        if (p == 1) {
            if (m == 18446744073709551615ull)
                cout << "18446744073709551616\n";
            else
                cout << m + 1 << "\n";
            continue;
        }
        //1.i! <= phi(p)
        int i = 0, phi = get_phi(p);
        ull fac = 1, ans = 0;
        for (; i <= m && fac <= phi; ++i) {
            if (qpow(i, fac, p) == b)
                ++ans;
            fac *= i + 1;
        }
        //2.i! > phi(p) 且 i! % phi(p) != 0
        fac %= phi;
        for (; i <= m && fac; ++i) {
            if (qpow(i, fac + phi, p) == b)
                ++ans;
            fac = (fac * (i + 1)) % phi;
        }
        //3.i! % phi(p) = 0
        if (i <= m) {
            ull cnt = 0;//表示一个循环节中有多少个b
            vector<int> ring(p);
            for (int j = 0; j < p; ++j) {
                ring[j] = qpow(i + j, phi, p);
                if (ring[j] == b)
                    ++cnt;
            }
            ans += (ull)(m - i + 1) / p * cnt;
            ull remain = (m - i + 1) % p;
            for (int j = 0; j < remain; ++j)
                if (ring[j] == b)
                    ++ans;
        }
        cout << ans << "\n";
    }
    
}

2.hdu2837

 这道题我们直接使用递归并搭配欧拉降幂即可,比较奇怪的是这道题的边界判断(题主也不知道为什么会这样).

#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <map>
#include <unordered_map>
#include <vector>
#include <stack>
#include <time.h>
#include <set>
#define IOS ios::sync_with_stdio(false), cin.tie(0)
#define PII pair<int, int>
#define PDI pair<double, int>
#define int long long
#define ll long long
#define ull unsigned long long
#define eps 1e-6
#define MOD 1000000007
#define debug(a) cout << #a << "=" << a << endl;
#define all(x) (x).begin(), (x).end()
#define mem(x, y) memset((x), (y), sizeof(x))
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
using namespace std;
int n, m;
int qpow(int a, int k, int p) {
    if (a == 0)
        return k == 0;
    int res = 1;
    while (k) {
        if (k & 1)
            res = res * a % p;
        k >>= 1;
        a = a * a % p;
    }
    return res == 0 ? p : res;
}
int phi(int x) {
    int ans = x;
    for (int i = 2; i <= x / i; ++i)
        if (x % i == 0) {
            ans = ans / i * (i - 1);
            while (x % i == 0)
                x /= i;
        }
    if (x > 1)
        ans = ans / x * (x - 1);
    return ans;
}
int dfs(int n, int m) {
    if (n == 0)
        return 1;
    if (n < 10)
        return n;
    return qpow(n % 10, dfs(n / 10, phi(m)), m);
}
signed main() {
#ifdef LOCAL
    freopen("input.in", "r", stdin);
    freopen("output.out", "w", stdout);
#endif
    IOS;
    int T;
    cin >> T;
    while (T--) {
        cin >> n >> m;
        cout << dfs(n, m) << "\n";
    }
}

3.Luogu P4139

 相比上面一题,这题可以直接套用欧拉降幂.

#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <map>
#include <vector>
#include <unordered_map>
#include <stack>
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
#define IOS ios::sync_with_stdio(false), cin.tie(0)
#define PII pair<int, int>
#define PDI pair<double, int>
#define int long long
#define ll long long
#define ull unsigned long long
#define debug(a) cout << #a << "=" << a << endl;
#define all(x) (x).begin(), (x).end()
#define sd(x) scanf("%d", &(x))
#define sddd(x, y, z) scanf("%lld%lld%lld", &(x), &(y), &(z))
using namespace std;
const int N = 1e7 + 10;
int phi[N], primes[N], tot;
bool vis[N];
int qmi(int a, int k, int p) {
    int res = 1;
    while (k) {
        if (k & 1)
            res = res * a % p;
        k >>= 1;
        a = a * a % p;
    }
    return res;
}
void euler() {
    phi[1] = 1;
    for (int i = 2; i <= 10000000; ++i) {
        if (vis[i] == 0) {
            primes[++tot] = i;
            phi[i] = i - 1;
        }
        for (int j = 1; j <= tot && i * primes[j] <= 10000000; ++j) {
            vis[i * primes[j]] = 1;
            if (i % primes[j] == 0) {
                phi[i * primes[j]] = phi[i] * primes[j];
                break;
            }else
                phi[i * primes[j]] = phi[i] * (primes[j] - 1);
        }
    }
}
int dfs(int u, int p) {
    if (p == 1)
        return 0;
    return qmi(u, dfs(u, phi[p]) + phi[p], p);
}
void solve() {
    int p;
    cin >> p;
    cout << dfs(2, p) << "\n";
}
signed main() {
#ifdef LOCAL
    freopen("input.in", "r", stdin);
    freopen("output.out", "w", stdout);
#endif
    IOS;
    euler();
    int T;
    cin >> T;
    while (T--)
        solve();
}

野指针*
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
统计数字 k 在 0 到 n 中的出现的次数
李的博客
06-18 1001
题目： 计算数字 k 在 0 到 n 中的出现的次数，k 可能是 0~9 的一个值。 分析： 首先想到使用组合数学的方法，分析n的每一个位上的数字，进行组合运算。分析可知，这种方法依赖于每一位上的数字，程序判断较复杂。 使用从0到n的遍历，找到每一个数的每一位，进行比较。时间复杂度为O(n)O(n)O(n) 代码 class Solution { public: /** * ...
2021-09-23大数据学习日志——Python基础——选择和循环语句（附练习题）
tingbaobaoo的博客
09-23 743
学习目标: 能使用比较运算符与逻辑运算符编写判断条件，例如：年龄满18岁，分数大于70同时少于80 掌握 if 选择语句的基本使用，例如：if语句、if-else、if-elif-else 掌握 while 循环语句的基本使用，例如：语法格式、实现1~100的累加和 能够说出 break 、continue在循环语句中的作用 01_比较和逻辑运算符 学习目标: 能使用比较运算符与逻辑运算符编写判断条件 1.1 比较(即关系，条件)运算符 示例代码： num1 = 15 num2 =
指数循环节
weixin_51652698的博客
10-08 126
指数循环节 定义 若干证明 引入 当n很大时，我们需要降幂进行处理 此时，用到一个新的东西，指数循环节 练习题： Poj2837 套用上面公式（特别是快速幂的部分）即可 此题涉及到快速幂+欧拉函数（可见之前博客）欧拉函数 **注：**只有在B大于等于C的欧拉函数值时，才可以套用此公式。 Calculation（指数循环节，欧拉函数的常用公式） Code： #include <stdio.h> #include <math.h> #include <algorithm&g
MySQL练习题全部
m0_52254303的博客
08-23 2789
sql分析
python大量练习题
m0_74618370的博客
10-14 223
【代码】python写着玩。
【关于 A^x = A^(x % Phi(C) + Phi(C)) (mod C)（x>=phi(c)） 的若干证明】【指数循环节】
忆梦
11-07 682
以下内容全部原创，转载请注明作者 :AekdyCoin以及本文地址 http://hi.baidu.com/aekdycoin/item/e493adc9a7c0870bad092fd9 曾经看过如下一个公式: huicpc0207习题小结：http://blog.csdn.net/longshuai0821/article/details/7827475 以上的公式如果第一次
Python练习—while循环语句&if判断语句
weixin_44795219的博客
10-24 3771
while练习2：有一张厚度为0.1毫米的纸，假设它足够大，重复将其对折，问对折多少次之后，其厚度可以抵达（再对折一次就超过）珠穆朗玛峰的高度？ height=0.1 i=1 while height<=8844000: height=2**i*0.1 i+=1 print(i-1) 2、小明身高1.75，体重80.5kg。请根据BMI公式（体重除以身高的平方）帮小明计算他的BMI指数，并根据BMI指数给出对应的提示信息： 低于18.5：过轻 18.5-25：正常 25-28：过重 28-32：肥胖
第3讲：分支和循环
2301_79811170的博客
10-24 413
如果只想输出开头的部分，可以用%.[m]s指定输出的长度，其中[m]代表一个数字，表示所要输出的长度。上面的示例中，%12f表示输出的浮点数最少要占据12位。前言：下面我们先把上一讲未讲解完的内容先推进完，然后从本讲开始，每次都会在本讲最后面留下一些习题，便于大家更好的掌握每一讲的知识。36⽐较，0
MATLAB课程第4章的练习题：判断、循环语句相关的题目
数学建模学习交流的博客
01-08 2165
Q1. 填空题。
xunhuan.rar_循环节
09-24
给定一个分式A/D，A为分子，D为分母，将分式化解成最简形式，例如1/3 = 0.[3],化成循环小数形式,并用[]将循环节提取出来
C语言循环控制练习题.doc
11-15
C语言循环控制练习题 本文档提供了一组C语言循环控制练习题，涵盖了while循环、do-while循环、循环控制表达式、循环体语句执行次数等知识点。下面对每个题目的知识点进行详细解释。 1. 设有程序段 `int k=10; ...
循环系统疾病 练习题.pdf
03-28
循环系统疾病 练习题.pdf循环系统疾病 练习题.pdf循环系统疾病 练习题.pdf循环系统疾病 练习题.pdf循环系统疾病 练习题.pdf循环系统疾病 练习题.pdf
JAVA入门&循环练习题
03-04
适合刚入门JAVA循环不够熟练的JAVA人拿来练习循环结构，加深对循环的理解
C/C++经典循环章节练习题循环
11-22
C/C++经典循环章节,练习题循环
【C++】——Stack与Queue（含优先队列（详细解读）
m0_74092462的博客
06-09 985
优先队列是一种容器适配器，根据严格的弱排序标准，它的第一个元素总是它所包含的元素中最大的。stack是一种容器适配器，专门用在具有后进先出操作的上下文环境中，其删除只能从容器的一端进行元素的插入与提取操作。 队列是一种容器适配器，专门用于在FIFO上下文(先进先出)中操作，其中从容器一端插入元素，另一端提取元素。
C++——计算不同的非空子串个数
我发布的博客，快来抄答案吧
06-09 138
BCSP-X 11题 C++代码详解
LeetCode 算法：轮转数组c++
小金牛来了
06-09 663
LeetCode 算法：轮转数组c++
9.1 图片的分割处理（c++）
最新发布
weixin_71719718的博客
06-14 603
本文的图片处理分为图片分割、图像的亚像素坐标处理。亚像素处理的原理可以看论文，该论文的详解及c++的代码实现可以看博文。下面的内容很多来自以上博文的内容，我只是对部分代码做了稍微的修改。
python 循环节
09-19
循环节是指字符串中重复出现的部分。对于给定的长度为n的字符串，我们可以通过向左循环移动字符串，找到使得移动后的字符串与原字符串相等的次数。这个次数就是字符串的循环节数。 在Python中，我们可以通过以下...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值