YY的GCD,洛谷P2257,莫比乌斯反演+狄利克雷卷积

正题

      题目要求这个东西

      \sum_{p\in prime} \sum_{i=1}^n \sum_{j=1}^m [gcd(i,j)=p]

      其实就是求\sum_{p\in prime} \sum_{i=1}^{floor(n/p)} \sum_{j=1}^{floor(m/p)} [gcd(i,j)=1]

      根据\sum_{d\mid n}\mu(d)=[n=1]

      换出来变成\sum_{p\in prime} \sum_{i=1}^{floor(n/p)} \sum_{j=1}^{floor(m/p)}\sum_{d\mid gcd(i,j)} \mu(d)

      枚举d变成\sum_{p\in prime} \sum_{d=1}^{floor(n/p)} \mu(d) \frac{n}{dp}\frac{m}{dp}

      设T=dp,那么就变成\sum_{T=1}^n\frac{n}{T}\frac{m}{T}\sum_{p\in prime\ dp=T} \mu(d)

      发现做不了。

      看一下后面的东西,我们把它设成G(T)

      那么G(T)=\sum_{p\in prime \ dp=T}\mu(d)

      很明显发现G=prime*\mu

      prime函数指的是第i位是否为质数,是则为1,否则为0.

      然后整除分块一下,预处理狄利克雷卷积。

      然后n才1e7,按定义枚举O(n\ln n+T\sqrt n)

      

#include
#include
#include
#include
using namespace std;

int T;
const int maxn=1e7;
int mu[maxn+10],q[maxn+10],P[maxn+10];
bool vis[maxn+10];
int n,m;

int main(){
	scanf("%d",&T);
	mu[1]=1;vis[1]=true;
	for(int i=2;i<=maxn;i++){
		if(!vis[i]) {P[++P[0]]=i;mu[i]=-1;}
		for(int j=1;j<=P[0] && (long long) i*P[j]<=maxn;j++){
			int temp=i*P[j];
			vis[temp]=true;
			if(i%P[j]==0) break;
			mu[temp]=-mu[i];
		}
	}
	for(int i=1;i<=P[0];i++)
		for(int j=1;j<=maxn/P[i];j++) q[P[i]*j]+=mu[j];
	for(int i=1;i<=maxn;i++) q[i]+=q[i-1];
	while(T--){
		scanf("%d %d",&n,&m);
		if(m

 

你可能感兴趣的:(YY的GCD,洛谷P2257,莫比乌斯反演+狄利克雷卷积)