BZOJ_P4407 于神之怒加强版(数论+莫比乌斯反演)

最新推荐文章于 2020-03-30 16:06:39 发布

北屿小智障

最新推荐文章于 2020-03-30 16:06:39 发布

阅读量576

点赞数

分类专栏： BZOJ 数论莫比乌斯反演

本文链接：https://blog.csdn.net/qq_18455665/article/details/50813543

版权

BZOJ 同时被 3 个专栏收录

129 篇文章 0 订阅

订阅专栏

数论

54 篇文章 0 订阅

订阅专栏

莫比乌斯反演

6 篇文章 0 订阅

订阅专栏

BZOJ传送门
Time Limit: 80 Sec Memory Limit: 512 MB
Submit: 195 Solved: 94
[Submit][Status][Discuss]
Description
给下N,M,K.求
这里写图片描述

Input
输入有多组数据，输入数据的第一行两个正整数T,K，代表有T组数据，K的意义如上所示，下面第二行到第T+1行，每行为两个正整数N,M，其意义如上式所示。

Output
如题

Sample Input
1 2
3 3

Sample Output
20

HINT
1<=N,M,K<=5000000,1<=T<=2000

Source
命题人：成都七中张耀楠，鸣谢excited上传。

这里写图片描述

#include<cstdio>
#include<iostream>
using namespace std;
#define N 5000010
#define P 1000000007
#define min(a,b) a>b?b:a
inline int in(int x=0,char ch=getchar(),int v=1){
    while(ch!='-'&&(ch>'9'||ch<'0')) ch=getchar();if(ch=='-') v=-1,ch=getchar();
    while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();return x*v;}
int f[N],pr[N],pw[N];bool b[N];
int cnt,n,m,k,t,ans;
inline int Pow(int a,int b,int res=1){while(b){if(b&1) res=(long long)res*a%P;a=(long long)a*a%P,b>>=1;}return res;}
inline void GetMu(){
    f[1]=1;
    for(int i=2;i<N;i++){
        if(!b[i]) pr[cnt]=i,pw[cnt]=Pow(i,k),f[i]=pw[cnt]-1,cnt++;
        for(int j=0;j<cnt&&i*pr[j]<N;j++){
            b[i*pr[j]]=1;
            if(!(i%pr[j])){f[i*pr[j]]=(long long)f[i]*pw[j]%P;break;}   
            f[i*pr[j]]=(long long)f[i]*f[pr[j]]%P;
        }
    }
    for(int i=2;i<N;i++)f[i]=(f[i]+f[i-1])%P;
}
int main(){
    t=in(),k=in();GetMu();
    while(t--){
        n=in(),m=in();if(n>m) swap(n,m);ans=0;
        for(int i=1,j=0;i<=n;i=j+1){
            j=min(n/(n/i),m/(m/i));
            ans=((long long)(f[j]-f[i-1]+P)*(n/i)%P*(m/i)%P+ans)%P;
        }
        printf("%d\n",ans);
    }
    return 0;
}