結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
|
提出日時 | 2023-03-10 19:00:50 |
言語 | Go (1.23.4) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,333 bytes |
コンパイル時間 | 14,889 ms |
コンパイル使用メモリ | 226,868 KB |
実行使用メモリ | 34,172 KB |
最終ジャッジ日時 | 2024-09-18 03:21:13 |
合計ジャッジ時間 | 23,096 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | -- * 4 |
other | TLE * 1 -- * 40 |
ソースコード
package mainimport ("bufio""fmt""os")func main() {in := bufio.NewReader(os.Stdin)out := bufio.NewWriter(os.Stdout)defer out.Flush()var n intfmt.Fscan(in, &n)nums := make([]int, n)for i := 0; i < n; i++ {fmt.Fscan(in, &nums[i])}fmt.Fprintln(out, solve(nums))}func solve(nums []int) int {n := len(nums)res := 0seg := NewSegmentTree(nums)for left := 0; left < n; left++ {right := seg.MaxRight(left, func(x int) bool {return x != 1})res += n - right}return res}const INF int = 1e18// PointSetRangeMintype E = intfunc (*SegmentTree) e() E { return 0 }func (*SegmentTree) op(a, b E) E {for b != 0 {a, b = b, a%b}return a}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}type SegmentTree struct {n, size intseg []E}func NewSegmentTree(leaves []E) *SegmentTree {res := &SegmentTree{}n := len(leaves)size := 1for size < n {size <<= 1}seg := make([]E, size<<1)for i := 0; i < n; i++ {seg[i+size] = leaves[i]}for i := size - 1; i > 0; i-- {seg[i] = res.op(seg[i<<1], seg[i<<1|1])}res.n = nres.size = sizeres.seg = segreturn res}func (st *SegmentTree) Get(index int) E {if index < 0 || index >= st.n {return st.e()}return st.seg[index+st.size]}func (st *SegmentTree) Set(index int, value E) {if index < 0 || index >= st.n {return}index += st.sizest.seg[index] = valuefor index >>= 1; index > 0; index >>= 1 {st.seg[index] = st.op(st.seg[index<<1], st.seg[index<<1|1])}}// [start, end)func (st *SegmentTree) Query(start, end int) E {if start < 0 {start = 0}if end > st.n {end = st.n}if start >= end {return st.e()}leftRes, rightRes := st.e(), st.e()start += st.sizeend += st.sizefor start < end {if start&1 == 1 {leftRes = st.op(leftRes, st.seg[start])start++}if end&1 == 1 {end--rightRes = st.op(st.seg[end], rightRes)}start >>= 1end >>= 1}return st.op(leftRes, rightRes)}func (st *SegmentTree) QueryAll() E { return st.seg[1] }// 二分查询最大的 right 使得切片 [left:right] 内的值满足 predicatefunc (st *SegmentTree) MaxRight(left int, predicate func(E) bool) int {if left == st.n {return st.n}left += st.sizeres := st.e()for {for left&1 == 0 {left >>= 1}if !predicate(st.op(res, st.seg[left])) {for left < st.size {left <<= 1if predicate(st.op(res, st.seg[left])) {res = st.op(res, st.seg[left])left++}}return left - st.size}res = st.op(res, st.seg[left])left++if (left & -left) == left {break}}return st.n}// 二分查询最小的 left 使得切片 [left:right] 内的值满足 predicatefunc (st *SegmentTree) MinLeft(right int, predicate func(E) bool) int {if right == 0 {return 0}right += st.sizeres := st.e()for {right--for right > 1 && right&1 == 1 {right >>= 1}if !predicate(st.op(st.seg[right], res)) {for right < st.size {right = right<<1 | 1if predicate(st.op(st.seg[right], res)) {res = st.op(st.seg[right], res)right--}}return right + 1 - st.size}res = st.op(st.seg[right], res)if right&-right == right {break}}return 0}