BZOJ 2005 [Noi2010]能量采集(严格来说应该叫欧拉反演)

最新推荐文章于 2019-07-08 21:39:54 发布

原创最新推荐文章于 2019-07-08 21:39:54 发布 · 195 阅读

1 ·

CC 4.0 BY-SA版权

数论同时被 3 个专栏收录

157 篇文章

订阅专栏

性质分析

86 篇文章

订阅专栏

反演

2 篇文章

订阅专栏

本文探讨了使用LaTeX表达复杂数学公式的方法，并深入解析了一种数论算法，该算法涉及欧拉函数、质数筛法及整除分块技术。通过实际代码展示了算法的实现过程，对理解数学公式与算法设计具有一定的指导意义。

练习一发LaTeX
$ans=∑x=1n∑y=1m2gcd(x,y)−1=∑x=1n∑y=1m(−1+2∑p∣gcd(x,y)φ(p))=2∑x=1n∑y=1m∑p∣gcd(x,y)φ(p)−n∗m=2∑p=1min(n,m)φ(p)∑x=1n∑y=1m[p∣x][p∣y]−n∗m=2∑p=1min(n,m)φ(p)⌊np⌋⌊mp⌋−n∗m\begin{aligned}ans &= \sum_{x=1}^n\sum_{y=1}^m2gcd(x,y)-1\\ &=\sum_{x=1}^n\sum_{y=1}^m\left(-1+2\sum_{p|gcd(x,y)}\varphi(p)\right)\\ &=2\sum_{x=1}^n\sum_{y=1}^m\sum_{p|gcd(x,y)}\varphi(p) - n*m\\ &=2\sum_{p=1}^{min(n,m)}\varphi(p)\sum_{x=1}^n\sum_{y=1}^m[p|x][p|y] - n*m\\ &=2\sum_{p=1}^{min(n,m)}\varphi(p)\lfloor\frac np\rfloor \lfloor\frac mp\rfloor - n*m\end{aligned}$
这个搁今天得整除分块后杜教筛1e10+多组数据

AC Code:

#include<bits/stdc++.h>
#define maxn 100005
#define LL long long
using namespace std;

int n,m,phi[maxn],pr[maxn/8],cnt_pr;
bool vis[maxn];

int main()
{
	scanf("%d%d",&n,&m); if(n>m) swap(n,m);
	phi[1] = 1;
	LL ans = 1ll*n*m;
	for(int i=2;i<=n;i++)
	{
		if(!vis[i]) pr[cnt_pr++] = i , phi[i]=i-1;
		for(int j=0;pr[j]*i<=n;j++)
		{
			vis[pr[j] * i] = 1;
			if(i % pr[j] == 0){ phi[i*pr[j]]=phi[i]*pr[j];break; }
			phi[i*pr[j]] = phi[i] * phi[pr[j]];
		}
		ans += 1ll*(2*phi[i])*(n/i)*(m/i);
	}
	printf("%lld\n",ans);
}

upd:自己杠精了。
$∑i=1n∑j=1njφ(i)=∑j=1n∑i=1njφ(i)=∑j=1nΦ(nj)=∑i=1ni=n(n+1)2\begin{aligned} \sum_{i=1}^n\sum_{j=1}^{\frac nj} \varphi(i) = \sum_{j=1}^n\sum_{i=1}^{\frac nj} \varphi(i) = \sum_{j=1}^n \Phi(\frac nj) = \sum_{i=1}^n i = \frac {n(n+1)}2 \end{aligned}$
AC Code:

#include<bits/stdc++.h>
#define maxn 100005
#define LL long long
using namespace std;

int n,m,phi[maxn],pr[maxn/8],cnt_pr;
bool vis[maxn];
map<int,LL>mp;

LL solve(int now)
{
	if(now <= 2137) return phi[now];
	if(mp.count(now)) return mp[now];
	LL &ret=mp[now]=1ll*now*(now+1)/2;
	for(int i=2,nxt;i<=now;i=nxt+1)
	{
		int tmp = now/i;
		nxt = now / tmp;
		ret -= solve(tmp) * (nxt-i+1);
	}
	return ret;
}

int main()
{
	scanf("%d%d",&n,&m); if(n>m) swap(n,m);
	phi[1] = 1;
	LL ans = -1ll*n*m;
	for(int i=2;i<=min(n,2137);i++)
	{
		if(!vis[i]) pr[cnt_pr++] = i , phi[i]=i-1;
		for(int j=0;pr[j]*i<=n;j++)
		{
			vis[pr[j] * i] = 1;
			if(i % pr[j] == 0){ phi[i*pr[j]]=phi[i]*pr[j];break; }
			phi[i*pr[j]] = phi[i] * (pr[j]-1);
		}
		phi[i] += phi[i-1];
	}
	for(int i=1,nxt,psolve=0,tmp,div1,div2;i<=n;i=nxt+1)
	{
		div1 = n / i , div2 = m / i;
		nxt = min(n/div1,m/div2);
		ans += 2ll * ((tmp=solve(nxt))-psolve) * div1 * div2;
		psolve = tmp;
	}
	printf("%lld\n",ans);
}