codeforce 583D 矩阵快速幂

/**
codeforce 583D Once Again...
给定一个数列包含n个数，将它复制t次，得到一个新的数列，
求这个数列的最长非递减子序列，简称最长上升子序列（非严格）。

新技能GET！- 利用传递性构造矩阵，用快速幂解决。
构造一个矩阵A[i][j],每一项表示
从0 到 j中， 起始数大于等于a[i],以a[j]结尾, 的最长上升子序列的长度。
可以发现矩阵两项之间具有传递性，
A[i][j] = max(A[i][j], x[i][k] + y[k][j]).
那么直接求出初始矩阵然后t次幂之后的到最终状态矩阵，
取最终矩阵的最大值就是答案。
**/
#include <bits/stdc++.h>
using namespace std;

const int inf = 1<<28;
/******************
矩阵快速幂
*******************/
const int N = 110;
typedef long long LL;
struct Matrix{
    int n,m;
    int a[N][N];

    Matrix(){
        memset(a,0,sizeof(a));
    }

    Matrix(int x,int y):n(x),m(y){
        memset(a,0,sizeof(a));
    }

    void InitOne(int x,int y){
        n =x, m = y;
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++) if(i==j) a[i][j] = 1;
        }
    }
};

Matrix multi(Matrix a,Matrix b){
    Matrix ret(a.n,b.m);
    for(int i=0;i<a.n;i++){
        for(int j=0;j<b.m;j++){
            ret.a[i][j] = -inf;
            for(int k=0;k<b.n;k++){
                ret.a[i][j] = max(ret.a[i][j],a.a[i][k] + b.a[k][j]);
            }
        }
    }
    return ret;
}

Matrix pow_mod(Matrix x,int n,int s){
    Matrix res(s,s);
    //res.InitOne(s,s);
    while(n){
        if(n&1) res = multi(res,x);
        n >>= 1;
        x = multi(x,x);
    }
    return res;
}
//
int a[110];
int n,t;

int main(){
    cin >> n >> t;
    for(int i=0;i<n;i++) cin >> a[i];
    Matrix ans;
    ans.n = ans.m = n;
    for(int s=0;s<n;s++){
        for(int i=0;i<n;i++){
            if(a[i] < a[s]) {
                ans.a[s][i] = -inf;
                continue;
            }
            ans.a[s][i] = 1;
            for(int j=0;j<i;j++){
                if(a[j] <= a[i]){
                    ans.a[s][i] = max(ans.a[s][i], ans.a[s][j] + 1);
                }
            }
        }
    }
    ans = pow_mod(ans,t,n);
    int ret = 0;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            ret = max(ret,ans.a[i][j]);
        }
    }
    cout <<ret << endl;
    return 0;
}