习题日常第十三练

1递推（题目链接）

这个题是个矩阵快速幂的典型题，主要难在如何找出矩阵递推关系。具体代码如下（矩阵一定要开long long ）。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=1e9+9;
using namespace std;
ll a[4][4],b[4][4],c[4][4]={{2, 0, 0, 0}, {1, 1, 0, 0}, {2, 2, 1, 0}, {1, 1, 1, 1}};
void mypow(ll n)
{
    ll i,j,k,l,p[4][4];
    if(n>0)
    {
    for(i=0;i<4;i++)
        for(j=0;j<4;j++)
        b[i][j]=c[i][j];
    ll d[4][4]={{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}};
    while(n)
    {
        if(n&1)
        {
            memset(p,0,sizeof(p));
            for(i=0;i<4;i++)
                for(j=0;j<4;j++)
                  for(k=0;k<4;k++)
            {
                p[i][j]+=(d[i][k]*b[k][j])%mod;
            }
              for(i=0;i<4;i++)
                for(j=0;j<4;j++)
                  d[i][j]=p[i][j]%mod;
        }
        memset(p,0,sizeof(p));
        for(i=0;i<4;i++)
                for(j=0;j<4;j++)
                  for(k=0;k<4;k++)
            {
                p[i][j]+=(b[i][k]*b[k][j])%mod;
            }
          for(i=0;i<4;i++)
                for(j=0;j<4;j++)
                  b[i][j]=p[i][j]%mod;
            n>>=1;
    }
    memcpy(a,d,sizeof(d));
    }
}
int main()
{
    ios::sync_with_stdio(false);
    ll i,j,n,m;
    while(cin>>n)
    {
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        mypow(n);
        cout<<a[3][0]%mod<<endl;
    }
  return 0;
}

2哲哲的疑惑（题目链接）

这个题是一个组合数学的题，对于思维的转换还是有点挑战的。C(k,m)指的是从k个里面挑m个，但是这k个又是从n个里面挑出来的。那岂不是可以直接从n个里面挑m个，再将n-m个分配给l个球，此时n-m不需要全部用完，所以分配种数有(n-m)的l次方种。具体计算的时候需要用到快速幂和费马小定理。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=998244353;
using namespace std;
ll mypow(ll n,ll a)
{
    ll ans=1;
    while(n)
    {
        if(n&1)
            ans=(ans*a)%mod;
             a=(a*a)%mod;
             n>>=1;
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    ll i,j,n,m,l,mp,lp,cp;
    cin>>n>>m>>l;
    mp=cp=lp=1;
    for(i=1;i<=n;i++)
    {
        lp*=i;
        lp%=mod;
    }
    for(i=1;i<=n-m;i++)
    {
        mp*=i;
        mp%=mod;
    }
    for(i=1;i<=m;i++)
    {
        cp*=i;
        cp%=mod;
    }
    if(n<m)
    {
        cout<<0;
    }
    else
    {
        ll ans=(((lp*mypow(mod-2,cp))%mod)*mypow(mod-2,mp))%mod;
        ans=(ans*mypow(l,n-m))%mod;
        cout<<ans;
    }
  return 0;
}

3洋灰三角（题目链接）

这个题也是一个关于矩阵快速幂的题。这个题给出的是数列第n项的求法，求前n项和。这些题主要难在对于矩阵该有量的确定以及求关系矩阵。接下来的就可以用矩阵快速幂解决了。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=1e9+7;
using namespace std;
ll a[3][3],b[3][3],c[3][3],d[3][3];
void mypow(ll n)
{
    int i,j,k;
    memset(b,0,sizeof(b));
    for(i=0;i<3;i++)
        b[i][i]=1;
    while(n)
    {
        if(n&1)
        {
            memset(d,0,sizeof(d));
            for(i=0;i<3;i++)
                for(j=0;j<3;j++)
                    for(k=0;k<3;k++)
            {
                d[i][j]+=(b[i][k]*a[k][j])%mod;
            }
            for(i=0;i<3;i++)
                for(j=0;j<3;j++)
                b[i][j]=d[i][j]%mod;
        }
        memset(d,0,sizeof(d));
           for(i=0;i<3;i++)
                for(j=0;j<3;j++)
                    for(k=0;k<3;k++)
            {
                d[i][j]+=(a[i][k]*a[k][j])%mod;
            }
             for(i=0;i<3;i++)
                for(j=0;j<3;j++)
                a[i][j]=d[i][j]%mod;
                n>>=1;
    }
}
int main()
{
    ios::sync_with_stdio(false);
    ll i,j,n,m,l,p,k;
    cin>>n>>k>>p;
    memset(a,0,sizeof(a));
    a[0][1]=a[0][0]=k;
    a[1][1]=a[2][0]=a[2][1]=a[2][2]=1;
    mypow(n-1);
    ll ans=(b[1][1]+b[0][1]+b[2][1]*p)%mod;
    cout<<ans;
  return 0;
}