#51195 (Java21) No.283 スライドパズルと魔方陣

提出ソース

結果

問題	No.283 スライドパズルと魔方陣
ユーザー	uwi
提出日時	2015-09-26 17:52:49
言語	Java21 (openjdk 21)
結果	WA
実行時間	-
コード長	10,700 bytes
コンパイル時間	4,593 ms
コンパイル使用メモリ	81,132 KB
実行使用メモリ	50,740 KB
最終ジャッジ日時	2023-08-27 01:58:28
合計ジャッジ時間	11,600 ms
ジャッジサーバーID （参考情報）	judge12 / judge14

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	42 ms 49,272 KB
testcase_01	AC	41 ms 47,288 KB
testcase_02	AC	39 ms 49,276 KB
testcase_03	AC	40 ms 49,368 KB
testcase_04	WA	-
testcase_05	WA	-
testcase_06	WA	-
testcase_07	WA	-
testcase_08	AC	41 ms 49,296 KB
testcase_09	AC	41 ms 49,280 KB
testcase_10	WA	-
testcase_11	AC	45 ms 49,408 KB
testcase_12	AC	42 ms 49,368 KB
testcase_13	WA	-
testcase_14	WA	-
testcase_15	AC	44 ms 50,004 KB
testcase_16	WA	-
testcase_17	AC	45 ms 49,328 KB
testcase_18	AC	43 ms 49,356 KB
testcase_19	AC	44 ms 49,388 KB
testcase_20	WA	-
testcase_21	WA	-
testcase_22	WA	-
testcase_23	AC	53 ms 49,468 KB
testcase_24	WA	-
testcase_25	AC	53 ms 50,420 KB
testcase_26	WA	-
testcase_27	AC	71 ms 50,320 KB
testcase_28	WA	-
testcase_29	WA	-
testcase_30	WA	-
testcase_31	AC	60 ms 50,484 KB
testcase_32	AC	63 ms 50,536 KB
testcase_33	AC	63 ms 50,388 KB
testcase_34	WA	-
testcase_35	AC	58 ms 49,508 KB
testcase_36	WA	-
testcase_37	AC	60 ms 50,632 KB
testcase_38	WA	-

権限があれば一括ダウンロードができます

ソースコード

raw source code

package q6xx;

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

public class Q664_3 {
	InputStream is;
	PrintWriter out;
	String INPUT = "";
	
	void solve()
	{
//		for(int i = 1;i <= 12;i++){
//			tr(i);
//			int[][] tw = panmagicX4(i*4, 1);
////			tr(parity(tw));
//			tr(ismagic(tw));
//		}
//		if(true)return;
			
		int n = ni();
		int[] a = na(n*n);
		int[] f = new int[n*n+1];
		for(int i = 0;i < n*n;i++){
			f[a[i]]++;
			if(a[i] == 0){
				a[i] = n*n;
			}
		}
		for(int i = 0;i < n*n;i++){
			if(f[i] != 1)throw new RuntimeException();
		}
		int parity = parity(a);
		
		if(n == 1){
			out.println("possible");
			out.println(1);
			return;
		}
		if(n == 2){
			out.println("impossible");
			return;
		}
		out.println("possible");
		
		int[][] ms = null;
		if(n % 2 == 1){
			ms = what(n, parity);
		}else if(n % 4 == 2){
			ms = what2(n, parity);
		}else{
			ms = panmagicX4(n, parity);
		}
		if(!ismagic(ms)){
			throw new RuntimeException();
		}
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				if(j > 0)out.print(" ");
				out.print(ms[i][j]);
			}
			out.println();
		}
	}
	
	static int[][] what2(int n, int want)
	{
		n /= 2;
		int[] A = new int[n];
		for(int i = 0;i < n;i++)A[i] = i;
		int[] B = new int[n];
		for(int i = 0;i < n;i++)B[i] = i;
		do{
			int[][] a = new int[n][n];
			for(int i = 0;i < n;i++){
				for(int j = 0;j < n;j++){
					a[i][j] = A[(i+j+n-n/2)%n]*n+B[(i+n-j+n/2)%n]+1;
	//				a[i][j] = A[(i+j*0)%n];
	//				a[i][j] = B[(i+n-j+n/2)%n];
	//				a[i][j] = (i*(n-s)+j+2)%n;
				}
			}
			if(!ismagic(a))continue;
			a = constructByLUXMethod(a);
			if(!ismagic(a))throw new RuntimeException();
			int lp = 0;
			for(int i = 0;i < a.length;i++){
				for(int j = 0;j < a.length;j++){
					if(a[i][j] == 4*n*n){
						lp = (i^j)&1;
					}
				}
			}
			if(parity(a) == (want^lp))return a;
		}while(nextPermutation(B));
		return null;
	}
	
	static int[][] what(int n, int want)
	{
		int[] A = new int[n];
		for(int i = 0;i < n;i++)A[i] = i;
		int[] B = new int[n];
		for(int i = 0;i < n;i++)B[i] = i;
		do{
			int[][] a = new int[n][n];
			int lp = 0;
			for(int i = 0;i < n;i++){
				for(int j = 0;j < n;j++){
					a[i][j] = A[(i+j+n-n/2)%n]*n+B[(i+n-j+n/2)%n]+1;
					if(a[i][j] == n*n){
						if(a[i][j] == n*n)lp = (i^j)&1;
					}
	//				a[i][j] = A[(i+j*0)%n];
	//				a[i][j] = B[(i+n-j+n/2)%n];
	//				a[i][j] = (i*(n-s)+j+2)%n;
				}
			}
			if(!ismagic(a))continue;
//			if(!ismagic(a))throw new RuntimeException();
			if(parity(a) == (want^lp))return a;
		}while(nextPermutation(B));
		return null;
	}
	
	public static boolean nextPermutation(int[] a) {
		int n = a.length;
		int i;
		for (i = n - 2; i >= 0 && a[i] >= a[i + 1]; i--)
			;
		if (i == -1)
			return false;
		int j;
		for (j = i + 1; j < n && a[i] < a[j]; j++)
			;
		int d = a[i];
		a[i] = a[j - 1];
		a[j - 1] = d;
		for (int p = i + 1, q = n - 1; p < q; p++, q--) {
			d = a[p];
			a[p] = a[q];
			a[q] = d;
		}
		return true;
	}
	
	static boolean ismagic(int[][] a)
	{
		int n = a.length;
		BitSet bs = new BitSet();
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				if(a[i][j] <= 0 || a[i][j] > n*n)return false;
				bs.set(a[i][j]);
			}
		}
		if(bs.cardinality() != n*n)return false;
		
		// 1-base
		int S = n*(n*n+1)/2;
		for(int i = 0;i < n;i++){
			int ls = 0;
			for(int j = 0;j < n;j++)ls += a[i][j];
			if(ls != S)return false;
		}
		for(int i = 0;i < n;i++){
			int ls = 0;
			for(int j = 0;j < n;j++)ls += a[j][i];
			if(ls != S)return false;
		}
		{
			int ls = 0;
			for(int i = 0;i < n;i++)ls += a[i][i];
			if(ls != S)return false;
		}
		{
			int ls = 0;
			for(int i = 0;i < n;i++)ls += a[n-1-i][i];
			if(ls != S)return false;
		}
		return true;
	}
	
	static int parity(int[][] t)
	{
		int n = t.length;
		int[] u = new int[n*n];
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				u[i*n+j] = t[i][j];
			}
		}
		return parity(u);
	}
	
	public static int[][] constructBySiameseMethod2(int n)
	{
		if(n <= 0 || n % 2 == 0)throw new RuntimeException();
		int[][] ret = new int[n][n];
		int r = n/2-1, c = n/2;
		for(int i = 1;i <= n*n;i++){
			ret[r][c] = i;
			int nr = r-1, nc = c+1;
			if(nr < 0)nr += n;
			if(nc >= n)nc -= n;
			if(ret[nr][nc] != 0){
				nr = r-2; nc = c;
				if(nr < 0)nr += n;
			}
			r = nr; c = nc;
		}
		return ret;
	}
	
	public static int parity(int[] a)
	{
		int x = 0;
		int n = a.length;
		for(int i = 0;i < n;i++){
			if(i != a[i]-1){
				int j = a[i]-1;
				int d = a[i]; a[i] = a[j]; a[j] = d;
				x ^= 1;
			}
		}
		return x;
	}
	
	public static int[][] constructMagicSquare(int n)
	{
		if(n < 0)throw new RuntimeException();
		if(n % 4 == 0){
			// 4m
			return constructX4(n);
		}else if(n % 2 == 1){
			// 2m+1
			return constructBySiameseMethod(n);
		}else{
			// 4m+2
			return constructByLUXMethod(constructBySiameseMethod(n/2));
		}
	}
	
	public static int[][] panmagicX4(int n, int want)
	{
		int[][] ret = new int[n][n];
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				int A = i%2 == 0 ? 
						(j < n/2 ? 1+j : n-(j-n/2)) :
						(j < n/2 ? n-j : 1+(j-n/2));
				int B = j%2 == 0 ? 
						(i < n/2 ? n/2+1+i : n/2-(i-n/2)) :
						(i < n/2 ? n/2-i : n/2+1+(i-n/2));
				ret[i][j] = A + n * B - n;
			}
		}
		int par = parity(ret);
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				if(ret[i][j] == n*n){
					par ^= (i^j)&1;
				}
			}
		}
		if(par != want){
			for(int j = 1;j < n;j++){
				int[][] nret = new int[n][];
				for(int i = 0;i < n;i++){
					nret[(i+1)%n] = ret[i];
				}
				if(!ismagic(nret))throw new RuntimeException();
				if(par == (parity(nret)^(j&1))){
					ret = nret;
					return ret;
				}
			}
			throw new RuntimeException();
		}
		return ret;
	}
	
	public static int[][] constructX4(int n)
	{
		if(n < 0 || n % 4 != 0)throw new RuntimeException();
		int[][] ret = new int[n][n];
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				int p = (i^j)&3;
				if(p == 0 || p == 3){
					ret[i][j] = i*n+j+1;
				}else{
					ret[i][j] = n*n-(i*n+j);
				}
			}
		}
		return ret;
	}
	
	public static int[][] constructByLUXMethod(int[][] a)
	{
		int n = a.length;
		if(n <= 0 || n % 2 == 0)throw new RuntimeException();
		int[][] ret = new int[2*n][2*n];
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				int base = a[i][j]-1<<2;
				if(i > n/2+1){
					// X
					ret[i*2][j*2] = base + 1;
					ret[i*2][j*2+1] = base + 4;
					ret[i*2+1][j*2] = base + 3;
					ret[i*2+1][j*2+1] = base + 2;
				}else if(j != n/2 && i <= n/2 || j == n/2 && i <= n/2+1 && i != n/2){
					// L
					ret[i*2][j*2] = base + 4;
					ret[i*2][j*2+1] = base + 1;
					ret[i*2+1][j*2] = base + 2;
					ret[i*2+1][j*2+1] = base + 3;
				}else{
					// U
					ret[i*2][j*2] = base + 1;
					ret[i*2][j*2+1] = base + 4;
					ret[i*2+1][j*2] = base + 2;
					ret[i*2+1][j*2+1] = base + 3;
				}
			}
		}
		return ret;
	}
	
	// val(i,j)=n*((i+j+1+n/2))%n)+(i+2*j+1)%n+1 TODO
	public static int[][] constructBySiameseMethod(int n)
	{
		if(n <= 0 || n % 2 == 0)throw new RuntimeException();
		int[][] ret = new int[n][n];
		int r = 0, c = n/2;
		for(int i = 1;i <= n*n;i++){
			ret[r][c] = i;
			int nr = r-1, nc = c+1;
			if(nr < 0)nr += n;
			if(nc >= n)nc -= n;
			if(ret[nr][nc] != 0){
				nr = r+1; nc = c;
				if(nr >= n)nr -= n;
			}
			r = nr; c = nc;
		}
		return ret;
	}
	
	public static boolean check(int[][] a)
	{
		int n = a.length;
		int[] f = new int[n*n+1];
		for(int i = 0;i < n;i++){
			for(int j = 0;j < n;j++){
				if(a[i][j] < 0 || a[i][j] > n*n)return false;
				if(++f[a[i][j]] > 1)return false;
			}
		}
		int base = n*(n*n+1)/2;
		// row
		for(int i = 0;i < n;i++){
			int s = 0;
			for(int j = 0;j < n;j++)s += a[i][j];
			if(s != base)return false;
		}
		// col
		for(int i = 0;i < n;i++){
			int s = 0;
			for(int j = 0;j < n;j++)s += a[j][i];
			if(s != base)return false;
		}
		// diagonal
		{
			int s = 0;
			for(int j = 0;j < n;j++)s += a[j][j];
			if(s != base)return false;
		}
		{
			int s = 0;
			for(int j = 0;j < n;j++)s += a[n-1-j][j];
			if(s != base)return false;
		}
		return true;
	}
	
	void run() throws Exception
	{
		is = INPUT.isEmpty() ? System.in : new ByteArrayInputStream(INPUT.getBytes());
		out = new PrintWriter(System.out);
		
		long s = System.currentTimeMillis();
		solve();
		out.flush();
		if(!INPUT.isEmpty())tr(System.currentTimeMillis()-s+"ms");
	}
	
	public static void main(String[] args) throws Exception { new Q664_3().run(); }
	
	private byte[] inbuf = new byte[1024];
	private int lenbuf = 0, ptrbuf = 0;
	
	private 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 boolean isSpaceChar(int c) { return !(c >= 33 && c <= 126); }
	private int skip() { int b; while((b = readByte()) != -1 && isSpaceChar(b)); return b; }
	
	private double nd() { return Double.parseDouble(ns()); }
	private char nc() { return (char)skip(); }
	
	private String ns()
	{
		int b = skip();
		StringBuilder sb = new StringBuilder();
		while(!(isSpaceChar(b))){ // when nextLine, (isSpaceChar(b) && b != ' ')
			sb.appendCodePoint(b);
			b = readByte();
		}
		return sb.toString();
	}
	
	private char[] ns(int n)
	{
		char[] buf = new char[n];
		int b = skip(), p = 0;
		while(p < n && !(isSpaceChar(b))){
			buf[p++] = (char)b;
			b = readByte();
		}
		return n == p ? buf : Arrays.copyOf(buf, p);
	}
	
	private char[][] nm(int n, int m)
	{
		char[][] map = new char[n][];
		for(int i = 0;i < n;i++)map[i] = ns(m);
		return map;
	}
	
	private int[] na(int n)
	{
		int[] a = new int[n];
		for(int i = 0;i < n;i++)a[i] = ni();
		return a;
	}
	
	private 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 long nl()
	{
		long num = 0;
		int 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) { System.out.println(Arrays.deepToString(o)); }
}

yukicoder

結果

テストケース

ソースコード