結果
| 問題 |
No.1036 Make One With GCD 2
|
| コンテスト | |
| ユーザー |
ikd
|
| 提出日時 | 2020-04-25 17:44:02 |
| 言語 | Nim (2.2.0) |
| 結果 |
AC
|
| 実行時間 | 664 ms / 2,000 ms |
| コード長 | 1,482 bytes |
| コンパイル時間 | 4,422 ms |
| コンパイル使用メモリ | 65,964 KB |
| 実行使用メモリ | 35,840 KB |
| 最終ジャッジ日時 | 2024-09-16 13:50:56 |
| 合計ジャッジ時間 | 20,808 ms |
|
ジャッジサーバーID (参考情報) |
judge6 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 41 |
ソースコード
import strutils, sequtils, math
type SegmentTree[T] = object
n: int
dat: seq[T]
e: T
multiply: proc(a, b: T): T
proc initSegmentTree[T](n: int, e: T, f: proc(a, b: T): T): SegmentTree[T] =
let nn = nextPowerOfTwo(n)
var dat = newSeqWith(nn * 2 - 1, e)
return SegmentTree[T](n: nn, dat: dat, e: e, multiply: f)
proc get[T](this: SegmentTree[T], i: int): T =
return this.dat[i + this.n - 1]
proc update[T](this: var SegmentTree[T], i: int, x: T) =
var k = i + this.n - 1
this.dat[k] = x
while k > 0:
k = (k - 1) div 2
this.dat[k] = this.multiply(this.dat[k * 2 + 1], this.dat[k * 2 + 2])
proc query[T](this: SegmentTree[T], ql, qr, i, il, ir: int): T =
if ql <= il and ir <= qr:
return this.dat[i]
if qr <= il or ir <= ql:
return this.e
let m = (il + ir) div 2
return this.multiply(
this.query(ql, qr, i * 2 + 1, il, m),
this.query(ql, qr, i * 2 + 2, m, ir))
proc query[T](this: SegmentTree[T], ql, qr: int): T =
return this.query(ql, qr, 0, 0, this.n)
let read = iterator: string {.closure.} =
for s in stdin.readAll.split:
yield s
proc main() =
let
n = read().parseInt
a = newSeqWith(n, read().parseBiggestInt)
var seg = initSegmentTree[int64](
n,
0,
proc(x, y: int64): int64 = gcd(x, y))
for i in 0..<n:
seg.update(i, a[i])
var
ans: int64 = 0
le = 0
for ri in 1..n:
while le < ri and seg.query(le, ri) == 1:
ans += n - ri + 1
le += 1
echo ans
main()
ikd