
数论
文章平均质量分 74
自ran而然
这个作者很懒,什么都没留下…
展开
-
ACdream 1084 寒假安排 (算术基本定理)
题目传送门:点击打开链接这道题做之前,我觉得先懂得怎么求n!中某质因子的个数。。。如果不知道请看大神博客点击打开链接。在这里我把他的模版复制下来int count(int n,int k){ int num=0; while(n) { num+=n/k; n/=k; } return num;}原创 2016-11-03 20:35:31 · 345 阅读 · 0 评论 -
欧拉函数模版
ll phi(ll n) { ll i,rea=n; for(i=2;i*i<=n;++i) { if(n%i==0) { rea=rea-rea/i; while(n%i==0) n/=i; }原创 2017-08-19 19:58:26 · 394 阅读 · 0 评论 -
中国剩余定理模版
typedef long long LL;const int N = 1005;LL a[N], m[N];LL gcd(LL a,LL b){ return b? gcd(b, a % b) : a;}void exgcd(LL a, LL b, LL &x, LL &y){ if(b == 0) { x = 1; y = 0;原创 2017-08-16 19:41:13 · 360 阅读 · 0 评论 -
HDU 4569 Special equations 枚举
点击打开链接本题的意思:已知f(x) = a nx n+...+ a 1x +a 0 ,整数m,m一个素数的平方。求出x使f(x) mod (m)=0;给出a nto a 0,(0 = 3, otherwise abs(a i) 解题思路:直接枚举(1-1e8)的话会超时因为m=pri*pri,且 m|f(x),那么pri| f(x),所以我们从0-pre-1枚举x0是否符合原创 2017-06-02 12:35:46 · 534 阅读 · 0 评论 -
HDU 4565 So Easy!
点击打开链接本题题意:给出整数a,b,n,m求出Sn。l0解题思路:因为a,b,n很大所以如果直接求(a+sqrt(b))^n会爆double范围 所以只能对公式化简。由题中b的范围为(a-1) 2代码:#include #include #include using namespace std;typedef long long ll;const dou原创 2017-06-02 12:32:04 · 486 阅读 · 0 评论 -
HDU - 4810 Wall Painting(二进制+组合数学)
题意:输出第x天 从n中选x个数进行异或的 和思路:看到异或 一般与二进制有关系 剩下就是记录二进制的个数然后用组合数学就可以搞定了代码:#include #include #define rep(i,o,u) for(int i=o;i<=u;++i)using namespace std;typedef long long ll;const int mod=1e6+3;原创 2017-04-20 18:40:31 · 402 阅读 · 0 评论 -
hdu 4407,4059,3501 容斥原理
hdu 3501 传送门:点击打开链接题意:求小与整数n不互质的数的sum。思路:这题有2种解法 先说第一种;用容斥原理。与n不互质 ,即gcd(x,n)!=1,即 符合这样的数 都是n的素因子的倍数,所以我们先求出n的素因子,然后用容斥来求就可以了。这里我们对于 p(素数)的倍数的和 怎么求进行一下说明。我们知道 小于n并且是p的倍数 有 p,2*p,3*p,4*p......原创 2016-11-23 17:22:07 · 564 阅读 · 0 评论 -
hdu 1695 hdu 2841 容斥原理
1695 题目传送门: 点击打开链接题目大意:给你a,b,c,d,k, 找出 gcd(x,y)=k的 组数,其中 a思路:由最大公约数 性质 ,我们可以将gcd(x,y)=k,化简为 gcd(x/k,y/k)=1。那么 也就是说 我们要找出x/k 与y/k互质的个数。首先1代码:#include #include #define rep(i,a,b) for原创 2016-11-22 15:15:42 · 374 阅读 · 0 评论 -
LightOJ 1236 Pairs Forming LCM (算术基本定理)
题目传送门:点击打开链接这题要求a,b两个数的最小公倍数等于n的对数,其中a,b根据算术基本定理:n=p1^x1*p2^x2*p3^x3*p4^x4...a=p1^a1*p2^a2*p3^a3*p4^a4.....b=p1^b1*p2^b2*p3^b3*p4^b4.....其中pi为素数如果a,b的公倍数为n,那么n=p1^max(a1,b1)*p2^max(a2,b2)*原创 2016-11-02 18:10:59 · 396 阅读 · 0 评论 -
POJ 3904 HDU 5072 容斥原理
poj 3904 点击打开链接 hdu 5072 点击打开链接poj 3904 题意 :在给出的n个数中,找出4个数 符合gcd(a,b,c,d)=1的 组数。hdu 5072 题意:再给出n 个数中,找出 3个数 符合(a,b)=(b,c)=(a,c)=1或 (a,b)!=1,(b,c)!=1,(a,c)!=1 的组数。。。总的来说 :这两题的思路 都是把这n个数原创 2016-11-29 09:35:40 · 414 阅读 · 0 评论 -
hdu 5768 Lucky7 容斥原理 中国剩余定理
题目传送门:点击打开链接题意:给你区间[a,b],在这区间里 找出 既是7的倍数 x,并且 x%pi!=ai的个数 有多少?思路:看到这题 跟模有关系 大概跟 中国剩余定理有关,可是题目说的 只要 符合 x%pi=ai一项 就不能算个数,所以只是中国剩余定理 不可行的,那么我们只有一个一个 求出个数,然后再用7的倍数的个数-求出的个数,可是 一个一个求这会有重复的,所以我们 就要用到容原创 2016-11-28 19:17:56 · 577 阅读 · 0 评论 -
hdu 5317 (普通素数筛的应用)+(前缀和预处理)
RGCDQ Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 2917 Accepted Submission(s): 1126Problem DescriptionMr. Hdu is interest原创 2016-10-28 16:49:37 · 778 阅读 · 1 评论 -
hdu 2588(欧拉函数)
GCDTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1834 Accepted Submission(s): 917Problem DescriptionThe greatest common div原创 2016-10-26 16:33:46 · 381 阅读 · 0 评论 -
CodeForces - 615D Multipliers 费马小定理+组合数学
题目传送门:点击打开链接题意:求出整数n的所有因子的乘积对1e9+7的模。从题目中的给出的例子1·2·3·6 = 36 可以看出前面的式子可以写成 (1*6)*(2*3)=6*6=36,而n=6,所以我们可以想到答案可以是 (n)^m这种类型。m=n的因子个数除以2.那么n的因子个数怎么求呢?根据题目给的数据,我们可以知道n的素因子种类可相应的个数,那我们可以假设 n的一个素因子为p原创 2016-11-04 21:43:21 · 599 阅读 · 0 评论 -
伯努利求幂和模版
∑ni=1ik∑i=1nik\sum_{i=1}^n{i^k}=(n+B)k+1−Bk+1k+1(n+B)k+1−Bk+1k+1\frac{{(n+B)}^{k+1}-B^{k+1}}{k+1}原创 2018-07-20 17:15:49 · 375 阅读 · 0 评论