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

提出ソース

結果

問題	No.1320 Two Type Min Cost Cycle
ユーザー	草苺奶昔
提出日時	2023-02-16 12:13:13
言語	Go (1.22.1)
結果	AC
実行時間	839 ms / 2,000 ms
コード長	4,769 bytes
コンパイル時間	12,235 ms
コンパイル使用メモリ	212,788 KB
実行使用メモリ	8,836 KB
最終ジャッジ日時	2023-09-25 17:10:18
合計ジャッジ時間	23,718 ms
ジャッジサーバーID （参考情報）	judge11 / judge15

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 4,376 KB
testcase_01	AC	2 ms 4,380 KB
testcase_02	AC	1 ms 4,380 KB
testcase_03	AC	1 ms 4,376 KB
testcase_04	AC	2 ms 4,376 KB
testcase_05	AC	2 ms 4,376 KB
testcase_06	AC	1 ms 4,380 KB
testcase_07	AC	59 ms 8,288 KB
testcase_08	AC	79 ms 8,332 KB
testcase_09	AC	502 ms 8,796 KB
testcase_10	AC	103 ms 6,772 KB
testcase_11	AC	421 ms 8,744 KB
testcase_12	AC	195 ms 8,520 KB
testcase_13	AC	292 ms 7,232 KB
testcase_14	AC	14 ms 6,752 KB
testcase_15	AC	122 ms 6,796 KB
testcase_16	AC	71 ms 8,552 KB
testcase_17	AC	40 ms 8,304 KB
testcase_18	AC	34 ms 8,300 KB
testcase_19	AC	273 ms 8,780 KB
testcase_20	AC	12 ms 8,204 KB
testcase_21	AC	384 ms 8,836 KB
testcase_22	AC	1 ms 4,376 KB
testcase_23	AC	2 ms 4,376 KB
testcase_24	AC	2 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,376 KB
testcase_28	AC	32 ms 8,308 KB
testcase_29	AC	411 ms 6,752 KB
testcase_30	AC	2 ms 4,380 KB
testcase_31	AC	66 ms 8,240 KB
testcase_32	AC	2 ms 4,376 KB
testcase_33	AC	658 ms 8,324 KB
testcase_34	AC	297 ms 8,340 KB
testcase_35	AC	30 ms 8,644 KB
testcase_36	AC	26 ms 8,516 KB
testcase_37	AC	7 ms 8,244 KB
testcase_38	AC	48 ms 8,588 KB
testcase_39	AC	10 ms 8,452 KB
testcase_40	AC	1 ms 4,376 KB
testcase_41	AC	1 ms 4,380 KB
testcase_42	AC	1 ms 4,380 KB
testcase_43	AC	299 ms 8,312 KB
testcase_44	AC	153 ms 6,760 KB
testcase_45	AC	656 ms 7,332 KB
testcase_46	AC	75 ms 6,756 KB
testcase_47	AC	424 ms 8,320 KB
testcase_48	AC	304 ms 6,876 KB
testcase_49	AC	39 ms 8,348 KB
testcase_50	AC	2 ms 4,376 KB
testcase_51	AC	1 ms 4,376 KB
testcase_52	AC	68 ms 8,504 KB
testcase_53	AC	42 ms 8,452 KB
testcase_54	AC	409 ms 7,032 KB
testcase_55	AC	406 ms 7,092 KB
testcase_56	AC	404 ms 7,080 KB
testcase_57	AC	839 ms 7,120 KB
testcase_58	AC	833 ms 8,072 KB
testcase_59	AC	827 ms 7,116 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

// 返回最小环的权值和. 不存在时返回INF.
type Edge struct{ from, to, weight int }

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 = bfs01(n, gi, from)[to]
		} else {
			x = dijkstra(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 bfs01(n int, adjList [][]Edge, start int) (dist []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()
		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
}

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 dijkstra(n int, adjList [][]Edge, start int) (dist []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 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
}

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

結果

テストケース

ソースコード