【题解】Zap(莫比乌斯反演)

【题解】Zap(莫比乌斯反演)

裸题...

直接化吧

[P3455 POI2007]ZAP-Queries

所有除法默认向下取整

[Sigma_{i=1}^xSigma_{j=1}^y[(i,j)=k] \ =Sigma_{i=1}^{x/k}Sigma_{j=1}^{y/k}[(i,j)=1] \ =Sigma_{i=1}^{x/k}Sigma_{j=1}^{y/k}Sigma_{d|(i,j)}mu(d) \ =Sigma_{d=1}^{min(x,y)}Sigma_{i=1}^{x/k}Sigma_{j=1}^{y/k}mu(d) imes[d|(i,j)] \ =Sigma_{d=1}^{min(x,y)}(frac x {dk})(frac y {dk})mu(d) ]

整除分块直接做...

有一个细节，可能有疑惑：

		r=min(x/(x/l),y/(y/l));
		ans+=1ll*(x/(l*k))*(y/((l*k)))*(sum[r]-sum[l-1]);

整除分块为什么是这样的？为什么r=min(x/(x/l),y/(y/l));中的"(l)"和ans+=1ll*(x/(l*k))*(y/((l*k)))*(sum[r]-sum[l-1]);不统一，为什么是(x/(l*k))*(y/(l*k))？这不是整除分块正常的套路啊？

可以这样理解，整除分块利用了(lfloor frac x l floor)在一定范围内不变的性质，所以我们同样也会有(lfloorfrac {lfloor frac x l floor} k floor)在一定范围内不变化，并且前面那个式子包括的(l)的范围一定小于后面的那个(l)的范围，所以我们按照(lfloor frac x l floor)整除分块即可。

至于如何按照(lfloorfrac {lfloor frac x l floor} k floor=lfloor frac x {lk} floor)分块，我也不知道怎么办，希望有高手指点一下QAQ

#include<bits/stdc++.h>

using namespace std;typedef long long ll;
template < class ccf >
 inline ccf qr(ccf b){
    register char c=getchar();register int q=1;register ccf x=0;
    while(c<48||c>57)q=c==45?-1:q,c=getchar();
    while(c>=48&&c<=57)x=x*10+c-48,c=getchar();
    return q==-1?-x:x;}
inline int qr(){return qr(1);}
const int maxn=1e5+5;
bool usd[maxn];
int mu[maxn];
int sum[maxn];
vector < int > ve;
int x,y,k;
#define pb push_back
inline void gen(){
    mu[1]=sum[1]=usd[1]=1;
    for(register int t=2;t< maxn;++t){
	if(not usd[t])
	    ve.pb(t),mu[t]=-1;
	for(register auto p:ve)
	    if(1ll*p*t<maxn)
		if(usd[p*t]=1,t%p) mu[p*t]=-mu[t];
		else break;
	    else break;
	sum[t]=sum[t-1]+mu[t];
    }
}

int main(){
    gen();
    int T=qr();
    while(T--){
	x=qr();y=qr();k=qr();
	ll ans=0;
	for(register int l=1,r=0,edd=min(x,y)/k;l<=edd;l=r+1){
	    r=min(x/(x/l),y/(y/l));
	    ans+=1ll*(x/(l*k))*(y/((l*k)))*(sum[r]-sum[l-1]);
	}
	cout<<ans<<endl;
    }
    return 0;
}