#include<bits/stdc++.h> #define N 5000010 #define yql 1000000007 using namespace std; typedef long long ll; ],m[],maxn,vis[N],prime[N]; ,fac[N]; ll f[N],s[N]; inline int read(){ ,x=;char ch; ;}'); +ch-'); return f*x; } inline ll fpow(ll x,ll p){ x%=…

题解 非常裸的莫比乌斯反演. 但是反演完还需要快速计算一个积性函数(我直接用$nlogn$卷积被TLE了 推荐一个博客 我也不想再写一遍了 代码 #include<cstring> #include<cstdio> #include<algorithm> #define ll long long #define rd read() using namespace std; ; const int N = 5e6; int T, n, m, k; int mu[N], p…

Description 给下N,M,K.求 感觉好迷茫啊,很多变换看的一脸懵逼却又不知道去哪里学.一道题做一上午也是没谁了,, 首先按照套路反演化到最后应该是这个式子 $$ans = \sum_{d = 1}^n d^k \sum_{i = 1}^{\frac{n}{d}} \frac{n}{di} \frac{m}{di} \mu(i)$$ 这样就可以$O(n)$计算 继续往下推,考虑$\frac{n}{di} \frac{m}{di}$对答案的贡献 设$T = id$ $ans = \sum…

现在感觉反演好多都是套路QAQ…… #include<bits/stdc++.h> using namespace std; ; typedef long long ll; int n,cnt,prime[N],phi[N],mu[N],vis[N]; ll ans,s[N],f[N]; void calcmu(){ memset(prime,,; memset(phi,,,sizeof(mu)); memset(s,,,sizeof(f)); mu[]=;phi[]=;memset(vis,…

也是一个反演. 第一次手动推出一个简单的式子,激动.jpg 原题意思是求:$Ans=\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}\phi(gcd(i,j))$随意化几步试试:$Ans=\sum\limits_{d=1}^{min(n,m)}\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}\phi(d)[gcd(i,j)==d]$$Ans=\sum\limits_{d=1}^{min(n,m)}\sum\limits_{i…

快三个月没做反演题了吧…… 感觉高一上学期学的全忘了…… 所以还得从零开始学推式子. # bzoj1011 标签(空格分隔): 未分类 --- 原题意思是求以下式子:$Ans=\sum\limits_{i=1}^{a}\sum\limits_{i=1}^{b}[gcd(i,j)==k]$首先把k拿下来,得到$Ans=\sum\limits_{i=1}^{a/k}\sum\limits_{i=1}^{b/k}[gcd(i,j)==1]$然后考虑mobius函数的性质:$\sum\limits_{d…

BZOJ挂了... 先把程序放上来,如果A了在写题解吧. #include<cstdio> #include<algorithm> #define N 5000010 #define ll long long #define mod (int)(1e9+7) using namespace std; inline int read() { ,f=;char ch=getchar(); ;ch=getchar();} +ch-';ch=getchar();} return x*f;…

同bzoj3930. (日常盗题图) #include<bits/stdc++.h> #define N 1000010 #define yql 1000000007 #define ll long long using namespace std; int n; ll m,l,r,k;ll f[N]; ll Pow(ll x,int y){ ll ans=; ) ans=ans*x%yql; x=x*x%yql; y>>=; } return ans; } ll read(){…

出题人怎么这么不认真啊==明明官方译名是爱蜜莉雅…… 而且我们爱蜜莉雅碳是有英文名哒!是Emilia.你那个aimiliya我实在是无力吐槽…… 不过抱图跑23333首先这很像约数个数和函数诶!但是唯一的不同是,D(1)=1,F(1)=0.那么如果就是D,我们怎么做?原题意思是求以下式子:$Ans\sum\limits_{i=1}^{n}d(gcd(i,n))$我们知道:$D(n)=\sum\limits_{d|n}1$所以:$Ans=\sum\limits_{d|n}\sum\limits_{…

膜拜cdc……他的推导详细到我这种蒟蒻都能看得懂! 膜拜的传送门 所以我附一下代码就好了. #include<bits/stdc++.h> #define N 10000005 #define yql 20101009 using namespace std; typedef long long ll; int mu[N],prime[N],cnt,s[N],vis[N]; ll n,m,ans,maxn; inline ll sum(ll x,ll y){ )/)%yql)*((y*(y+)…