Description
忙碌了一个学期的 Q老师 决定奖励自己 N 天假期。
假期中不同的穿衣方式会有不同的快乐值。
已知 Q老师 一共有 M 件衬衫,且如果昨天穿的是衬衫 A,今天穿的是衬衫 B,则 Q老师 今天可以获得 f[A][B] 快乐值。
在 N 天假期结束后,Q老师 最多可以获得多少快乐值?
Input
输入文件包含多组测试样例,每组测试样例格式描述如下:
第一行给出两个整数 N M,分别代表假期长度与 Q老师 的衬衫总数。(2 ≤ N ≤ 100000, 1 ≤ M ≤ 100)
接下来 M 行,每行给出 M 个整数,其中第 i 行的第 j 个整数,表示 f[i][j]。(1 ≤ f[i][j] ≤ 1000000)
测试样例组数不会超过 10。
Output
每组测试样例输出一行,表示 Q老师 可以获得的最大快乐值。
Sample Input
3 2
0 1
1 0
4 3
1 2 3
1 2 3
1 2 3
Sample Output
2
9
#include<iostream>
#include<string.h>
using namespace std;
long long L[102][102];
int m=10007;
const int N = 102;
struct Matrix {
long long x[N][N];
Matrix operator*(const Matrix& t) const {
Matrix rs;
for (int i = 0; i < N; ++i) {
for (int j = 0; j < N; ++j) {
rs.x[i][j] = 0;
for (int k = 0; k < N; ++k) {
rs.x[i][j] =max(rs.x[i][j],x[i][k]+t.x[k][j] );
}
}
}
return rs;
}
Matrix() {
memset(x, 0, sizeof(x));
}
Matrix(const Matrix& t) {
memcpy(x, t.x, sizeof(x));
}
};
Matrix qpow(Matrix m, int n) {
Matrix ret;
for (int i = 0; i < 102; i++) {
ret.x[i][i] = 1;
}
while (n) {
if (n & 1) {
ret = m * ret;
}
m = m * m;
n >>= 1;
}
return ret;
}
int main() {
int m,n;
while (scanf("%d%d", &n, &m) != EOF) {
Matrix m1;
for(int i=1;i<=m;i++){
for(int j=1;j<=m;j++){
cin>>m1.x[i][j];
}
}
Matrix mm = qpow(m1, n-1);
long long rs = 0;
for(int i=0;i<=m;i++){
for(int j=0;j<=m;j++){
rs=max(rs,mm.x[i][j]);
}
}
cout << rs-1 << endl;
}
return 0;
}