#856036 (Go) No.1320 Two Type Min Cost Cycle

提出ソース

結果

問題	No.1320 Two Type Min Cost Cycle
ユーザー	草苺奶昔
提出日時	2023-04-02 00:58:31
言語	Go (1.22.1)
結果	AC
実行時間	702 ms / 2,000 ms
コード長	5,063 bytes
コンパイル時間	13,244 ms
コンパイル使用メモリ	225,376 KB
実行使用メモリ	8,356 KB
最終ジャッジ日時	2024-09-24 14:43:47
合計ジャッジ時間	21,355 ms
ジャッジサーバーID （参考情報）	judge1 / judge3

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,816 KB
testcase_01	AC	1 ms 6,812 KB
testcase_02	AC	1 ms 6,944 KB
testcase_03	AC	2 ms 6,940 KB
testcase_04	AC	2 ms 6,940 KB
testcase_05	AC	1 ms 6,940 KB
testcase_06	AC	1 ms 6,944 KB
testcase_07	AC	55 ms 8,096 KB
testcase_08	AC	70 ms 8,196 KB
testcase_09	AC	507 ms 8,316 KB
testcase_10	AC	90 ms 8,148 KB
testcase_11	AC	440 ms 6,940 KB
testcase_12	AC	165 ms 6,944 KB
testcase_13	AC	270 ms 6,940 KB
testcase_14	AC	12 ms 8,064 KB
testcase_15	AC	100 ms 6,944 KB
testcase_16	AC	63 ms 8,268 KB
testcase_17	AC	33 ms 8,076 KB
testcase_18	AC	31 ms 6,944 KB
testcase_19	AC	268 ms 8,324 KB
testcase_20	AC	11 ms 8,084 KB
testcase_21	AC	370 ms 8,272 KB
testcase_22	AC	2 ms 6,944 KB
testcase_23	AC	2 ms 6,940 KB
testcase_24	AC	1 ms 6,940 KB
testcase_25	AC	1 ms 6,940 KB
testcase_26	AC	1 ms 6,940 KB
testcase_27	AC	1 ms 6,944 KB
testcase_28	AC	27 ms 8,068 KB
testcase_29	AC	265 ms 6,940 KB
testcase_30	AC	1 ms 6,940 KB
testcase_31	AC	46 ms 8,064 KB
testcase_32	AC	1 ms 6,944 KB
testcase_33	AC	474 ms 8,204 KB
testcase_34	AC	211 ms 8,204 KB
testcase_35	AC	25 ms 8,276 KB
testcase_36	AC	22 ms 8,268 KB
testcase_37	AC	6 ms 7,880 KB
testcase_38	AC	38 ms 8,280 KB
testcase_39	AC	8 ms 7,880 KB
testcase_40	AC	1 ms 6,940 KB
testcase_41	AC	1 ms 6,940 KB
testcase_42	AC	2 ms 6,940 KB
testcase_43	AC	254 ms 8,204 KB
testcase_44	AC	133 ms 6,940 KB
testcase_45	AC	702 ms 8,304 KB
testcase_46	AC	55 ms 8,076 KB
testcase_47	AC	332 ms 8,188 KB
testcase_48	AC	225 ms 8,184 KB
testcase_49	AC	27 ms 8,152 KB
testcase_50	AC	2 ms 6,944 KB
testcase_51	AC	1 ms 6,944 KB
testcase_52	AC	58 ms 8,236 KB
testcase_53	AC	32 ms 8,232 KB
testcase_54	AC	347 ms 8,344 KB
testcase_55	AC	350 ms 8,248 KB
testcase_56	AC	349 ms 8,336 KB
testcase_57	AC	512 ms 8,328 KB
testcase_58	AC	516 ms 8,256 KB
testcase_59	AC	514 ms 8,356 KB

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

ソースコード

raw source code

package main

import (
	"bufio"
	"fmt"
	"os"
)

func main() {
	in := bufio.NewReader(os.Stdin)
	out := bufio.NewWriter(os.Stdout)
	defer out.Flush()

	var directed, n, m int
	fmt.Fscan(in, &directed, &n, &m)
	edges := make([]Edge, 0, m)
	for i := 0; i < m; i++ {
		var u, v, w int
		fmt.Fscan(in, &u, &v, &w)
		u--
		v--
		edges = append(edges, Edge{u, v, w})
	}
	res := MincostCycle(n, edges, directed == 1)
	if res == INF {
		res = -1
	}
	fmt.Fprintln(out, res)
}

const INF int = 1e18

type Edge struct{ from, to, weight int }

// 返回最小环的权值和. 不存在时返回INF.
//  !把每条边断开,然后求从断开的边的两个端点到另一端点的最短路.
func MincostCycle(n int, edges []Edge, directed bool) int {
	m := len(edges)
	max_ := 0
	for i := range edges {
		if cand := edges[i].weight; cand > max_ {
			max_ = cand
		}
	}
	res := INF

	for i := 0; i < m; i++ {
		e := edges[i]
		cost := e.weight
		from, to := e.to, e.from
		gi := make([][][2]int, n)
		for j := 0; j < m; j++ {
			if i != j {
				e := edges[j]
				gi[e.from] = append(gi[e.from], [2]int{e.to, e.weight})
				if !directed {
					gi[e.to] = append(gi[e.to], [2]int{e.from, e.weight})
				}
			}
		}

		var x int
		if max_ <= 1 {
			x = bfs01Point(n, gi, from, to)
		} else {
			x = dijkstraPoint(n, gi, from, to)
		}
		res = min(res, cost+x)
	}

	return res
}

func min(a, b int) int {
	if a < b {
		return a
	}
	return b
}

func bfs01Point(n int, adjList [][][2]int, start, target int) int {
	dist := make([]int, n)
	for i := range dist {
		dist[i] = INF
	}

	queue := Deque{}
	queue.Append(start)
	dist[start] = 0
	for !queue.Empty() {
		cur := queue.PopLeft()
		if cur == target {
			return dist[cur]
		}
		for _, e := range adjList[cur] {
			next, cost := e[0], e[1]
			if dist[next] > dist[cur]+cost {
				dist[next] = dist[cur] + cost
				if cost == 0 {
					queue.AppendLeft(next)
				} else {
					queue.Append(next)
				}
			}
		}
	}

	return INF
}

type D = int
type Deque struct{ l, r []D }

func (q Deque) Empty() bool {
	return len(q.l) == 0 && len(q.r) == 0
}

func (q Deque) Size() int {
	return len(q.l) + len(q.r)
}

func (q *Deque) AppendLeft(v D) {
	q.l = append(q.l, v)
}

func (q *Deque) Append(v D) {
	q.r = append(q.r, v)
}

func (q *Deque) PopLeft() (v D) {
	if len(q.l) > 0 {
		q.l, v = q.l[:len(q.l)-1], q.l[len(q.l)-1]
	} else {
		v, q.r = q.r[0], q.r[1:]
	}
	return
}

func (q *Deque) Pop() (v D) {
	if len(q.r) > 0 {
		q.r, v = q.r[:len(q.r)-1], q.r[len(q.r)-1]
	} else {
		v, q.l = q.l[0], q.l[1:]
	}
	return
}

func (q Deque) Front() D {
	if len(q.l) > 0 {
		return q.l[len(q.l)-1]
	}
	return q.r[0]
}

func (q Deque) Back() D {
	if len(q.r) > 0 {
		return q.r[len(q.r)-1]
	}
	return q.l[0]
}

// 0 <= i < q.Size()
func (q Deque) At(i int) D {
	if i < len(q.l) {
		return q.l[len(q.l)-1-i]
	}
	return q.r[i-len(q.l)]
}

func dijkstraPoint(n int, adjList [][][2]int, start, target int) int {
	dist := make([]int, n)
	for i := range dist {
		dist[i] = INF
	}
	dist[start] = 0

	pq := nhp(func(a, b H) int {
		return a.dist - b.dist
	}, []H{{start, 0}})

	for pq.Len() > 0 {
		curNode := pq.Pop()
		cur, curDist := curNode.node, curNode.dist
		if cur == target {
			return curDist
		}
		if curDist > dist[cur] {
			continue
		}

		for _, edge := range adjList[cur] {
			next, weight := edge[0], edge[1]
			if cand := curDist + weight; cand < dist[next] {
				dist[next] = cand
				pq.Push(H{next, cand})
			}
		}
	}

	return INF
}

type H = struct{ node, dist int }

// Should return a number:
//    negative , if a < b
//    zero     , if a == b
//    positive , if a > b
type Comparator func(a, b H) int

func nhp(comparator Comparator, nums []H) *Heap {
	nums = append(nums[:0:0], nums...)
	heap := &Heap{comparator: comparator, data: nums}
	heap.heapify()
	return heap
}

type Heap struct {
	data       []H
	comparator Comparator
}

func (h *Heap) Push(value H) {
	h.data = append(h.data, value)
	h.pushUp(h.Len() - 1)
}

func (h *Heap) Pop() (value H) {
	if h.Len() == 0 {
		return
	}

	value = h.data[0]
	h.data[0] = h.data[h.Len()-1]
	h.data = h.data[:h.Len()-1]
	h.pushDown(0)
	return
}

func (h *Heap) Peek() (value H) {
	if h.Len() == 0 {
		return
	}
	value = h.data[0]
	return
}

func (h *Heap) Len() int { return len(h.data) }

func (h *Heap) heapify() {
	n := h.Len()
	for i := (n >> 1) - 1; i > -1; i-- {
		h.pushDown(i)
	}
}

func (h *Heap) pushUp(root int) {
	for parent := (root - 1) >> 1; parent >= 0 && h.comparator(h.data[root], h.data[parent]) < 0; parent = (root - 1) >> 1 {
		h.data[root], h.data[parent] = h.data[parent], h.data[root]
		root = parent
	}
}

func (h *Heap) pushDown(root int) {
	n := h.Len()
	for left := (root<<1 + 1); left < n; left = (root<<1 + 1) {
		right := left + 1
		minIndex := root

		if h.comparator(h.data[left], h.data[minIndex]) < 0 {
			minIndex = left
		}

		if right < n && h.comparator(h.data[right], h.data[minIndex]) < 0 {
			minIndex = right
		}

		if minIndex == root {
			return
		}

		h.data[root], h.data[minIndex] = h.data[minIndex], h.data[root]
		root = minIndex
	}
}

yukicoder

結果

テストケース

ソースコード