利用降幂公式。。呃,还是自己去搜题解吧。知道降幂公式后,就不难了。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define LL unsigned long long
using namespace std; bool mod[100005]; LL PHI(LL P){
LL ret=1;
for(LL i=2;i*i<=P;i++)
if(P%i==0){
ret*=i-1;
P/=i;
while(P%i==0){
P/=i;
ret*=i;
}
}
if(P>1)
ret*=P-1;
return ret;
} LL quick(LL a,LL k,LL m){
LL ret=1; LL t=a%m;
while(k){
if(k&1) ret=(ret*t)%m;
k>>=1;
t=(t*t)%m;
}
return ret;
} int main(){
int T;
LL b,P,M;
scanf("%d",&T);
for(int kase=1;kase<=T;kase++){
scanf("%I64u%I64u%I64u",&b,&P,&M);
printf("Case #%d: ",kase);
if(P==1){
if(M==18446744073709551615ULL)
printf("18446744073709551616\n");
else
printf("%I64u\n",M+1);
continue;
}
LL phi=PHI(P);
LL i,ans;
LL fac=1; ans=0;
for(i=0;i<=M&&fac<=phi;i++){
if(quick(i,fac,P)==b){
ans++;
}
fac*=(i+1);
}
fac%=phi;
for(;i<=M&&fac;i++){
if(quick(i,fac,P)==b)
ans++;
fac=(fac*(i+1))%phi;
}
if(i<=M){
memset(mod,false,sizeof(mod));
LL cn=0;
for(LL k=0;k<P;k++){
if(quick(i+k,phi,P)==b){
mod[k]=true;
cn++;
}
}
ans+=((M-i+1)/P)*cn;
LL e=(M-i+1)%P;
for(LL k=0;k<e;k++)
if(mod[k])
ans++;
}
printf("%I64u\n",ans);
}
return 0;
}

  

HDU 4335 Contest 4的更多相关文章

  1. 数论 + 公式 - HDU 4335 What is N?

    What is N? Problem's Link: http://acm.hdu.edu.cn/showproblem.php?pid=4335 Mean: 给你三个数b.P.M,让你求有多少个n满 ...

  2. HDU 5045 Contest(状压DP)

    Problem Description In the ACM International Collegiate Programming Contest, each team consist of th ...

  3. hdu - 5045 - Contest(国家压缩dp)

    意甲冠军:N个人M通过主打歌有自己的期望,每个问题发送人玩.它不能超过随机播放的次数1,追求最大业绩预期 (1 ≤ N ≤ 10,1 ≤ M ≤ 1000). 主题链接:pid=5045" ...

  4. [ACM] hdu 5045 Contest (减少国家Dp)

    Contest Problem Description In the ACM International Collegiate Programming Contest, each team consi ...

  5. HDU–5988-Coding Contest(最小费用最大流变形)

    Coding Contest Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)To ...

  6. hdu 4335 What is N?

    此题用到的公式:a^b%c=a^(b%phi(c)+phi(c))%c (b>=phi(c)). 1.当n!<phi(p)时,直接暴力掉: 2.当n!>=phi(p) &&a ...

  7. hdu 5045 Contest(状态压缩DP)

    题解:我们使用一个二位数组dp[i][j]记录进行到第i个任务时,人组合为j时的最大和(这里的j我们用二进制的每位相应一个人). 详细见代码: #include <iostream> #i ...

  8. HDU 5045 Contest

    pid=5045">主题链接~~> 做题感悟:比赛时这题后来才写的,有点小尴尬.两个人商议着写写了非常久才写出来,I want to Powerful ,I believe me ...

  9. HDU 4339 Contest 4

    树状数组,主要是抓住要求连续1的个数.这样,初始时,相同的加1,不同的加0. 查询时,用二分搜索右边界.就是比较当前mid-l+1的值与他们之间1的个数(这可以通过树状数组求区间和得出),记录右边界即 ...

随机推荐

  1. 计算机-award BIOS全教程

  2. js+ canvas 实现人物走动

    在网上看了一篇管道工玛利亚走动的图片,感觉人物走动的太生涩了,就写了一下代码改动一下: js 代码: //定义数组图片集合 var marios = new Array("image/QQ截 ...

  3. python的range()函数使用方法

    python的range()函数使用非常方便.它能返回一系列连续添加的整数,它的工作方式类似于分片.能够生成一个列表对象. range函数大多数时常出如今for循环中.在for循环中可做为索引使用.事 ...

  4. [JavaEE]Hibernate 所有缓存机制详解

    Hibernate提供的一级缓存 hibernate是一个线程对应一个session,一个线程可以看成一个用户.也就是说session级缓存(一级缓存)只能给一个线程用,别的线程用不了,一级缓存就是和 ...

  5. oracle数据泵备份与还原

    完整的常用的一套oracle备份以及还原方案 --在新库中新建数据目录,我没有特别说明在哪执行的语句都可在plsql中执行 CREATE OR REPLACE DIRECTORY dump_dir A ...

  6. 浅谈javascript的面向对象思想

    面向对象的三大基本特性 封装(把相关的信息(无论数据或方法)存储在对象中的能力) 继承(由另一个类(或多个类)得来类的属性和方法的能力) 多态(一个对象在不同情况下的多种形态) 定义类或对象 第一种: ...

  7. 关于requestAnimationFrame()和cancelAnimationFrame()与定时器之间的比较

    在渲染页面动画的时候,其实也没有必要用定时器(setInterval),其实requestAnimationFrame()和cancelAnimationFrame()也能达到相应的效果,他是HTML ...

  8. <错误>

    1. package com.multak.cookaraclient.adapter; import android.content.Context; import android.support. ...

  9. CentOS 7.1 下载,安装,配置

    CentOS(全称 Community Enterprise Operating System)是 RHEL 源代码再编译的产物. 2014年初,CentOS 宣布加入 Red Hat. CentOS ...

  10. Mac-O文件加载的全过程(一)

    在Mac的开发中, 有没有想过当我们点击可执行文件之后,Mac究竟做了什么事情才让我们的程序运行起来? 对于应用层开发人员或者普通的用户而言, 其实无需知道的这么详细:但是对于内核开发人员而言, 如果 ...