码迷,mamicode.com
首页 > 其他好文 > 详细

Exponial (欧拉定理+指数循环定理+欧拉函数+快速幂)

时间:2018-04-18 21:54:29      阅读:163      评论:0      收藏:0      [点我收藏+]

标签:hat   int   mon   and   struct   cal   frame   链接   代码   

题目链接:http://acm.csu.edu.cn/csuoj/problemset/problem?pid=2021

Description

Everybody loves big numbers (if you do not, you might want to stop reading at this point). There are many ways of constructing really big numbers known to humankind, for instance:

  • Exponentiation: 422016=4242...422016 times422016=42⋅42⋅...⋅42?2016 times.
  • Factorials: 2016!=2016?⋅?2015?⋅?...?⋅?2?⋅?1.

技术分享图片

In this problem we look at their lesser-known love-child the exponial, which is an operation defined for all positive integers n? as
exponial(n)=n(n?−?1)(n?−?2)?21
For example, exponial(1)=1 and exponial(5)=54321?≈?6.206?⋅?10183230 which is already pretty big. Note that exponentiation is right-associative: abc?=?a(bc).

Since the exponials are really big, they can be a bit unwieldy to work with. Therefore we would like you to write a program which computes exponial(n) mod m (the remainder of exponial(n) when dividing by m).

Input

There will be several test cases. For the each case, the input consists of two integers n (1?≤?n?≤?109) and m (1?≤?m?≤?109).

Output

Output a single integer, the value of exponial(n) mod m.

Sample Input

2 42
5 123456789
94 265

Sample Output

2
16317634
39

思路:本题是一道经典的指数循环定理简记e(n)=exponial(n)e(n)=exponial(n),利用欧拉定理进行降幂即可,不过要注意会爆int。指数循环公式为指数循环公式为A^B = A^(B %  φ(C) +  φ(C)) % C,其中 φ(C)为1~C-1中与C互质的数的个数。
技术分享图片

 


代码如下:
 1 #include <cstdio>
 2 #include <cstring>
 3 
 4 typedef long long ll;
 5 int n, m;
 6 
 7 ll euler(int n) {
 8     ll ans = n;
 9     for(int i = 2; i * i <= n; i++) {
10         if(n % i == 0) {
11             ans = ans / i * (i - 1);
12             while(n % i == 0) n /= i;
13         }
14     }
15     if(n > 1) ans = ans / n * (n - 1);
16     return ans;
17 }
18 
19 ll ModPow(int x, int p, ll mod) {
20     ll rec = 1;
21     while(p > 0) {
22         if(p & 1) rec = (ll)rec * x % mod;
23         x = (ll)x * x % mod;
24         p >>= 1;
25     }
26     return rec;
27 }
28 
29 ll slove(int n, ll m) {
30     if(m == 1) return 0;
31     if(n == 1) return 1 % m;
32     if(n == 2) return 2 % m;
33     if(n == 3) return 9 % m;
34     if(n == 4) return (1 << 18) % m;
35     return (ll)ModPow(n, euler(m), m) * ModPow(n, slove(n - 1, euler(m)), m) % m;
36 }
37 
38 int main() {
39     while(~scanf("%d%d", &n, &m)) {
40         printf("%lld\n",slove(n, m));
41     }
42     return 0;
43 }

 

有不懂的请私聊我QQ(右侧公告里有QQ号)或在下方回复

 

Exponial (欧拉定理+指数循环定理+欧拉函数+快速幂)

标签:hat   int   mon   and   struct   cal   frame   链接   代码   

原文地址:https://www.cnblogs.com/Dillonh/p/8877711.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!