HDU4675 莫比乌斯反演

最新推荐文章于 2019-10-15 11:16:05 发布

Gipsyu

最新推荐文章于 2019-10-15 11:16:05 发布

阅读量252

点赞数

分类专栏：数论组合数学

本文链接：https://blog.csdn.net/Gipsy_Danger/article/details/83002385

版权

数论同时被 2 个专栏收录

49 篇文章 0 订阅

订阅专栏

组合数学

11 篇文章 0 订阅

订阅专栏

卡了半天原来是卡在预处理上了

以为n根号n的预处理不碍大事

后来发现可以用nlog来处理

当时没想到。。

其实原理很简单。。。

看了一眼别人的代码

其他的都是自己推的

一道简单的组合数加中档莫比乌斯的题

其实莫比乌斯的本质就是容斥

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=3e5,M=1e9+7;
ll U[N+10],V[N+10],A[N+10],B[N+10],n,m,k,JC[N+10],INV[N+10],FF[N+10];
ll C(ll a,ll b,ll m){
    return JC[a]*INV[b]%m*INV[a-b]%m;
}
ll quickmod(ll a,ll b,ll m){
    ll ans=1;
    for(;b;b>>=1){
        if(b&1)ans=ans*a%m;
        a=a*a%m;
    }
    return ans;
}
ll F(ll x){
	if(B[x]<n-k)return 0;
    return C(B[x],n-k,M)*quickmod(m/x,n-B[x],M)%M*quickmod(m/x-1,B[x]-n+k,M)%M;
}
inline void read(ll &x){
	char ch;bool ok;
    for(ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if(ch=='-') ok=1;
    for(x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
    if(ok) x=-x;
}
int main(){
	//freopen("D://abc.txt","r",stdin);
    //freopen("D://1.txt","w",stdout);
    JC[0]=1,INV[0]=1;
    for(ll i=1;i<=N;++i)JC[i]=JC[i-1]*i%M,INV[i]=quickmod(JC[i],M-2,M);
    for(int i=1;i<=N;++i)U[i]=1;
    for(int i=2;i<=N;++i){
        if(V[i])continue;
        for(int j=i;j<=N;j+=i){
            V[j]=1;
            if(j%(i*i)==0)U[j]=0;
            else U[j]*=-1;
        }
    }
    while(cin>>n>>m>>k){
    	memset(B,0,sizeof B);
	    for(int i=1;i<=n;++i){
	        read(A[i]);
	        ++B[A[i]];
	    }
	    for(int i=1;i<=m;++i){
	    	for(int j=i*2;j<=m;j+=i){
	    		B[i]+=B[j];
			}
		}
	    for(int i=1;i<=m;++i)FF[i]=F(i);
	    for(int i=1;i<=m;++i){
	        ll ans=0;
	        for(int j=i;j<=m;j+=i){
	            ans=(ans+(U[j/i]*FF[j])%M+M)%M;
	        }
	    	printf("%lld",ans);
	    	if(i!=m)printf(" ");
	    	else printf("\n");
	    }
	}
}