看了Po神的题解一下子就懂了A了! 不过Po神的代码出锅了-solve中"d-temp"并没有什么用QwQQwQQwQ-应该把模数除以p^temp次方才行. 来自BZOJ讨论板的hack数据 hack data 1 5 3125 7812 正确输出应该是625, 但是很多人输出3125- CODE #include<bits/stdc++.h> using namespace std; typedef long long LL; const LL INF = 1e15; i…

题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2219 N次剩余+CRT... 就是各种奇怪的分类讨论.. #include<cstring> #include<iostream> #include<cstdio> #include<map> #include<cmath> #include<algorithm> #define rep(i,l,r) for (int i=l;i…

中国剩余定理+原根+扩展欧几里得+BSGS 题解:http://blog.csdn.net/regina8023/article/details/44863519 新技能get√: LL Get_yuangen(LL p,LL phi){ ; ;i*i<=phi;i++) ) f[++c]=i,f[++c]=phi/i; ;;g++){ int j; ;j<=c;j++) ) break; ) return g; } ; } 求原根 void Split(int x){ num=; ;i*i&…

#include <iostream> #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; typedef long long int64; ]; int64 A,B,K,pi,pk,ans,ti,q; int64 ksm(int64 x,int64 y){ ) ; ) return x; int64 d=k…

题目描述 在ACM_DIY群中,有一位叫做“傻崽”的同学由于在数论方面造诣很高,被称为数轮之神!对于任何数论问题,他都能瞬间秒杀!一天他在群里面问了一个神题: 对于给定的3个非负整数 A,B,K 求出满足 (1) X^A = B(mod 2*K + 1) (2) X 在范围[0, 2K] 内的X的个数!自然数论之神是可以瞬间秒杀此题的,那么你呢? 输入 第一行有一个正整数T,表示接下来的数据的组数( T <= 1000) 之后对于每组数据,给出了3个整数A,B,K (1 <= A, B <…

原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ2219.html 题目传送门 - BZOJ2219 题意 求同余方程 $x^A\equiv B \pmod{C}$ 的解的个数,其中 $C$ 为一个奇数. $1\leq A,B\leq 10^9,1\leq \lfloor C/2 \rfloor \leq 5\times 10^8$ 题解 UPD(2018-09-10): 详见数论总结. 传送门 - https://www.cnblogs.com/z…

题目描述 终于活成了自己讨厌的样子. 这是她们都还没长大的时候发生的故事.那个时候,栗子米也不需要为了所谓的爱情苦恼. 她们可以在夏日的午后,花大把的时间去研究生活中一些琐碎而有趣的事情,比如数论. 有一天西柚柚问了栗子米一个题,她想知道中有多少不同的数,这些不同的数字里面第k大的是多少. 输入描述: 第一行一个整数T(T≤ 10 5 ),表示数据组数.每组数据第一行两个整数,表示n,k(1≤ n≤ 10 18 ),保证k不会超过不同的数字个数. 输出描述: 对于每组数据输出,输出两个整数,表示…

今年的重庆省选? 具体就是,对于每次修改,A[p,q]这个位置,  设d=gcd(p,q) ,则 gcd为d的每一个格子都会被修改,且他们之间有个不变的联系 A[p,q]/p/q==A[k,t]/k/t   所以只要记录对于gcd为d的所有格子,只要保存A[d][d]的值就可以了. 那么求前k行k列的值ans,则所有gcd(p,q)==d的A[p,q]对答案的贡献就是    { 设k'=k/d;  (下取整)  f[k']*A[p,q]/(p/d)/(q/d) } 首先有个基本结论(当n>1时)…

这题的DP很难想,定义dp[i][j][a][b]表示用了i个男生,j个女生,任一连续的后缀区间内,男生比女生最多多a人,女生比男生最多多b人. 转移就是显然了. # include <cstdio> # include <cstring> # include <cstdlib> # include <iostream> # include <vector> # include <queue> # include <stack&…

首先问题的意思就是在找出n以内的所有x^2%n=1的数,那么我们可以得到(x+1)(x-1)=y*n,那么我们知道n|(x+1)(x-1),我们设n=a*b,那么我们对于任意的a,我们满足n%a==0,我们可以求出b,我们可以使x+1|a,x-1|b,然后我们可以构造所有满足被b整除的数,然后判断是否能被a整除, 然后再枚举x+1|b,x-1|a的情况,假设一组合法解不能拆开后被a,b分别整除,那么对于另外的a,b我们肯定可以再次枚举出这个解,然后对于相同的解用set去下重就可以了. 反思:手残…