【莫比乌斯反演】

Visible Lattice Points 莫比乌斯反演

A-莫比乌斯反演TimeLimit:1368MS    MemoryLimit:1572864KB    64bitIOFormat:%lld&%lluSubmitStatusDescriptionConsideraN
zp___waj·2016-04-19 14:00

【欧拉函数，莫比乌斯反演】

【普通版】首先是一道单纯的容斥题目（hdu4135）：求[a,b]内与n互质的数的个数，n #definelllonglong usingnamespacestd; intmain(){ intt,cnt=0; cin>>t; while(t--){ lla,b,n; scanf("%I64d%I64d%I64d",&a,&b,&n); vectorv; for(inti=2;i*i1) v.pu
cacyth·2016-04-19 00:00

PE 439 【莫比乌斯反演】【杜教筛】

垃圾题，坑了我一天QAQ题目大意：求∑Ni=1∑Nj=1σ(ij),N=1011∑Ni=1∑Nj=1σ(ij)=∑Ni=1∑Nj=1∑x|i∑y|jxjy[gcd(x,y==1)]=∑Ni=1∑Nj=1∑x|i∑y|jxjy∑d|x,yμ(d)=∑Nd=1μ(d)∑⌊nd⌋i=1id⌊nid⌋∑⌊nd⌋j=1j⌊nid⌋=∑Nd=1dμ(d)(∑⌊nd⌋i=1σ(i))2这貌似跟SDOI2015的
beginendzrq·2016-04-18 14:00

hdu1695莫比乌斯反演

#include #include #include #include #include #include #include #include #include typedeflonglongLL; constintmaxn=100000; intprime[maxn+8],primeSize; boolis[maxn+8]; intmu[maxn+8]; voidmobius(){ mem
u013491262·2016-04-15 10:00

老逗的gcd 莫比乌斯反演

这个题一看就是莫比乌斯反演，怎么处理无平方因子数呢？
liuguangzhe1999·2016-04-14 09:00

[置顶] CQOI2016省选游记

莫比乌斯反演？连通性状压？二项式定理？统统SMG……然后没有等成绩直接回来，记得那天在讲置换群，现在那几道题还没做完……哦对了我好像为全场创造了垫底……   这一次，CQ第一次NOI规则、细致
liuyunhui246·2016-04-12 15:00

莫比乌斯反演之入门

前言很久以前就学了莫比乌斯反演，然而一直都木有来写一个总结，省选完后今日来补坑…姿势其实莫比乌斯反演就是一个公式…F(n)=∑d|nf(d)⇒f(d)=∑d|nμ(d)∗F(nd)证明如下：∑d|nμ(
cqbztsy·2016-04-11 20:00

bzoj3601 一个人的数论 高斯消元&莫比乌斯反演

AC代码如下：#include #include #include #definelllonglong #definemod1000000007 #defineinv(x)ksm(x,1000000005) usingnamespacestd; intm,n,a[105][105],p[1005],q[1005]; intksm(intx,inty){ intt=1;if(y>=1,x=(ll)
lych_cys·2016-04-11 15:00

[BZOJ2005][NOI2010]能量采集（莫比乌斯反演）

题目描述传送门题解首先证明对于某个点（x,y）,k=gcd(x,y)-1：设gcd(x,y)=t，令x=at，y=bt，那么在这条直线上的整数点可以表示为(a,b)(2a,2b)(3a,3b)……(x,y)，由于不算x,y，则答案为gcd(x,y)-1那么总损耗2k+1=2×gcd(x,y)-1。我们最终要求的式子为：∑i=1n∑j=1m(gcd(i,j)∗2−1)=2∗∑i=1n∑j=1mgcd
Clove_unique·2016-04-07 20:43

bzoj3529 数表 莫比乌斯反演&树状数组

    然后按照f[i]和a分别排序之后离线插入用树状数组维护前缀和即可。AC代码如下：#include #include #include #defineN100005 usingnamespacestd; intcnt,pt,mx,p[N],mu[N],c[N],prm[N],ans[N]; structnode{intx,y,k,id;}a[N],f[N]; intread(){ intx
lych_cys·2016-04-04 20:00

ACM

1、数论素数测试（筛法），扩展欧几里得算法，同余模运算，高斯消元，中国剩余定理，莫比乌斯反演2、图论最短路，最小生成树，拓扑排序，二分图，最大团，最大流，强连通分量，最近公共祖先，次小生成树，欧拉回路，
mymilkbottles·2016-04-04 12:00

BZOJ2671: Calc

重新推了一下发现就是莫比乌斯反演暴力可过。。。
liutian429073576·2016-04-03 18:00

hdu1695 GCD（莫比乌斯反演）

题意：求（1，b）区间和（1，d）区间里面gcd(x,y)=k的数的对数（1 #include #include usingnamespacestd; constintN=100000+10; intu[N],prime[N]; boolvis[N]; voidinit() { memset(vis,0,sizeof(vis)); u[1]=1; intcnt=0; for(inti=2;i>t
strangedbly·2016-04-03 13:00

HDU 5656 (莫比乌斯反演)

CALovesGCDTimeLimit:6000/3000MS(Java/Others)    MemoryLimit:262144/262144K(Java/Others)TotalSubmission(s):61    AcceptedSubmission(s):14ProblemDescriptionCAisafinecomradewholovesthepartyandpeople;inev
morejarphone·2016-04-02 21:00

【bzoj3994】[SDOI2015]约数个数和 线性筛法+莫比乌斯反演+数论分块

这题有个奇怪的性质，发现这个性质才可做。为什么呢？考虑d是p的倍数，i是d的倍数，d可以是p的1~n/p倍，当d=kp时，i可以选n/(pk)个数。f数组怎么处理呢？O(n√n)就可以吧。#include #include #include #include #include #include #definemaxn50010 usingnamespacestd; longlongf[maxn
u012288458·2016-04-01 08:00

【bzoj3930】[CQOI2015]选数 莫比乌斯反演

首先问题是从l~r中选n个数，最大公约数为k的方案数再转化一下，l/k~r/k中选n个数，最大公约数为1的方案数n=2时很容易，我们看一下n=3的情况那么，其实选n个数也是同理分块，枚举r/kd和l-1/kd的取值，然后快速幂计算就可以了，问题在于如何处理出μ函数的前缀和？还是参考PoPoQQQ大爷的吧：http://blog.csdn.net/popoqqq/article/details/44
qingdaobaibai·2016-03-31 14:13

【bzoj3930】[CQOI2015]选数 莫比乌斯反演

首先问题是从l~r中选n个数，最大公约数为k的方案数再转化一下，l/k~r/k中选n个数，最大公约数为1的方案数n=2时很容易，我们看一下n=3的情况那么，其实选n个数也是同理分块，枚举r/kd和l-1/kd的取值，然后快速幂计算就可以了，问题在于如何处理出μ函数的前缀和？还是参考PoPoQQQ大爷的吧：http://blog.csdn.net/popoqqq/article/details/44
u012288458·2016-03-31 14:00

【bzoj3529】[Sdoi2014]数表 线性筛法+树状数组+莫比乌斯反演+数论分块

后面那一块暴力做出来就可以了，复杂度是O(nlogn)的考虑a的限制，因为f(i) #include #include #include #include #include #definemaxn100010 usingnamespacestd; structyts { intn,m,a,id; }q[20100]; structyts1 { intid,data; }seq[maxn];
u012288458·2016-03-31 10:00

【bzoj2820】YY的GCD 线性筛法+莫比乌斯反演+数论分块

枚举每个质数，然后暴力算，TLE换一种思路，改变枚举顺序这样可以枚举n/T的取值，只需要预处理的前缀和就可以了。因为1~n中大概有n/lnn个质数，每个质数平均会更新lnn次，所以暴力处理即可，总复杂度为O(n)。#include #include #include #include #include #include #definemaxn10000010 usingnamespacestd;
u012288458·2016-03-30 19:00

【bzoj2301】[HAOI2011]Problem b 莫比乌斯反演+线性筛法+数论分块

首先可以容斥一下，问题变为1~n和1~m中有多少对数最大公约数为k，再转换一下就是1~n/k和1~m/k中有多少对互质的数，这个问题的答案就是最后那个式子。这个东西是可以分块做的，枚举n/kd的取值，预处理出μ(d)的前缀和即可。#include #include #include #include #include #include #definemaxn50010 usingnamespac
u012288458·2016-03-30 14:00

[置顶] 计划

不断更新中计划打的博客1、莫比乌斯反演复习小记2、把以后要注意的东西里面进行算法分类计划刷的题目1、【GDOI2103模拟3.17】数树数2、【SDOI2013】方程3、【SDOI2014】旅行4、【HNOI2013
doyouseeman·2016-03-30 12:00

bzoj 2820 莫比乌斯反演

链接：戳这里2820:YY的GCDTimeLimit:10Sec MemoryLimit:512MB[Submit][Status][Discuss]Description神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1 #include #include #include #include #include #include #include #include #include #incl
libin66·2016-03-29 16:00

hdu 1695 莫比乌斯反演

链接：戳这里GCDTimeLimit:6000/3000MS(Java/Others)  MemoryLimit:32768/32768K(Java/Others)ProblemDescriptionGiven5integers:a,b,c,d,k,you'retofindxina...b,yinc...dthatGCD(x,y)=k.GCD(x,y)meansthegreatestcommond
libin66·2016-03-28 21:00

bzoj 2301 莫比乌斯反演

链接：戳这里2301:[HAOI2011]ProblembTimeLimit:50Sec MemoryLimit:256MB[Submit][Status][Discuss]Description对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y)=k，gcd(x,y)函数为x和y的最大公约数。Input第一行一个整数n，接下来n行每行五个整数，分别表示
libin66·2016-03-28 20:00

rqnoj 607 莫比乌斯反演应用

链接：戳这里PID607/[NOI2010]能量采集☆题目标签类型数论/数值 来源NOI2001 题目描述栋栋有一块长方形的地，他在地上种了一种能量植物，这种植物可以采集太阳光的能量。在这些植物采集能量后，栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起。栋栋的植物种得非常整齐，一共有n列，每列有m棵，植物的横竖间距都一样，因此对于每一棵植物，栋栋可以用一个坐标(x,y)来表示，其中x的
libin66·2016-03-28 18:00

bzoj 2440 莫比乌斯反演应用

链接：戳这里2440:[中山市选2011]完全平方数TimeLimit:10Sec MemoryLimit:128MB[Submit][Status][Discuss]Description小X自幼就很喜欢数。但奇怪的是，他十分讨厌完全平方数。他觉得这些数看起来很令人难受。由此，他也讨厌所有是完全平方数的正整数倍的数。然而这丝毫不影响他对其他数的热爱。 这天是小X的生日，小W想送一个数给他作为生日
libin66·2016-03-28 15:00

BZOJ-2154 && 2693 Crash的数字表格&&jzptab 莫比乌斯反演

2154:Crash的数字表格TimeLimit:20SecMemoryLimit:259MBSubmit:2174Solved:816[Submit][Status][Discuss]Description今天的数学课上，Crash小朋友学习了最小公倍数(LeastCommonMultiple)。对于两个正整数a和b，LCM(a,b)表示能同时被a和b整除的最小正整数。例如，LCM(6,8)=2
DaD3zZ·2016-03-24 19:55

莫比乌斯函数筛法 & 莫比乌斯反演

模板：intp[MAXN],pcnt=0,mu[MAXN]; boolnotp[MAXN]; voidshai(intn){ mu[1]=1; for(inti=2;i<=n;++i){ if(notp[i]==0){ p[++pcnt]=i; mu[i]=-1; } for(intj=1,t=p[j]*i;j<=pcnt&&t<=n;++j,t=p[j]*i){ notp[t
abclzr·2016-03-21 19:00

3529: [Sdoi2014]数表|莫比乌斯反演|树状数组

以下式子均设n≤m先设M(x)=∑d|xd在不考虑a的限制下Ans=∑i=1n∑j=1nM(gcd(i,j))=∑i=1nM(i)∗sum(i)sum(i)表示gcd(x,y)=i的(x,y)的个数显然这个反演一下就可以得到sum(i)=∑x=1⌊ni⌋u(i)∗⌊ni∗x⌋∗⌊mi∗x⌋Ans=∑i=1nM(i)∗∑x=1⌊ni⌋u(i)∗⌊ni∗x⌋∗⌊mi∗x⌋=∑T=1n⌊nT⌋∗⌊mT⌋
ws_yzy·2016-03-17 20:00

BZOJ3739: DZY loves math VIII

莫比乌斯反演可以推出来不会Markdown。。。。
liutian429073576·2016-03-17 20:00

bzoj2154 Crash的数字表格

    首先得到朴素递推式，然后运用莫比乌斯反演得到：    然后可以发现后面的值之和N/d和M/d有关，只有O(N^0.5)种；而后面的值同样可以分段O(N^0.5)解决，因此时间复杂度O(N)。
lych_cys·2016-03-16 09:00

hdu 1695 GCD 莫比乌斯反演入门

GCD题意:输入5个数a,b,c,d,k;(a=c=1,0 #include #include #include #include #include #include #include #include #include #include #include #include usingnamespacestd; #definerep0(i,l,r)for(inti=(l
hxer·2016-03-15 23:00

POJ 1394 (莫比乌斯反演)

SkyCodeTimeLimit: 1000MS MemoryLimit: 65536KTotalSubmissions: 1964 Accepted: 629DescriptionStanculikesspacetravelsbutheisapoorsoftwaredeveloperandwillneverbeabletobuyhisownspacecraft.Thatiswhyheisprep
morejarphone·2016-03-11 20:00

莫比乌斯反演学习笔记

莫比乌斯反演的形式：另一种描述是：一种是和所有的约数有关一种是和所有的倍数有关，解题的时候要根据题目选择合适的表达形式，感觉第二种用的比较多。
morejarphone·2016-03-09 20:00

[bzoj2956]模积和

学过莫比乌斯反演的都知道，对于任意i，⌊ni⌋只有根号n种取值。然后本题就可做了。#include #include #include #definefo(i,a,b)for(i=a;i<=b;i++
WerKeyTom_FTD·2016-03-09 19:00

bzoj 2154 Crash的数字表格（莫比乌斯反演及优化）

 Description今天的数学课上，Crash小朋友学习了最小公倍数(LeastCommonMultiple)。对于两个正整数a和b，LCM(a,b)表示能同时被a和b整除的最小正整数。例如，LCM(6,8)=24。回到家后，Crash还在想着课上学的东西，为了研究最小公倍数，他画了一张N*M的表格。每个格子里写了一个数字，其中第i行第j列的那个格子里写着数为LCM(i,j)。一个4*5的表格
hahalidaxin·2016-03-07 17:00

bzoj 3529 [Sdoi2014]数表（莫比乌斯反演+BIT）

 Description   有一张N×m的数表，其第i行第j列（1 2#include 3#include 4usingnamespacestd; 5 6constintN=1e5+10; 7 8structNode{ 9intn,m,a,id; 10booloperatorm)swap(n,m); 69que[i]=(Node){n,m,a,i}; 70mx=max(mx
hahalidaxin·2016-03-07 11:00

bzoj 2820 YY的GCD（莫比乌斯反演）

 Description神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1 2#include 3usingnamespacestd; 4 5typedeflonglongll; 6constintN=1e7+5; 7 8llmu[N],sum[N],su[N],sz,np[N]; 9 10voidget_mu() 11{ 12inti,j; 13mu[1]=1; 1
hahalidaxin·2016-03-07 10:00

bzoj 2301 [HAOI2011]Problem b（莫比乌斯反演）

 Description 对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y)=k，gcd(x,y)函数为x和y的最大公约数。Input第一行一个整数n，接下来n行每行五个整数，分别表示a、b、c、d、k Output共n行，每行一个整数表示满足要求的数对(x,y)的个数 SampleInput22515115152SampleOutput143HINT
hahalidaxin·2016-03-06 21:00

莫比乌斯反演

友情提醒：这篇文章中的大部分东西都出自popoqqq的课件《莫比乌斯反演》和hzwer的博客，orz首先我们来看一个函数，。这个函数还是十分常见的。
fjzzq2002·2016-03-04 23:00

2693: jzptab|莫比乌斯反演

作为公式恐惧症晚期患者..就继续写一发题解补救一下题目让求∑i=1n∑j=1mlcm(i,j)=∑i=1n∑j=1mi∗jgcd(i,j)可以枚举gcdAns=∑d=1n∑i=1n∑j=1mi∗jd(gcd(i,j)=d)=∑d=1n∑i=1⌊nd⌋∑j=1⌊md⌋d2∗i∗jd(gcd(i,j)=1)=∑d=1nd∑i=1⌊nd⌋∑j=1⌊md⌋i∗j(gcd(i,j)=1)然后只需要继续化简这
ws_yzy·2016-03-04 10:00

UVA 11014Make a Crystal

我们要求的就是g(1)啦，这就是莫比乌斯反演啦。(3)的情况只要令f(x)=[n/x]*[n/x]即可。
Fsss_7·2016-03-03 23:00

BZOJ 2154: Crash的数字表格

Orzwty大爷我冬令营的时候一定没学过莫比乌斯反演，没错，就是没学过。不然我怎么会什么都不记得了TAT然后就默默地复习了一遍，顺便写了这道题。PPT里最后一步实在推不下去了。。。。。。
nlj1999·2016-02-28 14:00

hdu 1695 莫比乌斯反演

GCDTimeLimit:6000/3000MS(Java/Others)    MemoryLimit:32768/32768K(Java/Others)TotalSubmission(s):8482    AcceptedSubmission(s):3152ProblemDescriptionGiven5integers:a,b,c,d,k,you'retofindxina...b,yinc.
zjw6463·2016-02-27 19:00

【bzoj3994】[SDOI2015]约数个数和 莫比乌斯反演

Description设d(x)为x的约数个数，给定N、M，求Input输入文件包含多组测试数据。第一行，一个整数T，表示测试数据的组数。接下来的T行，每行两个整数N、M。OutputT行，每行一个整数，表示你所求的答案。SampleInput2 74 56 SampleOutput110 121 HINT1 #include #include #include usingnamespace
LOI_DQS·2016-02-23 21:00

BZOJ 3994: [SDOI2015]约数个数和|莫比乌斯反演

有一个非常神的东西:d(n∗m)=∑i|n∑j|mgcd(i,j)==1并不知道怎么证明……然后有了这个式子后面就很简单了！然后Ans=∑i=1n∑j=1m⌊ni⌋⌊mj⌋(gcd(i,j)==1)=∑i=1n∑j=1m∑d|i,d|ju(d)⌊ni⌋⌊mj⌋=∑d=1min（n,m)u(d)∑i=1⌊nd⌋∑j=1⌊md⌋⌊nd∗i⌋⌊md∗j⌋=∑d=1min（n,m)u(d)(∑i=1⌊nd
ws_yzy·2016-02-21 09:00

还有哪些东西要写눈_눈

눈_눈感觉数论写不完先记下来哪些要写  容斥原理抽屉原理置换群与Polya定理母函数二项式反演莫比乌斯反演积性函数狄利克雷卷积偏序关系理论高斯消元概率与期望问题博弈论之极大极小过程博弈论之Nim鸽巢原理
镜外之主·2016-02-19 19:00

BZOJ2190: [SDOI2008]仪仗队

有两种做法一个是欧拉函数的前缀和另一个是莫比乌斯反演貌似数据水能让莫比乌斯反演过去我用的是莫比乌斯反演#include #include #include usingnamespacestd; constintmaxn
liutian429073576·2016-02-17 09:00

BZOJ4407: 于神之怒加强版

莫比乌斯反演发现搞出来的式子和昨天的差不多#include #include #include usingnamespacestd; const intmod=1000000007; const intmaxn
liutian429073576·2016-02-16 09:00

BZOJ 4407: 于神之怒加强版|莫比乌斯反演

不会搞数学公式很苦恼！！flag：会写数学公式之后一定好好写一发题解非常感谢龙爷（sd第一男选手！！可惜神犇都不写blog）提供线性筛做法2.16————————————————-一下均设n #include #include #include #include #include #include #include #include #include #defineN5000005 #define
ws_yzy·2016-02-15 20:00
上一页 14 15 16 17 18 19 20 21 下一页
按字母分类： ABCDEFGHIJKLMNOPQRSTUVWXYZ其他