看题解一开始还有地方不理解,果然是我的组合数学思维比较差
然后理解了之后自己敲了一个果断TLE。。。。
我以后果然还得多练啊
好巧妙的思路啊
知识1:
费马小定理是数论中的一个重要定理,其内容为:假如p是质数,且(a,p)=1,那么 a^(p-1) ≡1(mod p)假如p是质数,且a,p互质,那么 a的(p-1)次方除以p的余数恒等于1
对于除法取模还需要用到费马小定理: a ^ (p - 1) % p = 1; -> a ^ (p - 2) % p = (1 / a) % p;
巧妙1:
for(int i=1;i<=n;i++)
{ int temp; scanf("%d",&temp); sum1[temp]++; }
for(int j=i;j<=m;j+=i) sum+=sum1[j];
直接判断倍数是否有无。ORZ!!!
用这一块代码,这样再从1遍历到m的时候,速度增加非常快,然后就不会超时。
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <string>
#include <queue>
#include <set>
#include <map>
#include <vector>
#include <assert.h>
using namespace std;
#define lowbit(i) (i&-i)
#define sqr(x) ((x)*(x))
#define enter printf("\n")
#define is_sqr(x) (x&(x-1))
#define pi acos(-1.0)
#define clr(x) memset(x,0,sizeof(x))
#define fp1 freopen("in.txt","r",stdin)
#define fp2 freopen("out.txt","w",stdout)
#define pb push_back
typedef __int64 LL;
const double eps = 1e-7;
const double DINF = 1e100;
const int INF = 1000000006;
const LL LINF = 1000000000000000005ll;
const int MOD = (int) 1e9 + 7;
const int maxn=300005;
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
template<class T> inline T Min(T a,T b,T c){return min(min(a, b),c);}
template<class T> inline T Max(T a,T b,T c){return max(max(a, b),c);}
LL f[maxn],e[maxn],a[maxn],ans[maxn],sum1[maxn];
LL quick_pow(LL a,LL b)//a的b次方,快速幂取模
{
LL ret=1;
while(b)
{
if(b&1) ret=(ret*a)%MOD;
b/=2;
a=(a*a)%MOD;
}
return ret%MOD;
}
LL cal(LL n,LL k)
{
if(k==0||n==k) return 1;
return (f[n]*e[k]%MOD)*e[n-k]%MOD;//注意运算顺序
}
//以后某些变量还是不采用C99的写法了
int main()
{
f[0]=e[0]=1;
for(int i=1;i<=maxn;i++)
{
f[i]=f[i-1]*i%MOD;
e[i]=quick_pow(f[i],MOD-2);
}
int n,m,k;
while(scanf("%d%d%d",&n,&m,&k)!=EOF)
{
clr(sum1);
for(int i=1;i<=n;i++)
{
int temp;
scanf("%d",&temp);
sum1[temp]++;
}
for(int i=m;i>=1;i--)//倒着写不至于每次都m/i次循环
{
int sum=0;
for(int j=i;j<=m;j+=i)
sum+=sum1[j];
if(sum<n-k)//k个不一样的,n-k个一样的。
{
ans[i]=0;continue;
}
ans[i]=((cal(sum,n-k)*quick_pow(m/i-1,sum-(n-k))%MOD)*quick_pow(m/i,n-sum))%MOD;
for(int j=2*i;j<=m;j+=i)
ans[i]=(ans[i]-ans[j]+MOD)%MOD;
}
for(int i=1;i<m;i++) printf("%lld ",ans[i]);
printf("%lld\n",ans[m]);
}
return 0;
}