HDU 6372 2018HDU多校赛 第六场 sacul（Lucas定理+矩阵）

 

 

大致题意：HMBB(海绵宝宝矩阵)的长宽为，p表示第c个质数。然后每一个位置(i,j)，如果 那么这个位置为1，否则为0。现在令F[k][n]表示一个长度为的一个HMBB矩阵的k次方的每个位置和。现在让你求。

乍一看毫无下手之地。但是不要自闭，静下心来，发现这个 p为质数，不正好就是Lucas吗？然后你看看题目标题，不正是sacul，Lucas反过来吗 :)

根据Lucas定理，要使得一个位置为1，当且仅当的 i 在p进制下的每一位都要大于等于 j 在p进制下的每一位。那么，我们再来考虑这个矩阵的k次方怎么做。由于HMBB矩阵是一个01的邻接矩阵，所以可以有图论的性质，它的平方表示从起点出发经过1个点到达终点，n次方则是经过n-1个点。所以这个k次方后矩阵每个位置求和，相当于是从任意点出发，经过k-1个点能够到达任意点的方案数。

那么问题就变成了一个组合计数的问题。我们仔细观察可以发现，如果存在长度为k的一条路，那么表示任意相邻两点可以到达，也就是说任意相邻两点对应HMBB矩阵的位置是1，当然前后关系不能换。所以说，任意相邻两个数字，满足一个数字在p进制下的每一位都小于另一个数字在p进制下的每一位。这样，我们就可以套用经典的组合数学问题，显然每一位是独立的，假设有i位，那么答案就是。于是，可以发现，后面哪个是一个等比数列，利用等比数列公式求和即可。具体见代码：

#include <bits/stdc++.h>
#define mod 1000000007
#define LL long long
#define N 15000010
using namespace std;
#define IO ios::sync_with_stdio(0);cin.tie(0);

bool isp[N];
int p[N];

int qpow(int x,int n)
{
    int res=1;
    while(n)
    {
        if (n&1) res=(LL)res*x%mod;
        x=(LL)x*x%mod; n>>=1;
    }
    return res;
}

void init()
{
    int sz=0;
    for(int i=2;i<N;i++)
    {
        if(!isp[i])p[++sz]=i;
        for(int j=1;j<=sz&&(LL)i*p[j]<N;j++)
        {
            isp[i*p[j]]=1;
            if (i%p[j]==0) break;
        }
    }
    assert(sz>=100000);
}

int sum(int q,int n)
{
    if (q==1) return n;
    return (LL)(qpow(q,n+1)+mod-q)*qpow(q-1,mod-2)%mod;
}

int main()
{
    IO; init();
    int T;
    cin>>T;
    while(T--)
    {
        int ans=0;
        int c,n,k;
        cin>>c>>n>>k;
        int p=::p[c],q=p;
        for(int i=1;i<=k;i++)
        {
            q=(LL)q*(i+p)%mod*qpow(i+1,mod-2)%mod;
            ans=(ans+sum(q,n))%mod;
        }
        cout<<ans<<endl;
    }
    return 0;
}