『专题』指数循环节问题

最新推荐文章于 2020-05-30 14:44:27 发布
最新推荐文章于 2020-05-30 14:44:27 发布
阅读量684 收藏
点赞数
分类专栏： 『数据结构与算法』

转载链接:http://blog.csdn.net/acdreamers/article/details/8236942

个人感觉:
最近有点忙,今天才把几道指数循环节题目给刷完.而且过程有点坑,- -在nefuOJ上交题一旦有重定向都给WA了…我搞了一个下午才发现的,我就说我写得没错,略坑…浪费了好多时间.
基本的东西我就不说了,还缺一题是zju的,但OJ崩了没去刷,我只上我做了的代码吧。

Problem 1759 Super A^B mod C

分析:敲门题,难度不高

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e6+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

 ll A;
 char B[maxn];
 ll C;
 ll phi;

 ll ouler(ll c)
 {
     ll res=c;
     for(ll i=2;i*i<=c;i++)
     {
         if(c%i==0)
         {
             res=res/i*(i-1);
             while(c%i==0) c/=i;
         }
     }
     if(c>1) res=res/c*(c-1);

     return res;
 }

 ll getZ()
 {
     ll res=0;
     for(int i=0;B[i]!='\0';i++)
     {
         res=res*10+B[i]-'0';
         res%=phi;
     }

     return res+phi;
 }
 ll qmod(ll a,ll b)
 {
     ll res=1;
     while(b)
     {
        if(b&1)res=(res*a)%C;
        a=(a*a)%C;
        b>>=1;
     }
     return res%C;
 }
int main()
{
#ifndef ONLINE_JUDGE
    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);
#endif
    while(scanf("%I64d%s%I64d",&A,B,&C)!=EOF)
    {
        phi=ouler(C);
        //cout<<phi<<endl;
        ll z=getZ();
        //cout<<"--->"<<z<<endl;

        printf("%I64d\n",qmod(A%C,z));

    }

    return 0;
}

B.Calculation
分析:一丢丢的思维深度,开始理解指循环节!,注意如果你忽略了**（B>=phi§）这个条件,必定让你WA成狗,这个条件不能忽视,**
虽然第一题可以这样但是只是数据有点水,我也不太会举这个样例,

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

int T;
ll n,m;
 ll eulor(ll c)
 {
     ll res=c;
     for(ll i=2;i*i<=c;i++)
     {
         if(c%i==0)
         {
             res=res/i*(i-1);
             while(c%i==0)c/=i;
         }
     }

     if(c>1) res=res/c*(c-1);
     return res;
 }


 ll qmod(ll a,ll b,ll M)
 {
     ll res=1;
     while(b)
     {
         if(b&1)
         {
             res=res*a;
             if(res>=M) res=res%M+M;
         }
         a=(a*a);
         if(a>=M){a=a%M+M;}
         b>>=1;
     }
     return res;
 }
ll f(ll n,ll m)
{
    if(n<10) return n;
    ll p=f(n/10,eulor(m));
    return qmod(n%10,p,m);
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);
#endif
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d%I64d",&n,&m);
        printf("%I64d\n",f(n,m)%m);
    }
    return 0;
}

D:Remainder Calculator

分析:千万别留重定向,肯定会wa了… 你做了第2题这个就很深刻了.

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

ll n,m;
ll a[105];
int T;
bool f;

ll eulor(ll c)
{
    ll res=c;
    for(ll i=2;i*i<=c;i++)
    {
        if(c%i==0)
        {
            res=res/i*(i-1);
            while(c%i==0)c/=i;
        }
    }
    if(c>1) res=res/c*(c-1);
    return res;
}

ll qmod(ll a,ll b,ll M)
{
    ll res=1;
    while(b)
    {
        if(b&1)
        {
            res=res*a;
            if(res>=M) f=true;
            if(f)res%=M;
        }
        a=(a*a);
        if(a>=M) f=true;
        if(f) a%=M;
        b>>=1;
    }
    return res;
}
ll  solve(int step,ll m)
{

    if(step==n)
    {
        ll res=1;

        f=false;
        for(ll i=1;i<=a[step];i++)
        {
            res=(res*i);
            if(res>=m) f=true;
            if(f) res%=m;
        }
        if(f) return res%m+m;
        return res;
    }

    if(a[step]>=m)
    {
        return m;
    }

    ll z=solve(step+1,eulor(m));

    ll d=1;
    for(ll i=1;i<=a[step];i++)
    {
           d*=i;
           if(d>=m)f=true;
           if(f)d%=m;
    }

    ll res=qmod(d,z,m);
    if(f) return res%m+m;
    return res;
}
int main()
{

   scanf("%d",&T);
   while(T--)
   {
       cin>>n>>m;

       for(int i=1;i<=n;i++)cin>>a[i];

       if(m==1)
       {
           printf("0\n");
           continue;
       }

       cout<<(solve(1,m)%m)<<endl;
   }


    return 0;
}

E：What is N

分析:做完之后你会发现,这题就水题了.

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef unsigned long long ull;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }
 int T;
 ull b,p,M;

 ull phi(ull c)
 {
     ull res=c;
     for(ull i=2;i*i<=c;i++)
     {
         if(c%i==0)
         {
             res=res/i*(i-1);
             while(c%i==0) c/=i;
         }
     }
     if(c>1) res=res/c*(c-1);
     return res;
 }
 ull qmod(ull a,ull b,ull M)
 {
     ull res=1;
     while(b)
     {
         if(b&1)res=(res*a)%M;
         a=(a*a)%M;
         b>>=1;
     }
     return res%M;
 }
int main()
{
#ifndef ONLINE_JUDGE
    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);
#endif
    scanf("%d",&T);
    for(int cas=1;cas<=T;cas++)
    {
        scanf("%I64u%I64u%I64u",&b,&p,&M);
        if(b==0&&p==1&&M == 18446744073709551615ull)
        {
           printf("Case #%d: 18446744073709551616\n",cas);
           continue;
        }

        ull ans=0;
        ull mod=phi(p);
        if(b==0)ans++;
        bool flag=false;
        ull z=1;
        ull i;
        for(i=1;i<=M;i++)
        {
              z=z*i;

              if(z>=mod)
              {
                  z%=mod;
                  flag=true;

                  if(z==0)break;
              }

              if(flag)
              {
                  if(qmod(i,z+mod,p)==b) ans++;
              }
              else
              {
                   if(qmod(i,z,p)==b) ans++;
              }
        }
        for(ull k=0;i<=M&&k<p;i++,k++)
        {
            if(qmod(i,mod,p)==b)
            {
               ans+=1+(M-i)/p;
            }
        }


        printf("Case #%d: %I64u\n",cas,ans);


    }
    return 0;
}
Gavinjou大笨象
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【关于 A^x = A^(x % Phi(C) + Phi(C)) (mod C)（x>=phi(c)） 的若干证明】【指数循环节】
08-03 1454
以下内容全部原创，转载请注明作者 :AekdyCoin以及本文地址 http://hi.baidu.com/aekdycoin/item/e493adc9a7c0870bad092fd9 曾经看过如下一个公式: huicpc0207习题小结：http://blog.csdn.net/longshuai0821/article/details/7827475 以上的公式如果第一次
指数循环节
weixin_33972649的博客
12-15 103
$$求证a^b\equiv a^{b\%\varphi(m)+\varphi(m)}(\%m),其中b\geq \varphi(m)$$ 我们模$m$最多只有$m$种结果，所以根据鸽巢原理，在 $$a^0,a^1,...,a^m(\%m)$$ 这$m+1$个数中，一定存在最小的$r$和最小的$s$，满足$a^r\equiv a^{r+s}(\%m)$且$r+s\leq m.$$......(...
关于指数循环节的证明
dolphin198455的博客
04-19 169
关于指数循环节某位大神给了一种简单的证明。 不过前提是知道欧拉定理：A^phi(C)=1(mod C) 就有循环节出现了：A=A^phi(c)*A; 就有 A^x = A^(x % Phi(C) + Phi(C)) (mod C)。 转载于:https://www.cnblogs.com/SXia/p/6747777.html...
Calculation（指数循环节，欧拉函数的常用公式）
C_Dreamy的博客
05-30 1095
Assume that f(0) = 1 and 0^0=1. f(n) = (n%10)^f(n/10) for all n bigger than zero. Please calculate f(n)%m. (2 ≤ n , m ≤ 10^9, x^y means the y th power of x). Input The first line contains a single positive integer T. which is the number of test cases....
HDU2837 Calculation 指数循环节 欧拉函数+快速幂
紫芝的博客
09-19 267
Calculation Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3848    Accepted Submission(s): 925 Problem Description Assume that f(0) = 1 and...
xunhuan.rar_循环节
09-24
在数学领域，循环节是一个与有理数相关的概念...在实际编程中，这些问题可以通过设计算法和数据结构来解决，以识别和表示循环小数的循环节。理解这些概念对于学习计算机科学，特别是数值计算和算法设计，是至关重要的。
xiaoshu.rar_求循环节c_程序循环节
09-22
本程序“xiaoshu.rar”是一个专门用于求解循环小数循环节的C语言实现，对于理解和掌握C语言编程以及循环结构的运用具有很高的实践价值。 首先，我们要理解什么是循环小数。循环小数是小数部分有无限重复数字的小数...
C语言实现小数转分数(包括带循环节小数)
04-29
C语言实现小数转分数，包括带循环节的小数
如何解决Java的循环引用问题
最新发布
10-12
在Java开发中，循环引用是指两个或多个对象相互引用，导致无法被垃圾回收器回收，从而引发内存泄漏的问题。这种问题在大型应用程序中经常会出现，因此解决循环引用问题是非常重要的。 例子及解决方法看附件
单向循环链表解决约瑟夫问题
05-12
用单向循环链表解决约瑟夫问题。使用c++语言，结构体，链表的操作。
CSU2021 Exponial 指数循环节定理
Hagtaril
08-16 374
CSU2021 Exponial 题面: 上图为Exponial(5) 。 给一个 n 和 m 求形如 Exponial(5) 的 Exponial(n) % m 的值 。 解题思路： 指数降幂公式 ： 公式在 K &gt; phi（m）时成立，所以前 5 个 n 要特判。 AC代码: #include&lt;iostream&gt; #incl...
hdu -2837(指数循环节+欧拉函数)
wust_cyl的博客
08-20 605
Problem DescriptionAssume that f(0) = 1 and 0^0=1. f(n) = (n%10)^f(n/10) for all n bigger than zero. Please calculate f(n)%m. (2 ≤ n , m ≤ 10^9, x^y means the y th power of x).InputThe first line cont...
2016多校训练一 PowMod，hdu5728（欧拉函数+指数循环节）
努力
07-20 3414
Declare: k=∑mi=1φ(i∗n) mod 1000000007 n is a square-free number. φ is the Euler's totient function. find: ans=kkkk...k mod p There are infinite number of k   Input Multipl
指数循环节 处理A^B 问题 Super A^B mod C + Calculation
ZJLORD
08-01 967
指数循环节：用于计算   A^B    ; 例子：http://acm.fzu.edu.cn/problem.php?pid=1759 Given A,B,C, You should quickly calculate the result of A^B mod C.             (1 A和B都贼大，所以就需要用一个公式来减少计算； 需要注意的是这个公式有个
（指数循环节降幂）牛客训练赛4.Problem E
weixin_40588429的博客
01-29 106
https://ac.nowcoder.com/acm/contest/330/E 分析题意，求2^n，但n很大。 根据上图公式，因为mod=1e9+7为质数，欧拉函数即为mod-1. #include &lt;iostream&gt; #include &lt;algorithm&gt; #include &lt;cstdlib&gt; #include &lt;cstdio&gt; #i...
Exponial （欧拉定理+指数循环定理+欧拉函数+快速幂）
weixin_30752699的博客
04-18 128
题目链接：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 cons...
bzoj-3884(指数循环节)
wust_cyl的博客
08-21 339
Description根据一些书上的记载，上帝的一次失败的创世经历是这样的：第一天， 上帝创造了一个世界的基本元素，称做“元”。第二天， 上帝创造了一个新的元素，称作“α”。“α”被定义为“元”构成的集合。容易发现，一共有两种不同的“α”。第三天， 上帝又创造了一个新的元素，称作“β”。“β”被定义为“α”构成的集合。容易发现，一共有四种不同的“β”。第四天， 上帝创造了新的元素“γ”，“γ”被定...
指数循环节题集
whai的专栏
01-05 761
指数循环节题集 有个公式 x >= Phi(C), A^x = A ^ (x%Phi(C) + Phi(C)) (mod C) 不会证 1.hdu 3221 Brute-force Algorithm 题目链接 点击打开链接 题意： Function Find(int n,function func)   If n=1   For i=1 to a do func()
解决python循环指数矩阵运算耗时问题
07-28
你可以尝试使用NumPy库来优化Python循环指数矩阵运算的耗时问题。NumPy是一个功能强大的科学计算库，它提供了高效的数组操作和数学函数。 在NumPy中，你可以使用指数函数`np.exp()`来对矩阵进行指数运算。它能够...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值