数论
Lcomyn
这个作者很懒,什么都没留下…
展开
-
【BZOJ 1013】【JSOI 2008】球形空间产生器sphere
高斯消元裸题,设球的一般方程为a^2+b^2+c^2+...+z^2+Aa+Bb+Cc+...+Zz=常数,求解即可code:球心坐标为(A/-2,B/-2,C/-2,...,Z/-2)#include#include#include#includeusing namespace std;double a[15][15],b[15],x[15];int n;int mai原创 2015-02-15 19:54:12 · 696 阅读 · 0 评论 -
莫比乌斯反演学习笔记
莫比乌斯反演,之前做过一些题,一直没有太理解,膜了下faebdc学长的姿势,终于搞懂了一些。 首先我们有两个式子: 1:∑d|nϕ(d)=n\sum_{d|n} \phi(d)=n2:∑d|nμ(d)=e(n)\sum_{d|n} \mu(d)=e(n) 1式证明:对于nn的质因数xx对ϕ(n)\phi(n)贡献了(x−1)∗xt−1(x-1)*x^{t-1} 单独对于xx而言约数可以为x0原创 2015-08-04 20:07:04 · 1640 阅读 · 1 评论 -
【SDOI 2015】【BZOJ 3994】约数个数和
送个结论d(n∗m)=∑d|n∑k|ne(gcd(d,k))d(n*m)=\sum_{d|n}\sum_{k|n} e(gcd(d,k)) 记g(n)=∑ni=1⌊ni⌋g(n)=\sum_{i=1}^{n} \lfloor{n \over i} \rfloor 然后我们开始推 ∑ni=1∑mj=1d(i∗j)=∑ni=1∑mj=1⌊ni⌋⌊mj⌋e(gcd(i,j))=∑ni=1∑mj=1⌊原创 2015-10-22 17:27:21 · 833 阅读 · 0 评论 -
【BZOJ 2671】Calc
题目要求求出∑Na=1∑Nb=a+1[a+b|ab]\sum_{a=1}^{N}\sum_{b=a+1}^{N}[a+b|ab] 我们设d=gcd(a,b),a=nd,b=mdd=gcd(a,b),a=nd,b=md则转化为 [nd+md|nmd2]=[n+m|nmd][nd+md|nmd^2]=[n+m|nmd]又有n,mn,m互质,则显然有n+m|dn+m|d 原式可推为 ∑Nd=1∑⌊原创 2015-09-10 17:28:50 · 1021 阅读 · 0 评论 -
【BZOJ 4253】密码箱加强版
首先无视掉这题题面,您就全当机房的神犇在D蒟蒻我好了。 进入正题,首先这题的弱化版还是蛮好想的。设: n=pk11∗pk22∗...∗pknnn=p_1^{k_1}*p_2^{k_2}*...*p_n^{k_n} 对于x2=1(modn)x^2=1 (mod n),一定有 x2=1(modpk11)x^2=1 (mod p_1^{k_1}) x2=1(modpk22)x^2=1 (mod原创 2015-09-03 17:50:13 · 1334 阅读 · 0 评论 -
【BZOJ 3884】上帝与集合的正确用法
在2与模数pp互质的情况下,我们考虑超级幂,对指数模ϕ(p)\phi(p),那么我们发现这样是可以递归的,那如果pp不与2互质的话,我们可以提出公因子来,就可以转化为互质的情况了。 code:#include<iostream>#include<cstdio>#include<cstring>#include<cmath>using namespace std;int t;int p;原创 2015-09-03 21:06:42 · 721 阅读 · 0 评论 -
【BZOJ 2154】Crash的数字表格
我们定义sum(n,m)=∑ni=1∑mj=1i∗jsum(n,m)=\sum_{i=1}^{n}\sum_{j=1}^{m}i*j ∑ni=1∑mj=1lcm(i,j)=∑ni=1∑mj=1i∗jgcd(i,j)=∑min(n,m)d=1d∑⌊nd⌋i=1∑⌊md⌋j=1di∗djde(gcd(i,j))=∑min(n,m)d=1d∑⌊nd⌋i=1∑⌊md⌋j=1di∗djde(gcd(i,j)原创 2015-08-04 20:29:27 · 1342 阅读 · 1 评论 -
【NOI 2010】能量采集&&超级钢琴
其实这两个题一点关系都没有,同一天做的,写在一起= = 能量采集 ans=∑x=1n∑y=1m2∗gcd(x,y)+1ans=\sum_{x=1}^{n}\sum_{y=1}^{m} 2*gcd(x,y)+1 ∑x=1n∑y=1m[d|gcd(x,y)]=⌊nd⌋∗⌊md⌋\sum_{x=1}^{n}\sum_{y=1}^{m}[d|gcd(x,y)]=\lfloor \frac {n}{d}原创 2015-07-11 17:03:36 · 985 阅读 · 0 评论 -
【BZOJ 2818】Gcd
这题一开始我竟然想要用与能量采集差不多的思路去做= =(no zuo no die,why you try?) 有个显然的转化 ∑nx=1∑ny=1[gcd(x,y)==P]=∑⌊n/P⌋x=1∑⌊n/P⌋y=1[gcd(x,y)==1]=2∗∑⌊n/P⌋i=2ϕ(i)+1\sum_{x=1}^{n}\sum_{y=1}^{n} [gcd(x,y)==P]= \sum_{x=1}^{\lfl原创 2015-07-11 16:48:57 · 845 阅读 · 0 评论 -
【BZOJ 4128】 matrix
题目其实很简单,BSGS可以一眼看出来,但这就需要矩阵求逆,然而我并不会。。。 于是发现了一种BSGS的非求逆方法,借此介绍一下。 aix−y=b(modP)a^{ix-y}=b (mod P) 推出 aix=b∗ay(modP)a^{ix}=b*a^y(mod P) 因此我们只需要对每个b∗ayb*a^yhash一下就好了,然后穷举ii,验证是否存在即可。code:#include<cst原创 2015-07-01 20:45:30 · 1140 阅读 · 0 评论 -
数论学习笔记
素数:#include<iostream>#include<cstdio>#include<cstring>using namespace std;int p[10000000],a[100000001];int main(){ int i,t=0,j,n; scanf("%d",&n); a[1]=1; for (i=2;i<=n;++i)原创 2015-05-28 14:18:08 · 3057 阅读 · 8 评论 -
【HAOI 2007】【BZOJ 1053】反素数ant
虽然这题很水,但蒟蒻还是想了很久。 首先由一个很显然的结论,这题只会用到10个质数。一个感性的证明: 设有一数x=pk11+pk22+...+pkmmx=p_1^{k_1}+p_2^{k_2}+...+p_m^{k_m} 那么我们现在要给xx乘上一个质数,使它的约数数目最多且相对较小。 显然我们我们要乘第1≤i≤m+11\le i \le m+1 个质数,如果乘再靠后的质数,显然不如乘第m+原创 2015-06-03 18:53:32 · 804 阅读 · 0 评论 -
【BZOJ 2190】【SDOI 2008】仪仗队
求出2~n-1内所有数的欧拉函数,求和之后*2再加上3((1,1),(0,1),(1,0)三个点),输出即可code:#include#include#includeusing namespace std;int n;long long sum;long long phi(int n){ long long ans; int m,i; m=floor(sqrt(n));原创 2015-02-16 15:35:19 · 750 阅读 · 0 评论 -
扩展欧几里德与中国剩余定理学习笔记&&模版
蒟蒻这几天比较闲,于是自学了早就学了很久的扩欧与中国剩余定理。。。。。。--------------------------------------------------------------------------------------------------------------扩展欧几里德----------------------------------------------原创 2014-12-02 18:58:39 · 874 阅读 · 0 评论 -
【BZOJ 3560】DZY loves Math V
∑i1|a1∑i2|a2...∑in|anϕ(i1i2..in)\sum_{i_1|a_1}\sum_{i_2|a_2}...\sum_{i_n|a_n} \phi(i_1i_2..i_n) 设共有质因子tottot个,第ii个质因子在第jj个数中指数为kijk_ij则 ans=Πtoti=1∑ji1=0ki1∑ji2=0ki2...∑jin=0kinϕ(prime∑np=1jipi)ans=\原创 2015-11-25 15:36:11 · 1078 阅读 · 0 评论