2019.6.1BZOJ3028:食物生成函数题，母函数乘起来就好了BZOJ3544:[ONTAK2010]CreativeAccounting嗯，就是可以用set维护前缀和，取后继或最小数贪心就好啦BZOJ2820

题目：BZOJ2820.题目大意：求有多少对x,yx,yx,y满足1≤x≤n,1≤y≤m1\leqx\leqn,1\leqy\leqm1≤x≤n,1≤y≤m且gcd(x,y)gcd(x,y)gcd(x,y)为质数.1≤n,m≤1071\leqn,m\leq10^71≤n,m≤107，数据组数=104=10^4=104.设素数集为P，那么题目要求即为：∑p∈P∑i=1n∑j=1m[gcd(i,j)=

bzoj1101\2301：入门题，利用μ(i)前缀和的性质bzoj2440：不能算反演吧，只是利用了μ的性质bzoj2226：可以用跟欧拉函数有关的一个公式，也可以反演bzoj2820：基础练习题bzoj2154

题目链接：YY的GCD公式和上一题一样地推最后得到答案为令T=pd'，化简公式得到后面那个sigma可以线性筛筛出来，枚举素数贡献答案即可#include#include#include#include#include#include#include#include#include#include#defineLLlonglongusingnamespacestd;constintmaxn=100

点击这里查看原题经过一系列反演，得到ans(n)=∑T=1n⌊nT⌋∗⌊nT⌋∗∑i|Tφ(i)∗μ(Ti)（看不懂的建议先去做BZOJ2820）于是接下来就是要对∑i|Tφ(i)∗μ(Ti)求前缀和，

BZOJ2820-YY的GCD（莫比乌斯反演）题目链接BZOJ2820-YY的GCD题意TTT组查询，每次给定N,MN,MN,M，求1<=x<=N,1<=y<=M1<=x<=N,1<=y<=M1usingnamespacestd;typedeflonglongll;constintmaxn=10000005;intvis[maxn];intprime

DarkBZOJ传送门洛谷传送门解析：应教练要求开始准备数论讲义，现在把开通博客之前的写一些数论题目给弄上来。显然题目要求的是这个东西：∑i=1n∑j=1m[gcd(i,j)isaprime]\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)\mathrm{is\text{}a\text{}prime}]i=1∑nj=1∑m[gcd(i,j)isaprime]接下来以P\mathb

BZOJ2820莫比乌斯反演模板证明：莫比乌斯反演的两种形式：1.F(n)=∑d∣nf(d)⇒f(n)=∑d∣nμ(d)F(nd)F(n)=\sum_{d|n}{f(d)}\Rightarrowf(n)

题意给定N,M,求1#include#includeusingnamespacestd;constintmaxn=10000005,N=10000000;intn,m,_test,mu[maxn],p[maxn],sum[maxn];longlongans;boolvis[maxn];voidget_mu(){memset(vis,1,sizeof(vis));mu[1]=1;for(inti=2

莫比乌斯反演洛谷题目传送门BZOJ题目传送门和NOI2010能量采集很像。同样设f(x)=∑ni=1∑mj=1[(i,j)==x],F(x)=∑ni=1∑mj=1[x|(i,j)]=⌊nx⌋⌊mx⌋f(x)=∑i=1n∑j=1m[(i,j)==x],F(x)=∑i=1n∑j=1m[x|(i,j)]=⌊nx⌋⌊mx⌋设pp为质数，n#include#include#include#defineN10

======这里放传送门======题解记得当初懵（fei）懂（chang）无（sha）知（bi）的时候把这个题当一个无脑水题写了一发结果T的死惨死惨的。。。当时直接O(n)回答询问了也不想想能不能过。。然后现在这个题就变成了ATP的反演入门题_(:з」∠)_题目要求∑p为质数∑i=1N∑j=1M[(i,j)=p]设F(k)为满足1≤i≤N，1≤j≤M并且gcd(i,j)=k的数对个数，f(k)为

求gcd(i,j)为质数的个数。即∑p∑i=1n∑j=1mgcd(i,j)==p首先我们可以像bzoj2301一样把后一部分化简，得到(假定nusingnamespacestd;#definelllonglong#defineN10000010intT,n,m,mu[N],prime[N],tot=0;llf[N];boolnotprime[N];inlineintread(){intx=0;ch

传送门题意给x,y,求1≤i≤x,1≤j≤y且gcd(i,j)为质数的(i,j)对数。题解不妨设nusingnamespacestd;constintMaxn=1e7;typedeflonglongll;inlineintread(){charch=getchar();inti=0,f=1;while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}whil

题目描述传送门题解这题和上一题差不多的…令p(i)表示第i个质数，假设n#include#include#include#includeusingnamespacestd;#defineN10000005#defineLLlonglong#defineMod1000000007intT,n,m;intp[N],prime[N],mu[N],f[N];LLans;voidget(intn){mu[1

很好的一道莫比乌斯反演看了题解才把式子推出来去看ioi爷的博客吧http://www.cnblogs.com/iwtwiioi/p/4132095.html#include #include #include usingnamespacestd; const intmaxn=10000001; intmu[maxn],G[maxn]; longlongPresum[maxn]; boolcheck

mobius反演bzoj2820 yy的gcdbzoj3529sdoi2014  http://blog.csdn.net/iamzky/article/details/40376189神奇数论bzoj3884

2820:YY的GCDTimeLimit: 10Sec  MemoryLimit: 512MBSubmit: 1296  Solved: 672[Submit][Status][Discuss]Description神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1 #include #include #include #include #include #defineF(i,j,n)for(

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

YY的GCDDescription神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1 usingnamespacestd; #definerep(i,l,r)for(inti=(l);i=(l);i--) #defineXfirst #defineYsecond #defineMS(_)memset(_,0,sizeof(_)) #definePBpush_back #defineMPmake

链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2820题意：中文题。。。分析：同bzoj2301，论文题。代码：#include #include #include #include #include #include #include #include #include #include #include #pragmacomment(li

题目描述：求有多少数对(x,y)(1 #include #include usingnamespacestd; constintMAXN=10000000; boolnprime[MAXN+10]; intpcnt,prime[MAXN+10],sum[MAXN+10],mu[MAXN+10]; voidInit(intup=MAXN-1){ inttmp;mu[1]=1; for(inti=2;