結果
問題 | No.1720 Division Permutation |
ユーザー |
|
提出日時 | 2023-02-27 14:39:53 |
言語 | Go (1.23.4) |
結果 |
AC
|
実行時間 | 370 ms / 4,000 ms |
コード長 | 5,817 bytes |
コンパイル時間 | 14,414 ms |
コンパイル使用メモリ | 222,144 KB |
実行使用メモリ | 46,720 KB |
最終ジャッジ日時 | 2024-09-14 22:10:34 |
合計ジャッジ時間 | 29,607 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 60 |
ソースコード
package mainimport ("bufio""fmt""os""strings")const MDO int = 998244353func main() {in := bufio.NewReader(os.Stdin)out := bufio.NewWriter(os.Stdout)defer out.Flush()var n, k intfmt.Fscan(in, &n, &k)perm := make([]int, n)for i := 0; i < n; i++ {fmt.Fscan(in, &perm[i])perm[i]--}tree := NewPermutationTree()root := tree.Build(perm)dp := make([][]int, k+1) // dp[k][i] 表示将perm[0:i]分成k段的方案数for i := 0; i <= k; i++ {dp[i] = make([]int, n+1)}dp[0][0] = 1var dfs func(*Node)dfs = func(node *Node) {if node.IsCut() || node.IsLeaf() {for i := 0; i < k; i++ {dp[i+1][node.right] += dp[i][node.left]dp[i+1][node.right] %= MDO}}sum := make([]int, k)for _, child := range node.children {dfs(child)if node.IsJoin() {for i := 0; i < k; i++ {dp[i+1][child.right] += sum[i]dp[i+1][child.right] %= MDOsum[i] += dp[i][child.left]sum[i] %= MDO}}}}dfs(root)for i := 1; i <= k; i++ {fmt.Fprintln(out, dp[i][n])}}type NodeType intconst (JOIN_ASC NodeType = iotaJOIN_DESCLEAFCUT)type Node struct {kind NodeTypeleft, right int // [left, right) => 顶点对应的区间minV, maxV int // [minV, maxV) => 顶点对应的区间中的最小值和最大值children []*Node}func (n *Node) String() string {mp := map[NodeType]string{JOIN_ASC: "JOIN_ASC",JOIN_DESC: "JOIN_DESC",LEAF: "LEAF",CUT: "CUT",}sb := []string{"Node:\n"}sb = append(sb, fmt.Sprintf("kind : %v\n", mp[n.kind]))sb = append(sb, fmt.Sprintf("left && right : [%v, %v)\n", n.left, n.right))sb = append(sb, fmt.Sprintf("minV && maxV : [%v, %v)\n", n.minV, n.maxV))sb = append(sb, fmt.Sprintf("childrenCount : %v", len(n.children)))return strings.Join(sb, " ")}func (n *Node) Size() int { return n.right - n.left }func (n *Node) IsJoin() bool { return n.kind == JOIN_ASC || n.kind == JOIN_DESC }func (n *Node) IsLeaf() bool { return n.kind == LEAF }func (n *Node) IsCut() bool { return n.kind == CUT }type PermutationTree struct{ Root *Node }func NewPermutationTree() *PermutationTree { return &PermutationTree{} }func (pt *PermutationTree) Build(nums []int) *Node {n := len(nums)desc, asc := []int{-1}, []int{-1}st := []*Node{}seg := newRangeAddRangeMin(make([]int, n))for i := 0; i < n; i++ {for desc[len(desc)-1] != -1 && nums[i] > nums[desc[len(desc)-1]] {seg.Update(desc[len(desc)-2]+1, desc[len(desc)-1]+1, nums[i]-nums[desc[len(desc)-1]])desc = desc[:len(desc)-1]}for asc[len(asc)-1] != -1 && nums[i] < nums[asc[len(asc)-1]] {seg.Update(asc[len(asc)-2]+1, asc[len(asc)-1]+1, nums[asc[len(asc)-1]]-nums[i])asc = asc[:len(asc)-1]}desc = append(desc, i)asc = append(asc, i)t := &Node{kind: LEAF, left: i, right: i + 1, minV: nums[i], maxV: nums[i] + 1}for {kind := CUTif len(st) > 0 {if st[len(st)-1].maxV == t.minV {kind = JOIN_ASC} else if t.maxV == st[len(st)-1].minV {kind = JOIN_DESC}}if kind != CUT {r := st[len(st)-1]if kind != r.kind {r = &Node{kind: kind, left: r.left, right: r.right, minV: r.minV, maxV: r.maxV, children: []*Node{r}}}pt.addChild(r, t)st = st[:len(st)-1]t = r} else if seg.Query(0, i+1-t.Size()) == 0 {t = &Node{kind: CUT, left: t.left, right: t.right, minV: t.minV, maxV: t.maxV, children: []*Node{t}}// do whilefor {pt.addChild(t, st[len(st)-1])st = st[:len(st)-1]if t.maxV-t.minV == t.Size() {break}}for i, j := 0, len(t.children)-1; i < j; i, j = i+1, j-1 {t.children[i], t.children[j] = t.children[j], t.children[i]}} else {break}}st = append(st, t)seg.Update(0, i+1, -1)}pt.Root = st[0]return st[0]}func (pt *PermutationTree) addChild(t, c *Node) {t.children = append(t.children, c)t.left = min(t.left, c.left)t.right = max(t.right, c.right)t.minV = min(t.minV, c.minV)t.maxV = max(t.maxV, c.maxV)}func min(a, b int) int {if a < b {return a}return b}func max(a, b int) int {if a > b {return a}return b}const INF int = 1e18type rangeAddRangeMin struct {n, head intrangeMin, rangeAdd []int}func newRangeAddRangeMin(nums []int) *rangeAddRangeMin {res := &rangeAddRangeMin{}res.n = len(nums)res.head = 1for res.head < res.n {res.head <<= 1}res.rangeMin, res.rangeAdd = make([]int, res.head*2), make([]int, res.head*2)for i := range res.rangeMin {res.rangeMin[i] = INF}copy(res.rangeMin[res.head:], nums)for i := res.head - 1; i > 0; i-- {res.merge(i)}return res}func (r *rangeAddRangeMin) Update(start, end, add int) {r.update(start, end, 1, 0, r.head, add)}func (r *rangeAddRangeMin) Query(start, end int) int {return r.query(start, end, 1, 0, r.head)}func (r *rangeAddRangeMin) Get(pos int) int { return r.Query(pos, pos+1) }func (r *rangeAddRangeMin) query(start, end, pos, left, right int) int {if right <= start || end <= left {return INF}if start <= left && right <= end {return r.rangeMin[pos] + r.rangeAdd[pos]}return min(r.query(start, end, pos*2, left, (left+right)/2), r.query(start, end, pos*2+1, (left+right)/2, right)) + r.rangeAdd[pos]}func (r *rangeAddRangeMin) update(start, end, pos, left, right, add int) {if right <= start || end <= left {return}if start <= left && right <= end {r.rangeAdd[pos] += addreturn}r.update(start, end, pos*2, left, (left+right)/2, add)r.update(start, end, pos*2+1, (left+right)/2, right, add)r.merge(pos)}func (r *rangeAddRangeMin) merge(pos int) {r.rangeMin[pos] = min(r.rangeMin[pos*2]+r.rangeAdd[pos*2], r.rangeMin[pos*2+1]+r.rangeAdd[pos*2+1])}