数论
文章平均质量分 66
thy_asdf
我也不知道说什么...
展开
-
bzoj1477: 青蛙的约会
传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1477思路:扩展欧几里德求解同余方程。设x步后碰面,初始位置为x0,y0。那么就有x0+mx=y0+nx(mod L)(m-n)x=y0-x0(mod L) //程序里要先把m-n变成正数,L已经是正数,不用变,负数模会有问题,这两步的等号是同余。于是就可以把这个同原创 2015-10-11 10:32:56 · 730 阅读 · 0 评论 -
bzoj1407: [Noi2002]Savage
传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1407思路:因为M山洞数量满足要求,就是对于任何两个野人组成的同余方程最小解大于寿命短的野人的寿命,或者没有解。代码:#include#include#include#define abs(a) (a<0?-a:a)const int maxn=100000原创 2015-10-11 10:44:38 · 770 阅读 · 0 评论 -
一道组合数取模题
题目大意:求长度为n且每项均在[1,n]的不上升数列与不下降数列的个数和。思路:总数就是不下降数列的个数*2-n(常数列的个数)然后考虑不下降数列的个数为了方便,把第0项设为0,把第n+1项设为n。差分,然后不下降数列就是差分数组a[i]每一项大于等于0,且Σa[i]=n。每项+1,就相当于在2n-1个空位(本来是2n+1,首尾不能放)放n个板子。于是答案就是C(2n-1,原创 2015-10-11 11:10:50 · 367 阅读 · 0 评论 -
bzoj2075[SDOI2012]Longge的问题
传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=2705思路:首先要求的是Σgcd(i,n)(1于是可以化成ΣΣd*[gcd(i,n/d)==1] (d|n,1把d提出来Σd*Σ[gcd(i,n/d)==1] (d|n,1右边的不就是phi(欧拉函数)吗于是有Σd*phi(n/d) (d|n)因为原创 2015-10-12 09:17:20 · 555 阅读 · 0 评论 -
spoj GCDEX - GCD Extreme
传送门:http://www.spoj.com/problems/GCDEX/思路:令g(n)=ΣΣgcd(i,j)(i很显然ans=g(n)-(n+1)*n/2那么g(n)=Σf(i)(if(n)在上一题已证为积性函数,线筛求出,就可以在O(n)时间打出表来。为了方便先设mindiv[i]为i的最小质因子,T[i]为i的最小质因子次数,d[i]=mindiv[i]^t[i]原创 2015-10-12 11:24:08 · 637 阅读 · 0 评论 -
bzoj2045: 双亲数&bzoj1101: [POI2007]Zap
传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1101思路:设a求的就是f(a,b,d)=ΣΣ[gcd(i,j)==d](1又有公式Σmiu(d) =[n==1] (d|n) //miu是莫比乌斯函数把n用gcd(i,j)代入得ΣΣΣmiu(d)(d|gcd(i,j),1=Σmiu(d)*(a/d)*(b/d原创 2015-10-12 17:04:59 · 515 阅读 · 0 评论 -
bzoj2820: YY的GCD
传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=2820思路:#include#include#includeconst int N=10000010;using namespace std;typedef long long ll;int mu[N+10],g[N+10],pri[原创 2015-10-13 11:25:12 · 438 阅读 · 0 评论 -
bzoj3667: Rabin-Miller算法
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3667思路:首先我们说说Miller_Rabin算法我们发现了费马小定理那它倒过来对不对呢如果a^(p-1)=1(mod p),那么p一定是素数吗?很不幸,是错的虽然出错概率很低,但是可以被卡于是我们就给它打补丁我们又找到了一个二次探测的方法如果p是质数,那原创 2016-05-08 22:39:04 · 2540 阅读 · 1 评论