#488984 (Java21) No.1069 電柱 / Pole (Hard)

提出ソース

結果

問題	No.1069 電柱 / Pole (Hard)
ユーザー	uwi
提出日時	2020-05-29 22:39:39
言語	Java21 (openjdk 21)
結果	AC
実行時間	1,402 ms / 2,000 ms
コード長	9,272 bytes
コンパイル時間	4,725 ms
コンパイル使用メモリ	84,500 KB
実行使用メモリ	93,124 KB
最終ジャッジ日時	2023-09-30 21:18:37
合計ジャッジ時間	22,126 ms
ジャッジサーバーID （参考情報）	judge14 / judge13

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	113 ms 55,768 KB
testcase_01	AC	113 ms 55,668 KB
testcase_02	AC	113 ms 55,840 KB
testcase_03	AC	114 ms 56,268 KB
testcase_04	AC	928 ms 77,568 KB
testcase_05	AC	1,402 ms 93,124 KB
testcase_06	AC	195 ms 64,888 KB
testcase_07	AC	197 ms 65,672 KB
testcase_08	AC	197 ms 58,988 KB
testcase_09	AC	185 ms 62,048 KB
testcase_10	AC	191 ms 65,012 KB
testcase_11	AC	117 ms 55,776 KB
testcase_12	AC	121 ms 55,884 KB
testcase_13	AC	180 ms 55,528 KB
testcase_14	AC	191 ms 57,528 KB
testcase_15	AC	183 ms 57,148 KB
testcase_16	AC	165 ms 58,860 KB
testcase_17	AC	131 ms 55,436 KB
testcase_18	AC	165 ms 58,796 KB
testcase_19	AC	181 ms 59,960 KB
testcase_20	AC	131 ms 56,048 KB
testcase_21	AC	196 ms 62,168 KB
testcase_22	AC	133 ms 55,972 KB
testcase_23	AC	270 ms 60,856 KB
testcase_24	AC	201 ms 69,604 KB
testcase_25	AC	200 ms 69,352 KB
testcase_26	AC	194 ms 69,540 KB
testcase_27	AC	196 ms 69,228 KB
testcase_28	AC	197 ms 70,080 KB
testcase_29	AC	122 ms 57,896 KB
testcase_30	AC	118 ms 55,648 KB
testcase_31	AC	116 ms 56,124 KB
testcase_32	AC	119 ms 55,616 KB
testcase_33	AC	117 ms 55,656 KB
testcase_34	AC	115 ms 55,684 KB
testcase_35	AC	115 ms 55,840 KB
testcase_36	AC	118 ms 55,600 KB
testcase_37	AC	131 ms 55,620 KB
testcase_38	AC	135 ms 55,828 KB
testcase_39	AC	119 ms 55,456 KB
testcase_40	AC	124 ms 55,880 KB
testcase_41	AC	139 ms 56,088 KB
testcase_42	AC	131 ms 55,316 KB
testcase_43	AC	144 ms 56,228 KB
testcase_44	AC	193 ms 63,200 KB
testcase_45	AC	166 ms 57,144 KB
testcase_46	AC	192 ms 60,904 KB
testcase_47	AC	198 ms 60,376 KB
testcase_48	AC	206 ms 66,108 KB
testcase_49	AC	198 ms 62,968 KB
testcase_50	AC	122 ms 55,688 KB
testcase_51	AC	204 ms 62,136 KB
testcase_52	AC	120 ms 55,376 KB
testcase_53	AC	197 ms 64,812 KB
testcase_54	AC	115 ms 56,392 KB
testcase_55	AC	116 ms 55,480 KB
testcase_56	AC	114 ms 56,104 KB
testcase_57	AC	118 ms 55,548 KB
testcase_58	AC	114 ms 56,380 KB
testcase_59	AC	120 ms 55,788 KB
testcase_60	AC	115 ms 55,644 KB
testcase_61	AC	118 ms 55,596 KB
testcase_62	AC	114 ms 55,652 KB
testcase_63	AC	112 ms 55,548 KB
testcase_64	AC	120 ms 55,776 KB
testcase_65	AC	113 ms 55,672 KB
testcase_66	AC	115 ms 55,440 KB
testcase_67	AC	115 ms 55,876 KB
testcase_68	AC	115 ms 55,868 KB
testcase_69	AC	198 ms 68,552 KB
testcase_70	AC	138 ms 65,444 KB
testcase_71	AC	149 ms 67,468 KB
testcase_72	AC	124 ms 57,612 KB
testcase_73	AC	194 ms 64,352 KB
testcase_74	AC	213 ms 69,472 KB
testcase_75	AC	202 ms 72,304 KB
testcase_76	AC	193 ms 64,428 KB
testcase_77	AC	194 ms 62,544 KB
testcase_78	AC	121 ms 55,872 KB
testcase_79	AC	200 ms 61,892 KB
testcase_80	AC	199 ms 62,788 KB
testcase_81	AC	183 ms 57,244 KB
testcase_82	AC	179 ms 56,652 KB

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

ソースコード

raw source code

package contest200529;
import java.io.ByteArrayInputStream;
import java.io.IOException;
import java.io.InputStream;
import java.io.PrintWriter;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.HashMap;
import java.util.HashSet;
import java.util.InputMismatchException;
import java.util.List;
import java.util.Map;
import java.util.Set;

public class G2 {
	InputStream is;
	PrintWriter out;
	String INPUT = "";
	
	void solve()
	{
		int n = ni(), m = ni(), K = ni();
		int X = ni()-1, Y = ni()-1;
		int[][] co = new int[n][];
		for(int i = 0;i < n;i++){
			co[i] = na(2);
		}
		int[] from= new int[m];
		int[] to = new int[m];
		int[] ws = new int[m];
		for(int i = 0;i < m;i++){
			from[i] = ni()-1;
			to[i] = ni()-1;
			ws[i] = (int)d2(co[from[i]], co[to[i]]);
		}
		
		int[][][] g = packWU(n, from, to, ws);
		List<Double> res = yen(g, X, Y, K);
		for(int i = 0;i < K;i++){
			if(i < res.size()){
				out.printf("%.14f\n", res.get(i));
			}else{
				out.println(-1);
			}
		}
	}
	
	public static int[][][] packWU(int n, int[] from, int[] to, int[] w) {
		int[][][] g = new int[n][][];
		int[] p = new int[n];
		for (int f : from)
			p[f]++;
		for (int t : to)
			p[t]++;
		for (int i = 0; i < n; i++)
			g[i] = new int[p[i]][2];
		for (int i = 0; i < from.length; i++) {
			--p[from[i]];
			g[from[i]][p[from[i]]][0] = to[i];
			g[from[i]][p[from[i]]][1] = w[i];
			--p[to[i]];
			g[to[i]][p[to[i]]][0] = from[i];
			g[to[i]][p[to[i]]][1] = w[i];
		}
		return g;
	}

	
	public static long d2(int[] a, int[] b) { return (long)(b[0]-a[0])*(b[0]-a[0])+(long)(b[1]-a[1])*(b[1]-a[1]); }

	
	public static class MinHeapD {
		public double[] a;
		public int[] map;
		public int[] imap;
		public int n;
		public int pos;
		public static double INF = Double.POSITIVE_INFINITY;
		
		public MinHeapD(int m)
		{
			n = Integer.highestOneBit((m+1)<<1);
			a = new double[n];
			map = new int[n];
			imap = new int[n];
			Arrays.fill(a, INF);
			Arrays.fill(map, -1);
			Arrays.fill(imap, -1);
			pos = 1;
		}
		
		public double add(int ind, double x)
		{
			int ret = imap[ind];
			if(imap[ind] < 0){
				a[pos] = x; map[pos] = ind; imap[ind] = pos;
				pos++;
				up(pos-1);
			}
			return ret != -1 ? a[ret] : x;
		}
		
		public double update(int ind, double x)
		{
			int ret = imap[ind];
			if(imap[ind] < 0){
				a[pos] = x; map[pos] = ind; imap[ind] = pos;
				pos++;
				up(pos-1);
			}else{
				a[ret] = x;
				up(ret);
				down(ret);
			}
			return x;
		}
		
		public double remove(int ind)
		{
			if(pos == 1)return INF;
			if(imap[ind] == -1)return INF;
			
			pos--;
			int rem = imap[ind];
			double ret = a[rem];
			map[rem] = map[pos];
			imap[map[pos]] = rem;
			imap[ind] = -1;
			a[rem] = a[pos];
			a[pos] = INF;
			map[pos] = -1;
			
			up(rem);
			down(rem);
			return ret;
		}
		
		public double min() { return a[1]; }
		public int argmin() { return map[1]; }
		public int size() {	return pos-1; }
		
		private void up(int cur)
		{
			for(int c = cur, p = c>>>1;p >= 1 && a[p] > a[c];c>>>=1, p>>>=1){
				double d = a[p]; a[p] = a[c]; a[c] = d;
				int e = imap[map[p]]; imap[map[p]] = imap[map[c]]; imap[map[c]] = e;
				e = map[p]; map[p] = map[c]; map[c] = e;
			}
		}
		
		private void down(int cur)
		{
			for(int c = cur;2*c < pos;){
				int b = a[2*c] < a[2*c+1] ? 2*c : 2*c+1;
				if(a[b] < a[c]){
					double d = a[c]; a[c] = a[b]; a[b] = d;
					int e = imap[map[c]]; imap[map[c]] = imap[map[b]]; imap[map[b]] = e;
					e = map[c]; map[c] = map[b]; map[b] = e;
					c = b;
				}else{
					break;
				}
			}
		}
	}
	
	static List<Integer> dijkpath(int[][][] g, int src, int sink, boolean[] ng)
	{
		int n = g.length;
		double I = Double.POSITIVE_INFINITY;
		double[] td = new double[n];
		Arrays.fill(td, I);
		MinHeapD h = new MinHeapD(n);
		h.add(sink, 0);
		td[sink] = 0;
		
		int[] prev = new int[n];
		while(h.size() > 0){
			int cur = h.argmin();
			h.remove(cur);
			
			for(int[] e : g[cur]){
				int nex = e[0];
				if(e[1] < 0)continue;
				if(ng[nex])continue;
				double nd = td[cur] + Math.sqrt(e[1]);
				if(nd < td[nex]){
					td[nex] = nd;
					prev[nex] = cur;
					h.update(nex, nd);
				}
			}
		}
		
		if(Double.isInfinite(td[src]))return null;
		
		List<Integer> ret = new ArrayList<>();
		for(int cur = src;cur != sink;cur = prev[cur]){
			ret.add(cur);
		}
		ret.add(sink);
		return ret;
	}
	
	static long h(List<Integer> a)
	{
		long h = -a.size() * a.size();
		for(int x : a){
			h =h * 1000000009 + x;
		}
		return h;
	}
	
	static double len(List<Integer> path, int[][] link)
	{
		double s = 0;
		for(int j = 0;j < path.size()-1;j++){
			s += Math.sqrt(link[path.get(j)][path.get(j+1)]);
		}
		return s;
	}
	
	static List<Double> yen(int[][][] g, int src, int sink, int K)
	{
		int n = g.length;
		List<List<Integer>> A = new ArrayList<>();
		A.add(dijkpath(g, src, sink, new boolean[n]));
		
		int[][] link = new int[n][n];
		for(int i = 0;i < n;i++){
			for(int[] e : g[i]){
				link[i][e[0]] = e[1];
			}
		}
		
		Map<Long, List<Integer>> B = new HashMap<>();
		Set<Long> done = new HashSet<>();
		
		List<Double> ret = new ArrayList<>();
		done.add(h(A.get(0)));
		ret.add(len(A.get(0), link));
		for(int k = 1;k <= K-1;k++){
			List<Integer> last = A.get(A.size()-1);
			for(int i = 0;i < last.size()-1;i++){
				int spurNode = last.get(i);
				
				inner:
				for(List<Integer> p : A){
					if(p.size() < i+2)continue;
					for(int j = 0;j <= i;j++){
						if(!p.get(j).equals(last.get(j))){
							continue inner;
						}
					}
					int rf = p.get(i);
					int rt = p.get(i+1);
					for(int[] e : g[rf]){
						if(e[0] == rt && e[1] >= 0)e[1] = -e[1]-1;
					}
					for(int[] e : g[rt]){
						if(e[0] == rf && e[1] >= 0)e[1] = -e[1]-1;
					}
				}
				boolean[] remed = new boolean[n];
				for(int j = 0;j < i;j++){
					remed[last.get(j)] = true;
				}
				
				List<Integer> spurPath = dijkpath(g, spurNode, sink, remed);
				if(spurPath != null){
					List<Integer> totalPath = new ArrayList(last.subList(0, i));
					totalPath.addAll(spurPath);
					
					if(done.add(h(totalPath))){
						B.put(h(totalPath), totalPath);
					}
				}
				
				// restore
				for(int[][] row : g){
					for(int[] u : row){
						if(u[1] < 0){
							u[1] = -u[1]-1;
						}
					}
				}
			}
			
			if(B.isEmpty())break;
			
			double min = Double.POSITIVE_INFINITY;
			List<Integer> arg = null;
			for(List<Integer> p : B.values()){
				double s = len(p, link);
				if(s < min){
					min = s;
					arg = p;
				}
			}
			
			ret.add(min);
			A.add(arg);
			B.remove(h(arg));
		}
		return ret;
	}
	
	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");
//		Thread t = new Thread(null, null, "~", Runtime.getRuntime().maxMemory()){
//			@Override
//			public void run() {
//				long s = System.currentTimeMillis();
//				solve();
//				out.flush();
//				if(!INPUT.isEmpty())tr(System.currentTimeMillis()-s+"ms");
//			}
//		};
//		t.start();
//		t.join();
	}
	
	public static void main(String[] args) throws Exception { new G2().run(); }
	
	private byte[] inbuf = new byte[1024];
	public 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 int[] na(int n)
	{
		int[] a = new int[n];
		for(int i = 0;i < n;i++)a[i] = ni();
		return a;
	}
	
	private long[] nal(int n)
	{
		long[] a = new long[n];
		for(int i = 0;i < n;i++)a[i] = nl();
		return a;
	}
	
	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[][] nmi(int n, int m) {
		int[][] map = new int[n][];
		for(int i = 0;i < n;i++)map[i] = na(m);
		return map;
	}
	
	private int ni() { return (int)nl(); }
	
	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

結果

テストケース

ソースコード