数论 + 公式 - 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满足下式。

analyse:

看到数据被吓到了,没半点思路,后来看了解题报告,方法竟然是暴力!

当然暴力是有条件的。

有这样一个公式:

A^x = A^(x % Phi(C) + Phi(C)) (mod C) (x>=Phi(C))

这个公式的具体证明原来在aekdycoin的百度空间有,但是随着百度空间被转移(百度作死,流失了好多优质的文章==),这篇文章的完整版也流失了。

我们就当这个公式是定理吧!

当n!<Phi(C)时,此时我们暴力解决就可。

当n!大于phi(P)的时候,就需要用上面的降幂公式了。

方法还是暴力,n!%phi(p)会出现0,这是必然的,至少n>=phi(p)为0,

那么(n+1)!%phi(p)也为0,这便出现了重复,转变为n^(phi(p))%p==b的问题了。

固定了指数,根据鸽巢原理,余数是循环的,那么只要找出p个的结果,之后通过循环节求解便可以了。

Trick:当P为1的时候,b为0,这时候答案是m+1,不过m可能为2^64-1,如果加1的话就会溢出,巨坑。

Time complexity: O(N)

Source code:

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-08-25-23.41
* Time: 0MS
* Memory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
typedef __int64(LL);
typedef unsigned __int64(ULL);
const double eps(1e-8);

LL get_eular(LL m)
{
     LL ret=1;
     for(LL i=2; i*i<=m; i++)
           if(m%i==0)
           {
                 ret*=i-1;
                 m/=i;
                 while(m%i==0)
                 {
                       m/=i;
                       ret*=i;
                 }
           }
     if(m>1) ret*=m-1;
     return ret;
}

long long Quickpow(long long a,long long b,long long m)
{
     long long ans=1;
     while(b)
     {
           if(b&1) { ans=(ans*a)%m,b--; }
           b/=2,a=a*a%m;
     }
     return ans;
}

LL b,p,m,ring[100010];
int main()
{
     int t,Cas=0;
     scanf("%d",&t);
     while(t--)
     {
           scanf("%I64u %I64u %I64u",&b,&p,&m);
           if(p==1)
           {
                 if(m==18446744073709551615ULL)
                       printf("18446744073709551616\n");
                 else
                       printf("%I64u\n",m+1);
                 continue;
           }
           LL i=0,phi=get_eular(p),fac=1,ans=0;
           for(i=0; i<=m&&fac<=phi; i++)
           {
                 if(Quickpow(i,fac,p)==b)
                       ans++;
                 fac*=i+1;
           }
           fac=fac%phi;
           for(; i<=m&&fac; i++)
           {
                 if(Quickpow(i,fac+phi,p)==b)
                       ans++;
                 fac=(fac*(i+1))%phi;
           }
           if(i<=m)
           {
                 LL cnt=0;
                 for(int j=0; j<p; j++)
                 {
                       ring[j]=Quickpow(i+j,phi,p);
                       if(ring[j]==b)
                             cnt++;
                 }
                 LL idx=(m-i+1)/p;
                 ans+=cnt*idx;
                 LL remain=(m-i+1)%p;
                 for(int j=0; j<remain; j++)
                       if(ring[j]==b)
                             ans++;
           }
           printf("Case #%d: %I64u\n",++Cas,ans);
     }
     return 0;
}/

时间: 2024-10-11 22:14:52

数论 + 公式 - HDU 4335 What is N?的相关文章

【中国剩余定理】【容斥原理】【快速乘法】【数论】HDU 5768 Lucky7

题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5768 题目大意: T组数据,求L~R中满足:1.是7的倍数,2.对n个素数有 %pi!=ai  的数的个数. 题目思路: [中国剩余定理][容斥原理][快速乘法][数论] 因为都是素数所以两两互素,满足中国剩余定理的条件. 把7加到素数中,a=0,这样就变成解n+1个同余方程的通解(最小解).之后算L~R中有多少解. 但是由于中国剩余定理的条件是同时成立的,而题目是或的关系,所以要用容斥原理叠加删

数论之 HDU 2866

http://acm.hdu.edu.cn/showproblem.php?pid=2866 水题 题意:找 [2,L]内有多少个 p 满足 该式子,其中 n m 可以为任意整数,p为素数. 别人的gcd思路: n^b + p*n^(b-1) = m^b   ==>   n^(b-1)*[n+p]=m^b 因为n里面要么有p因子,要么没有,所以gcd(n^(b-1),n+p)=1或(含有p因子的数) 当gcd(n^(b-1),n+p)== (含有p因子的数)的时候,显然无解,因为假设有解,那么n

推公式 HDU 2552

T 给你2个值 求另外一个 需要推一下 tan(a+b)=(tan(a)+tan(b))/(1-tan(a)*tan(b)); 等式左右取tan tan(atan(a))=a 1/s=tan(...)=(1/u+1/v)/(1-1/(uv)); 最后推出 那个式子等于1 #include<stdio.h> #include<string.h> #include<algorithm> #include<math.h> using namespace std;

HDU 4335 Party All the Time(三分|二分)

 #include<cstdio> #include<cstring> #include<cmath> #include<cstdlib> #include<iostream> #include<algorithm> #include<vector> #include<map> #include<queue> #include<stack> #include<string>

Humble Numbers数论(HDU 1055)

Humble Numbers Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 18028    Accepted Submission(s): 7845 Problem Description A number whose only prime factors are 2,3,5 or 7 is called a humble numb

数论公式总结

目录 1.中间式子&常用技巧 2.gcd相关 3.d相关(d是约数个数函数) 4.mu相关(mu是莫比乌斯函数) 5.sigma相关(sigma是约数和函数) 6.因子相关 更新中... 1.中间式子&常用技巧 \[[n==1]=\sum_{d|n}\mu (d)\] 这个式子用来替换条件式,从而降低复杂度 \[\sum_{i=1}^n\sum_{j=1}^{[\frac ni]}f(i)=\sum_{i=1}^n\sum_{j=1}^{[\frac ni]}f(j)\] 被加数可以任意以

HDU6706 CCPC 2019网络赛 huntian oy 推式子+杜教筛

CCPC 2019 网络赛 HDU 6706 huntian oy 标签 奇奇怪怪的数论结论 杜教筛 前言 我的csdn和博客园是同步的,欢迎来访danzh-博客园~ 简明题意 给定n,a,b,求: \[\sum_{i=1}^n\sum_{j=1}^igcd(i^a-j^a,i^b-j^b)[gcd(i,j)=1]\%(10^9+7)\] 思路 首先有一个结论: \[gcd(i^a-j^a,i^b-j^b)=i^{gcd(a,b)}-j^{gcd(a,b)}\] 上面的结论对于i,j互质是成立的

waiting

rand()函数学习 各OJ网址 数论相关 hdu分类相关 特别的DP教程

HDU 4861 Couple doubi(数论)

HDU 4861 Couple doubi 题目链接 题意:给定k,p,有k个球,每个球的值为1^i+2^i+...+(p-1)^i (mod p) (1 <= i <= k),现在两人轮流取球,最后球的值总和大的人赢,问先手是否能赢 思路:先手不可能输,非赢即平,那么只要考虑每种球的值, 利用费马小定理或欧拉定理,很容易得到该函数的循环节为p - 1, 那么i如果为p - 1的倍数,即为循环节的位置,那么每个值都为1,总和为p - 1 如果i不在循环节的位置,任取一个原根g,根据原根的性质,