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

提出ソース

結果

問題	No.1320 Two Type Min Cost Cycle
ユーザー	草苺奶昔
提出日時	2023-02-16 12:25:54
言語	Go (1.22.1)
結果	AC
実行時間	872 ms / 2,000 ms
コード長	4,885 bytes
コンパイル時間	11,525 ms
コンパイル使用メモリ	212,820 KB
実行使用メモリ	8,812 KB
最終ジャッジ日時	2023-09-25 17:26:52
合計ジャッジ時間	23,510 ms
ジャッジサーバーID （参考情報）	judge15 / judge14

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 4,376 KB
testcase_01	AC	1 ms 4,376 KB
testcase_02	AC	1 ms 4,380 KB
testcase_03	AC	1 ms 4,376 KB
testcase_04	AC	1 ms 4,376 KB
testcase_05	AC	1 ms 4,376 KB
testcase_06	AC	1 ms 4,376 KB
testcase_07	AC	59 ms 8,248 KB
testcase_08	AC	81 ms 8,336 KB
testcase_09	AC	477 ms 8,736 KB
testcase_10	AC	105 ms 6,712 KB
testcase_11	AC	410 ms 8,700 KB
testcase_12	AC	193 ms 8,572 KB
testcase_13	AC	283 ms 7,180 KB
testcase_14	AC	15 ms 6,756 KB
testcase_15	AC	123 ms 6,736 KB
testcase_16	AC	68 ms 8,560 KB
testcase_17	AC	39 ms 8,304 KB
testcase_18	AC	33 ms 8,304 KB
testcase_19	AC	261 ms 7,276 KB
testcase_20	AC	12 ms 8,212 KB
testcase_21	AC	326 ms 8,812 KB
testcase_22	AC	1 ms 4,376 KB
testcase_23	AC	1 ms 4,376 KB
testcase_24	AC	1 ms 4,376 KB
testcase_25	AC	1 ms 4,380 KB
testcase_26	AC	1 ms 4,380 KB
testcase_27	AC	1 ms 4,380 KB
testcase_28	AC	32 ms 8,308 KB
testcase_29	AC	403 ms 6,732 KB
testcase_30	AC	2 ms 4,376 KB
testcase_31	AC	67 ms 8,244 KB
testcase_32	AC	1 ms 4,380 KB
testcase_33	AC	648 ms 8,292 KB
testcase_34	AC	290 ms 8,340 KB
testcase_35	AC	26 ms 8,588 KB
testcase_36	AC	26 ms 8,512 KB
testcase_37	AC	7 ms 8,248 KB
testcase_38	AC	41 ms 8,564 KB
testcase_39	AC	9 ms 8,444 KB
testcase_40	AC	1 ms 4,376 KB
testcase_41	AC	1 ms 4,376 KB
testcase_42	AC	2 ms 4,380 KB
testcase_43	AC	299 ms 8,336 KB
testcase_44	AC	150 ms 6,800 KB
testcase_45	AC	652 ms 7,320 KB
testcase_46	AC	73 ms 6,756 KB
testcase_47	AC	419 ms 8,280 KB
testcase_48	AC	259 ms 8,364 KB
testcase_49	AC	38 ms 8,340 KB
testcase_50	AC	1 ms 4,376 KB
testcase_51	AC	1 ms 4,380 KB
testcase_52	AC	64 ms 8,472 KB
testcase_53	AC	42 ms 6,916 KB
testcase_54	AC	462 ms 7,084 KB
testcase_55	AC	460 ms 7,064 KB
testcase_56	AC	457 ms 7,076 KB
testcase_57	AC	861 ms 7,116 KB
testcase_58	AC	871 ms 7,116 KB
testcase_59	AC	872 ms 7,112 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 T, N, M int
	fmt.Fscan(in, &T, &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, T == 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
	res := INF

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

		var x int
		if max_ <= 1 {
			x = bfs01Point(n, gi, from, to)
		} else {
			x = dijkstraPoint(n, gi, from, to)
		}
		if x == -1 {
			x = INF
		}
		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 [][]Edge, 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.to, e.weight
			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 [][]Edge, start, target int) int {
	dist := make([]int, n)
	for i := range dist {
		dist[i] = INF
	}
	dist[start] = 0

	pq := NewHeap(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.to, edge.weight
			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 NewHeap(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

結果

テストケース

ソースコード