莫比乌斯反演初步2:BZOJ 2301: [HAOI2011]Problem b

题目描述:戳这里
题解:
这题看似和HDU1965很像,其实就是很像。(HDU1965题解
反演的过程是一样的,只是这题的数据更加大一点。所以我们一定要用除法分块。
除法的原理可以看这个blog,讲得很清楚。
还有一点,这题的a和c不是1,所以要容斥一下:
答案是:mu(b,d)-mu(a-1,d)-mu(b,c-1)+mu(a-1,c-1)

代码如下:

#include<cstdio>
#include<string>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=50005;
int T,m,a,b,c,d,k,su[maxn],f[maxn];
bool vis[maxn];
void make_p(){
    f[1]=vis[1]=1;
    for (int i=2;i<=50000;i++){
        if (!vis[i]) su[++m]=i,f[i]=-1;
        for (int j=1;j<=m&&i*su[j]<=50000;j++) {
            vis[i*su[j]]=1;
            if (i%su[j]==0) {f[i*su[j]]=0; break;}
            f[i*su[j]]=-f[i];
        }
    }
    for (int i=2;i<=50000;i++) f[i]+=f[i-1];
}
ll getmu(int x,int y){
    x/=k; y/=k;
    int minn=min(x,y);
    int r=0; ll ret=0;
    for (int i=1;i<=minn;i=r+1) {
        r=min(x/(x/i),y/(y/i));
        ret+=1ll*(f[r]-f[i-1])*(x/i)*(y/i);
    }
    return ret;
}
int main(){
    make_p();
    scanf("%d",&T);
    while (T--) {
        scanf("%d %d %d %d %d",&a,&b,&c,&d,&k);
        printf("%lld\n",getmu(b,d)-getmu(a-1,d)-getmu(b,c-1)+getmu(a-1,c-1));
    }
    return 0;
}
发布了141 篇原创文章 · 获赞 72 · 访问量 1万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 深蓝海洋 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览