結果
| 問題 | No.444 旨味の相乗効果 | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2016-11-11 23:53:44 | 
| 言語 | Java (openjdk 23) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 953 ms / 2,500 ms | 
| コード長 | 1,579 bytes | 
| コンパイル時間 | 4,134 ms | 
| コンパイル使用メモリ | 78,284 KB | 
| 実行使用メモリ | 46,964 KB | 
| 最終ジャッジ日時 | 2024-11-25 10:07:23 | 
| 合計ジャッジ時間 | 17,777 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 5 | 
| other | AC * 23 | 
ソースコード
import java.io.PrintStream;
import java.util.Scanner;
public class Y444 {
	Scanner in = new Scanner(System.in);
	PrintStream out = new PrintStream(System.out);
	long MOD = 1000000007;
	
	
	Y444() throws Exception {
		int n = in.nextInt();
		long c = in.nextLong();
		int[] a = new int[n];
		for (int i = 0; i < n; i++) {
			a[i] = in.nextInt();
		}
		
		Matrix m = new Matrix();
		for (int i = 0; i < n; i++) {
			m.a[i][i] = a[i];
			m.a[i+n][i+n] = a[i];
			for (int j = 0; j < i; j++) {
				m.a[i+n][j] = a[i];
				m.a[i+n][j+n] = a[i];
			}
		}
				
		Matrix power = m.pow(c - 1);
		long ans = 0;
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				ans = (ans + a[j] * power.a[i+n][j]) % MOD;
			}
		}
		power=m.mul(m);
		out.println(ans);
	}
	public static void main(String argv[]) throws Exception {
		new Y444();
	}
	
	
	class Matrix {
		long[][] a = new long[160][160];
		Matrix mul(Matrix that) {
			Matrix ret = new Matrix();
			for (int i = 0; i < a.length; i++) {
				for (int j = 0; j < a.length; j++) {
					for (int k = 0; k < a.length; k++) {
						ret.a[i][k] = (ret.a[i][k] + a[i][j] * that.a[j][k]);
					}
					
					if ((j & 7) == 7) {
						for (int k = 0; k < a.length; k++) {
							ret.a[i][k] %= MOD;
						}
					}
				}
			}
			return ret;
		}
		Matrix pow(long p) {
			Matrix ret = new Matrix();
			Matrix x = this;
			
			for (int i = 0; i < a.length; i++) {
				ret.a[i][i] = 1;
			}
			
			while (p > 0) {
				if ((p & 1) == 1) {
					ret = x.mul(ret);
				}
				x = x.mul(x);
				p >>= 1;
			}
			
			return ret;
		}
	}
}
            
            
            
        