多校第九场

1001：hdu4686 矩阵快速幂

这题一开始没有想到怎么构造矩阵，纠结了半天，最后ws搞出来了，思路还是蛮简单的。直接把ai*bi的递推式化简，然后看有哪些项，然后构造这些项，再多加一个求和的项。

还有注意一点就是，矩阵不能开太大，8*8的TLE了，5*5就行了。

代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<string>
#include<queue>
#include<cstring>
#define maxn 10005
#define INF 0xfffffff
#define mem(a,b) memset(a,b,sizeof(a))
#define FOR(i,s,t) for(int i=s;i<=t;i++)
#define ull unsigned long long
#define ll __int64
using namespace std;
ll Mod =1000000007;
ll a0,b0,ax,ay,bx,by,n;
struct Mat
{
    int n,m;
    ll a[6][6];
};
Mat A0,A,B,C;
Mat MUL(Mat x,Mat y,ll mod)
{
    Mat tp;
    mem(tp.a,0);
    int n=x.n,m=y.m,l=x.m;
    tp.n=n,tp.m=m;
    for(int i=0; i<n; i++)
    {
        for(int j=0; j<m; j++)
        {
            tp.a[i][j]=0;
            for(int k=0; k<l; k++)
            {
                tp.a[i][j]=(tp.a[i][j]+x.a[i][k]*y.a[k][j]%mod)%mod;
            }
        }
    }
    return tp;
}
Mat POW(Mat tp,Mat x,ll n,ll mod)
{
    while(n)
    {
        if(n&1)
        {
            tp=MUL(tp,x,mod);
        }
        n/=2;
        x=MUL(x,x,mod);
    }
    return tp;
}
void init()
{
    A0.m=A0.n=5;
    mem(A0.a,0);
    A0.a[0][0]=a0*b0%Mod,A0.a[0][1]=a0%Mod,A0.a[0][2]=b0%Mod,A0.a[0][3]=1,A0.a[0][4]=0;

    A.m=A.n=5;
    mem(A.a,0);
    A.a[0][0]=ax*bx%Mod,A.a[1][0]=ax*by%Mod,A.a[2][0]=ay*bx%Mod,A.a[3][0]=ay*by%Mod;
    A.a[1][1]=ax%Mod,A.a[3][1]=ay%Mod,A.a[2][2]=bx%Mod,A.a[3][2]=by%Mod,A.a[3][3]=1;
    A.a[0][4]=1,A.a[4][4]=1;
}
void solve()
{
    C=POW(A0,A,n,Mod);//A0*A^n
    ll ans=C.a[0][4];
    printf("%I64d\n",ans);
}
int main()
{
    while(scanf("%I64d",&n)!=EOF)
    {
        scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&a0,&ax,&ay,&b0,&bx,&by);
        if(!n)
        {
            printf("0\n");
            continue;
        }
        init();
        solve();
    }
    return 0;
}

1005：将输入的两个字符对应一个16进制数，每个16进制数对应另一个16进制数，输入输出即可。