[BZOJ 2301][HAOI 2011]Problemb(数論)

タイトルの説明
トランスファゲート
問題解
基本的にBZOJ 1101と同じです.しかし、ここには上下の境界があり、それはすべて1を求めることに転化して、それから加減します.
コード#コード#

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=5e4;
int prime[N+5],mu[N+5],p[N+5];
int T,a,b,c,d,k,ans;
inline void get_mu(){
    mu[1]=1;
    for (int i=2;i<=N;++i){
        if (!p[i]){
            prime[++prime[0]]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=prime[0]&&i*prime[j]<=N;++j){
            p[i*prime[j]]=1;
            if (i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            }
            else mu[i*prime[j]]=-mu[i];
        }
    }
}
inline int calc_ans(int n,int m){
    if (n>m) swap(n,m);
    int ans=0;
    for (int i=1,j;i<=n;i=j+1){
        j=min(n/(n/i),m/(m/i));
        ans+=(n/i)*(m/i)*(mu[j]-mu[i-1]);
    }
    return ans;
}
int main(){
    get_mu();
    for (int i=1;i<=N;++i) mu[i]+=mu[i-1];
    scanf("%d",&T);
    while (T--){
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);a--,c--;
        a/=k; b/=k; c/=k; d/=k;
        ans=calc_ans(a,c)+calc_ans(b,d)-calc_ans(a,d)-calc_ans(b,c);
        printf("%d
",ans);
    }
}