結果

問題 No.93 ペガサス
ユーザー uwi
提出日時 2014-12-10 14:16:10
言語 Java
(openjdk 23)
結果
AC  
実行時間 320 ms / 5,000 ms
コード長 3,643 bytes
コンパイル時間 2,145 ms
コンパイル使用メモリ 80,272 KB
実行使用メモリ 55,440 KB
最終ジャッジ日時 2024-06-11 20:31:46
合計ジャッジ時間 5,660 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.ByteArrayInputStream;
import java.io.IOException;
import java.io.InputStream;
import java.io.PrintWriter;
import java.util.Arrays;
import java.util.InputMismatchException;

public class Main {
	static InputStream is;
	static PrintWriter out;
	static String INPUT = "";
	
	// 2 2 8 28 152 952 7208
	static void solve()
	{
		out.println(go(ni()));
	}
	
	static long go(int n){
		if(n == 1)return 1;
		int mod = 1000000007;
		long[][][] pre = new long[n+1][2][2]; // [violateしている列の個数][2個上の行がviolateしているかどうか][1個上の行がviolateしているかどうか]
		long[][][] cur = new long[n+1][2][2];
		pre[0][0][0] = 2;
		for(int i = 2;i < n;i++){
			for(int j = 0;j <= i;j++){
				for(int k = 0;k < 2;k++){
					for(int l = 0;l < 2;l++){
						cur[j][k][l] = 0;
					}
				}
			}
			
			for(int j = 0;j <= i;j++){
				for(int k = 0;k < 2;k++){
					// 3個上以上:切断しない, 2-4:切断しない, 1-3:切断しない, 0-2:なにもない
					if(i+1-j-(2-k) >= 0){
						cur[j][0][0] += pre[j][k][0] * (i+1-j-(2-k));
						cur[j][0][0] %= mod;
						cur[j][1][0] += pre[j][k][1] * (i+1-j-(2-k));
						cur[j][1][0] %= mod;
					}
					
					// 3個上以上:切断する, 2-4:切断しない, 1-3:切断しない, 0-2:なにもない
					if(j-1 >= 0){
						cur[j-1][0][0] += pre[j][k][0] * (j-k-0);
						cur[j-1][0][0] %= mod;
						cur[j-1][1][0] += pre[j][k][1] * (j-k-1);
						cur[j-1][1][0] %= mod;
					}
					
					// 3個上以上:切断しない, 2-4:切断する, 1-3:切断しない, 0-2:結合する
					if(k == 1 && j >= 1){
						cur[j][0][1] += pre[j][k][0];
						if(cur[j][0][1] >= mod)cur[j][0][1] -= mod;
						cur[j][1][1] += pre[j][k][1];
						if(cur[j][1][1] >= mod)cur[j][1][1] -= mod;
					}
					
					// 3個上以上:切断しない, 2-4:切断しない, 1-3:切断しない, 0-2:結合する
					if(j+1 <= n){
						cur[j+1][0][1] += pre[j][k][0] * (2-k);
						cur[j+1][0][1] %= mod;
						cur[j+1][1][1] += pre[j][k][1] * (2-k);
						cur[j+1][1][1] %= mod;
					}
					
					// 3個上以上:切断しない, 2-4:切断しない, 1-3:切断する, 0-2:なにもない
					if(j-1 >= 0){
						cur[j-1][0][0] += pre[j][k][1];
						if(cur[j-1][0][0] >= mod)cur[j-1][0][0] -= mod;
					}
				}
			}
			long[][][] dum = pre; pre = cur; cur = dum;
		}
		return pre[0][0][0];
	}
	
	public static void main(String[] args) throws Exception
	{
		long S = System.currentTimeMillis();
		is = INPUT.isEmpty() ? System.in : new ByteArrayInputStream(INPUT.getBytes());
		out = new PrintWriter(System.out);
		
		solve();
		out.flush();
		long G = System.currentTimeMillis();
		tr(G-S+"ms");
	}
	
	private static byte[] inbuf = new byte[1024];
	static int lenbuf = 0, ptrbuf = 0;
	
	private static int readByte()
	{
		if(lenbuf == -1)throw new InputMismatchException();
		if(ptrbuf >= lenbuf){
			ptrbuf = 0;
			try { lenbuf = is.read(inbuf); } catch (IOException e) { throw new InputMismatchException(); }
			if(lenbuf <= 0)return -1;
		}
		return inbuf[ptrbuf++];
	}
	
	private static int ni()
	{
		int num = 0, b;
		boolean minus = false;
		while((b = readByte()) != -1 && !((b >= '0' && b <= '9') || b == '-'));
		if(b == '-'){
			minus = true;
			b = readByte();
		}
		
		while(true){
			if(b >= '0' && b <= '9'){
				num = num * 10 + (b - '0');
			}else{
				return minus ? -num : num;
			}
			b = readByte();
		}
	}
	
	private static void tr(Object... o) { if(INPUT.length() != 0)System.out.println(Arrays.deepToString(o)); }
}
0