BZOJ 3601 一个人的数论 ——莫比乌斯反演 高斯消元

http://www.cnblogs.com/jianglangcaijin/p/4033399.html

——lych_cys

我还是太菜了，考虑一个函数的值得时候，首先考虑是否积性函数，不行的话就强行展开，

如果是的话考虑最小因子的高次幂的情况

然后还要一点点猜想才行。

#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define ll long long
#define mp make_pair
#define md 1000000007

ll d,m;
ll h[110],a[110][110];

ll ksm(ll a,ll b)
{
    if (b<0) return ksm(ksm(a,md-2),-b);
    ll ret=1;
    while (b)
    {
        if (b&1) (ret*=a)%=md;
        (a*=a)%=md;
        b>>=1;
    }
    return ret;
}

ll pp[110],A[110];

void init()
{
    int i,j;
    for (i=0;i<=d+1;++i)
    {
        pp[i]=(ksm(i+1,d)+(i==0?0:pp[i-1]))%md;
        a[i][d+2]=pp[i];
        a[i][0]=1;
        ll pre=1;
        for (j=1;j<=d+1;++j)
        {
            (pre*=(i+1))%=md;
            a[i][j]=pre;
        }
    }
    int k;
    for (i=0;i<=d+1;++i)
    {
        for (j=i;j<=d+1;++j) if (a[j][i]) break;
        if (i!=j) for (k=0;k<=d+2;++k) swap(a[i][k],a[j][k]);
        for (j=0;j<=d+1;++j) if (j!=i&&a[j][i])
        {
            ll tmp=(a[j][i]*ksm(a[i][i],-1))%md;
            for (k=0;k<=d+2;++k) (a[j][k]-=tmp*a[i][k])%=md;
        }
    }
    for (i=0;i<=d+1;++i) A[i]=(a[i][d+2]*ksm(a[i][i],-1))%md;
}

ll dd[1220][3];

int main()
{
    scanf("%lld%lld",&d,&m);
    init();
    F(i,1,m) scanf("%lld%lld",&dd[i][0],&dd[i][1]);
    ll ans=0;
    F(i,0,d+1)
    {
        ll tmp=1;
        F(j,1,m)
        {
            tmp=tmp*ksm(dd[j][0],(ll)dd[j][1]*i)%md;
            tmp=tmp*(1-ksm(dd[j][0],d-i))%md;
        }
        ans+=A[i]*tmp%md;
        ans%=md;
    }
    ((ans%=md)+=md)%=md;
    printf("%lld\n",ans);
}