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

提出ソース

結果

問題	No.283 スライドパズルと魔方陣
ユーザー	uwi
提出日時	2015-09-26 17:07:37
言語	Java21 (openjdk 21)
結果	WA
実行時間	-
コード長	9,855 bytes
コンパイル時間	3,764 ms
コンパイル使用メモリ	79,824 KB
実行使用メモリ	51,368 KB
最終ジャッジ日時	2023-08-27 01:57:40
合計ジャッジ時間	11,850 ms
ジャッジサーバーID （参考情報）	judge13 / judge14

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	39 ms 49,192 KB
testcase_01	AC	41 ms 49,332 KB
testcase_02	AC	40 ms 49,312 KB
testcase_03	WA	-
testcase_04	AC	43 ms 49,120 KB
testcase_05	AC	41 ms 49,284 KB
testcase_06	WA	-
testcase_07	WA	-
testcase_08	AC	43 ms 49,404 KB
testcase_09	AC	41 ms 49,440 KB
testcase_10	WA	-
testcase_11	AC	41 ms 49,336 KB
testcase_12	WA	-
testcase_13	AC	41 ms 49,144 KB
testcase_14	WA	-
testcase_15	AC	45 ms 50,172 KB
testcase_16	AC	44 ms 50,144 KB
testcase_17	AC	44 ms 49,152 KB
testcase_18	WA	-
testcase_19	AC	45 ms 49,608 KB
testcase_20	WA	-
testcase_21	AC	53 ms 50,328 KB
testcase_22	AC	52 ms 48,184 KB
testcase_23	AC	52 ms 50,476 KB
testcase_24	AC	52 ms 49,932 KB
testcase_25	AC	53 ms 49,908 KB
testcase_26	WA	-
testcase_27	WA	-
testcase_28	WA	-
testcase_29	WA	-
testcase_30	AC	60 ms 49,744 KB
testcase_31	WA	-
testcase_32	AC	62 ms 50,836 KB
testcase_33	WA	-
testcase_34	AC	59 ms 50,320 KB
testcase_35	WA	-
testcase_36	AC	59 ms 50,544 KB
testcase_37	WA	-
testcase_38	AC	61 ms 49,352 KB

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

ソースコード

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_2 {
	InputStream is;
	PrintWriter out;
	String INPUT = "";
	
	int[][] KP(int[][] a, int[][] b)
	{
		int n = a.length, m = b.length;
		int[][] c = new int[n*m][n*m];
		for(int i = 0;i < n*m;i++){
			for(int j = 0;j < n*m;j++){
				c[i][j] = (a[i/m][j/m]-1)*m*m + b[i%m][j%m];
			}
		}
		return c;
	}
	
	void solve()
	{
//		{
////			int[][] tw = KP(constructX4(4), what(11, 1));
//			int[][] tw = what2(50, 1);
//			tr(parity(tw));
//			tr(ismagic(tw));
//		}
			
		int n = ni();
		int[] a = na(n*n);
		for(int i = 0;i < n*n;i++){
			if(a[i] == 0)a[i] = n*n;
		}
		int parity = parity(a);
		
		if(n == 1){
			out.println("possible");
			out.println(1);
			return;
		}
		if(n == 2){
			out.println("impossible");
			return;
		}
		if((n == 4 || n % 8 == 0) && parity == 1){
			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 if(n % 8 == 4 && parity == 1){
			ms = KP(constructX4(4), what(n/4, 1));
		}else{
			ms = constructMagicSquare(n);
		}
		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();
			if(parity(a) == want)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];
			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;
//			if(!ismagic(a))throw new RuntimeException();
			if(parity(a) == want)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[][] 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_2().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

結果

テストケース

ソースコード