結果
問題 | No.160 最短経路のうち辞書順最小 |
ユーザー |
|
提出日時 | 2020-10-31 15:18:15 |
言語 | Go (1.23.4) |
結果 |
AC
|
実行時間 | 10 ms / 5,000 ms |
コード長 | 10,508 bytes |
コンパイル時間 | 11,122 ms |
コンパイル使用メモリ | 234,060 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-22 05:02:14 |
合計ジャッジ時間 | 12,380 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 26 |
ソースコード
/*URL:https://yukicoder.me/problems/no/160*/package mainimport ("bufio""container/heap""errors""fmt""io""math""os""strconv")var (n, m, s, g intG [200 + 5][]Edge)func main() {defer stdout.Flush()n, m, s, g = readi4()for i := 0; i < m; i++ {a, b, c := readi3()G[a] = append(G[a], Edge{b, Weight{c}})G[b] = append(G[b], Edge{a, Weight{c}})}vinf := Value{INF_B60}less := func(l, r Value) bool { return l.v < r.v }transit := func(cv *Vertex, AG [][]Edge) []*Vertex {res := []*Vertex{}for _, e := range AG[cv.vid] {nv := &Vertex{e.to, Value{cv.v.v + e.w.v}}res = append(res, nv)}return res}SS := []StartPoint{{g, Value{0}}}ds := NewDijkstraSolver(vinf, less, transit)res := ds.Dijkstra(SS, n, G[:n])A := []int{}cid := sfor {A = append(A, cid)if cid == g {break}minId := INF_B60for _, e := range G[cid] {if res[e.to].v+e.w.v == res[cid].v {chmin(&minId, e.to)}}cid = minId}fmt.Println(PrintIntsLine(A...))}// type Value and Weight should be modified according to problems.// DP value typetype Value struct {v int}// weight of edgetype Weight struct {v int}// edge of graphtype Edge struct {to intw Weight}// for initializing start points of dijkstra algorithmtype StartPoint struct {vid intvzero Value}// Less returns l < r, and shared with pq.type Less func(l, r Value) bool// Transit calculates all possible transition.type Transit func(cv *Vertex, AG [][]Edge) []*Vertex// func NewDijkstraSolver(vinf Value, less Less, estimate Estimate) *DijkstraSolver {func NewDijkstraSolver(vinf Value, less Less, transit Transit) *DijkstraSolver {ds := new(DijkstraSolver)// shared with priority queue__less = lessds.vinf, ds.less, ds.transit = vinf, less, transitreturn ds}// verified by [ABC143-E](https://atcoder.jp/contests/abc143/tasks/abc143_e)func (ds *DijkstraSolver) Dijkstra(S []StartPoint, n int, AG [][]Edge) []Value {// initialize datadp, colors := ds._initAll(n)// configure about start points (some problems have multi start points)pq := ds._initStartPoint(S, dp, colors)// body of dijkstra algorithmfor pq.Len() > 0 {pop := pq.pop()colors[pop.vid] = BLACKif ds.less(dp[pop.vid], pop.v) {continue}// to next nodesestimates := ds.transit(pop, AG)for _, es := range estimates {if colors[es.vid] == BLACK {continue}if ds.less(es.v, dp[es.vid]) {dp[es.vid] = es.vpq.push(es)colors[es.vid] = GRAY}}}return dp}type DijkstraSolver struct {vinf Valueless Lesstransit Transit}const (WHITE = 0GRAY = 1BLACK = 2)// _initAll returns initialized dp and colors slices.func (ds *DijkstraSolver) _initAll(n int) (dp []Value, colors []int) {dp, colors = make([]Value, n), make([]int, n)for i := 0; i < n; i++ {dp[i] = ds.vinfcolors[i] = WHITE}return dp, colors}// _initStartPoint returns initialized priority queue, and update dp and colors slices.// *This function update arguments (side effects).*func (ds *DijkstraSolver) _initStartPoint(S []StartPoint, dp []Value, colors []int) *VertexPQ {pq := NewVertexPQ()for _, sp := range S {pq.push(&Vertex{vid: sp.vid, v: sp.vzero})dp[sp.vid] = sp.vzerocolors[sp.vid] = GRAY}return pq}// Less function is shared with a priority queue.var __less Less// Definitions of a priority queuetype Vertex struct {vid intv Value}type VertexPQ []*Vertexfunc NewVertexPQ() *VertexPQ {temp := make(VertexPQ, 0)pq := &tempheap.Init(pq)return pq}func (pq *VertexPQ) push(target *Vertex) {heap.Push(pq, target)}func (pq *VertexPQ) pop() *Vertex {return heap.Pop(pq).(*Vertex)}func (pq VertexPQ) Len() int { return len(pq) }func (pq VertexPQ) Less(i, j int) bool {return __less(pq[i].v, pq[j].v)}func (pq VertexPQ) Swap(i, j int) {pq[i], pq[j] = pq[j], pq[i]}func (pq *VertexPQ) Push(x interface{}) {item := x.(*Vertex)*pq = append(*pq, item)}func (pq *VertexPQ) Pop() interface{} {old := *pqn := len(old)item := old[n-1]*pq = old[0 : n-1]return item}/*******************************************************************//********** common constants **********/const (MOD = 1000000000 + 7// MOD = 998244353ALPH_N = 26INF_I64 = math.MaxInt64INF_B60 = 1 << 60INF_I32 = math.MaxInt32INF_B30 = 1 << 30NIL = -1EPS = 1e-10)/********** bufio setting **********/func init() {// bufio.ScanWords <---> bufio.ScanLinesreads = newReadString(os.Stdin, bufio.ScanWords)stdout = bufio.NewWriter(os.Stdout)}// mod can calculate a right residual whether value is positive or negative.func mod(val, m int) int {res := val % mif res < 0 {res += m}return res}// min returns the min integer among input set.// This function needs at least 1 argument (no argument causes panic).func min(integers ...int) int {m := integers[0]for i, integer := range integers {if i == 0 {continue}if m > integer {m = integer}}return m}// max returns the max integer among input set.// This function needs at least 1 argument (no argument causes panic).func max(integers ...int) int {m := integers[0]for i, integer := range integers {if i == 0 {continue}if m < integer {m = integer}}return m}// chmin accepts a pointer of integer and a target value.// If target value is SMALLER than the first argument,// then the first argument will be updated by the second argument.func chmin(updatedValue *int, target int) bool {if *updatedValue > target {*updatedValue = targetreturn true}return false}// chmax accepts a pointer of integer and a target value.// If target value is LARGER than the first argument,// then the first argument will be updated by the second argument.func chmax(updatedValue *int, target int) bool {if *updatedValue < target {*updatedValue = targetreturn true}return false}// sum returns multiple integers sum.func sum(integers ...int) int {var s ints = 0for _, i := range integers {s += i}return s}// abs is integer version of math.Absfunc abs(a int) int {if a < 0 {return -a}return a}// pow is integer version of math.Pow// pow calculate a power by Binary Power (二分累乗法(O(log e))).func pow(a, e int) int {if a < 0 || e < 0 {panic(errors.New("[argument error]: PowInt does not accept negative integers"))}if e == 0 {return 1}if e%2 == 0 {halfE := e / 2half := pow(a, halfE)return half * half}return a * pow(a, e-1)}/********** FAU standard libraries **********///fmt.Sprintf("%b\n", 255) // binary expression/********** I/O usage **********///str := reads()//i := readi()//X := readis(n)//S := readrs()//a := readf()//A := readfs(n)//str := ZeroPaddingRuneSlice(num, 32)//str := PrintIntsLine(X...)/*********** Input ***********/var (// reads returns a WORD string.reads func() stringstdout *bufio.Writer)func newReadString(ior io.Reader, sf bufio.SplitFunc) func() string {r := bufio.NewScanner(ior)r.Buffer(make([]byte, 1024), int(1e+9)) // for Codeforcesr.Split(sf)return func() string {if !r.Scan() {panic("Scan failed")}return r.Text()}}// readi returns an integer.func readi() int {return int(_readInt64())}func readi2() (int, int) {return int(_readInt64()), int(_readInt64())}func readi3() (int, int, int) {return int(_readInt64()), int(_readInt64()), int(_readInt64())}func readi4() (int, int, int, int) {return int(_readInt64()), int(_readInt64()), int(_readInt64()), int(_readInt64())}// readll returns as integer as int64.func readll() int64 {return _readInt64()}func readll2() (int64, int64) {return _readInt64(), _readInt64()}func readll3() (int64, int64, int64) {return _readInt64(), _readInt64(), _readInt64()}func readll4() (int64, int64, int64, int64) {return _readInt64(), _readInt64(), _readInt64(), _readInt64()}func _readInt64() int64 {i, err := strconv.ParseInt(reads(), 0, 64)if err != nil {panic(err.Error())}return i}// readis returns an integer slice that has n integers.func readis(n int) []int {b := make([]int, n)for i := 0; i < n; i++ {b[i] = readi()}return b}// readlls returns as int64 slice that has n integers.func readlls(n int) []int64 {b := make([]int64, n)for i := 0; i < n; i++ {b[i] = readll()}return b}// readf returns an float64.func readf() float64 {return float64(_readFloat64())}func _readFloat64() float64 {f, err := strconv.ParseFloat(reads(), 64)if err != nil {panic(err.Error())}return f}// ReadFloatSlice returns an float64 slice that has n float64.func readfs(n int) []float64 {b := make([]float64, n)for i := 0; i < n; i++ {b[i] = readf()}return b}// readrs returns a rune slice.func readrs() []rune {return []rune(reads())}/*********** Output ***********/// PrintIntsLine returns integers string delimited by a space.func PrintIntsLine(A ...int) string {res := []rune{}for i := 0; i < len(A); i++ {str := strconv.Itoa(A[i])res = append(res, []rune(str)...)if i != len(A)-1 {res = append(res, ' ')}}return string(res)}// PrintIntsLine returns integers string delimited by a space.func PrintInts64Line(A ...int64) string {res := []rune{}for i := 0; i < len(A); i++ {str := strconv.FormatInt(A[i], 10) // 64bit int versionres = append(res, []rune(str)...)if i != len(A)-1 {res = append(res, ' ')}}return string(res)}// Printf is function for output strings to buffered os.Stdout.// You may have to call stdout.Flush() finally.func printf(format string, a ...interface{}) {fmt.Fprintf(stdout, format, a...)}/*********** Debugging ***********/// debugf is wrapper of fmt.Fprintf(os.Stderr, format, a...)func debugf(format string, a ...interface{}) {fmt.Fprintf(os.Stderr, format, a...)}// ZeroPaddingRuneSlice returns binary expressions of integer n with zero padding.// For debugging use.func ZeroPaddingRuneSlice(n, digitsNum int) []rune {sn := fmt.Sprintf("%b", n)residualLength := digitsNum - len(sn)if residualLength <= 0 {return []rune(sn)}zeros := make([]rune, residualLength)for i := 0; i < len(zeros); i++ {zeros[i] = '0'}res := []rune{}res = append(res, zeros...)res = append(res, []rune(sn)...)return res}