版权声明:2333 https://blog.csdn.net/liangzihao1/article/details/84890025
题目大意:
给定
和
,求
其中,
,
,
为质数,
。
分析:
根据题意推式子得到
对于最后面的自然数幂和,可以通过高斯消元求出多项式,设多项式为
.
显然
,这一位就不用消了(避免出错)。
原式变成
把
移到前面,得到
我们可以枚举
,当
一定时,后面这些是一个积性函数,我们可以分别计算出每一个
的贡献,乘起来即可。而根据
的性质,只有
和
时有贡献,这个可以
算出。总复杂度是
的,另外高斯消元是
。
代码:
/**************************************************************
Problem: 3601
User: ypxrain
Language: C++
Result: Accepted
Time:812 ms
Memory:9240 kb
****************************************************************/
#include <iostream>
#include <cstdio>
#include <cmath>
#define LL long long
const int maxn=1007;
const LL mod=1e9+7;
using namespace std;
int n,k;
int p[maxn],q[maxn];
LL a[maxn][maxn],b[maxn],r[maxn];
LL sum,ans,num;
LL power(LL x,LL y)
{
if (y==1) return x;
if (y==0) return 1;
LL c=power(x,y/2);
c=(c*c)%mod;
if (y&1) c=(c*x)%mod;
return c;
}
void guass(int n)
{
for (int i=1;i<=n;i++)
{
for (int j=i+1;j<=n;j++)
{
if (a[j][i]>a[i][i])
{
for (int k=i;k<=n;k++) swap(a[i][k],a[j][k]);
swap(b[i],b[j]);
}
}
LL inv=power(a[i][i],mod-2);
for (int j=1;j<=n;j++)
{
if (i==j) continue;
LL rate=a[j][i]*inv%mod;
for (int k=i;k<=n;k++) a[j][k]=(a[j][k]+mod-a[i][k]*rate%mod)%mod;
b[j]=(b[j]+mod-b[i]*rate%mod)%mod;
}
}
for (int i=1;i<=n;i++) r[i]=b[i]*power(a[i][i],mod-2)%mod;
}
int main()
{
scanf("%d%d",&k,&n);
for (int i=1;i<=n;i++)
{
scanf("%d%d",&p[i],&q[i]);
q[i]=power(p[i],q[i]);
}
for (int i=1;i<=k+1;i++) a[i][1]=i;
for (int j=2;j<=k+1;j++)
{
for (int i=1;i<=k+1;i++) a[i][j]=(a[i][j-1]*(LL)i)%mod;
}
for (int i=1;i<=k+1;i++)
{
sum=(sum+power(i,k))%mod;
b[i]=sum;
}
guass(k+1);
for (int i=1;i<=k+1;i++)
{
num=1;
for (int j=1;j<=n;j++)
{
LL e=q[j]*power(p[j],mod-2)%mod;
num=num*((power(q[j],i)+mod-power(p[j],k)*power(e,i)%mod)%mod)%mod;
}
ans=(ans+r[i]*num%mod)%mod;
}
printf("%lld\n",ans);
}