vijos1049——送给圣诞夜的礼品

vijos1049
这题。。我已经不想说什么了

(╯—﹏—)╯（┷━━━┷
这题根据题目意思，暴力做法应该是模拟，模拟k次
然而100%数据k=maxlongint-1;
模拟显然炸飞

好吧我一开始真的不知道这题要怎么搞到log或者sqrt
看了下题解
矩阵快速幂=_=
突然无言以对
其实思路很简单
例如n=4
我们构造一个初始矩阵
{1}
{2}
{3}
{4}
{..}
然后比如我们要一次变化：3214
那么我们再构造一个只有0,1的n*n的矩阵
{0,0,1,0}
{0,1,0,0}
{1,0,0,0}
{0,0,0,1}
就是v[i][s[i]]=1
至于为什么这样呢？
模拟一下就可以理解了
毕竟我这么傻的人都想到了
我们把m次变化看作一整个变化c(a)
这是很容易理解的，如果不能理解。。拿两个变化的矩阵乘一下就知道了
因为k远大于m
所以我们可以使用快速幂来很快的求出c(a)^k/m
余下小于m次的变化我们再去for循环乘一遍就好了

这里注意一下！
本题快速幂一定要非递归！！
递归快速幂re…
【请吸取前人血的教训 正确率已一蹶不振

#include<iostream>#include<cstdio>#include<cmath>#include<algorithm>#include<queue>#include<string>#include<cstring>#define inf 1e9#define ll long long#define For(i,j,k) for(int i=j;i<=k;i++)#define Dow(i,j,k) for(int i=k;i>=j;i--)using namespace std;int n,m,k,tt;struct mat{    int v[101][101];    int size_x,size_y;    mat(){memset(v,0,sizeof v);size_x=size_y=0;}}tmp[201];inline mat mul(mat x,mat y){    mat rul;    rul.size_y=y.size_y;rul.size_x=x.size_x;    For(i,1,x.size_x)        For(j,1,y.size_y)        {            int tmp=0;            For(k,1,x.size_y)                tmp=tmp+x.v[i][k]*y.v[k][j];            rul.v[i][j]=tmp;        }    return rul;}inline void out(mat x){    For(i,1,x.size_x)    {        For(j,1,x.size_y)            cout<<x.v[i][j]<<' ';    }    cout<<endl;}inline mat ksm(mat x,int y){    mat sum=x;y--;for(;y;y>>=1){if(y&1)sum=mul(sum,x);x=mul(x,x);}return sum;}int main(){    scanf("%d%d%d",&n,&m,&k);    mat t;    For(i,1,n)      {        scanf("%d",&tt);        t.v[i][tt]=1;    }    mat b;    For(i,1,n)        b.v[i][1]=i;    b.size_x=n;b.size_y=1;    t.size_x=t.size_y=n;    tmp[1]=t;    For(i,2,m)    {        For(j,1,n)        {            scanf("%d",&tt);            tmp[i].v[j][tt]=1;        }        tmp[i].size_x=tmp[i].size_y=n;        t=mul(tmp[i],t);    }    if(k>=m)        b=mul(ksm(t,k/m),b);    For(i,1,k%m)        b=mul(tmp[i],b);    out(b);}