結果

問題 No.1661 Sum is Prime (Hard Version)
ユーザー 👑 箱星箱星
提出日時 2021-07-10 19:21:28
言語 Kotlin
(1.9.23)
結果
AC  
実行時間 1,452 ms / 3,000 ms
コード長 2,602 bytes
コンパイル時間 15,133 ms
コンパイル使用メモリ 448,064 KB
実行使用メモリ 127,816 KB
最終ジャッジ日時 2024-05-01 00:52:48
合計ジャッジ時間 34,701 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 306 ms
56,784 KB
testcase_01 AC 306 ms
56,892 KB
testcase_02 AC 1,147 ms
111,260 KB
testcase_03 AC 285 ms
54,332 KB
testcase_04 AC 307 ms
56,836 KB
testcase_05 AC 307 ms
56,788 KB
testcase_06 AC 305 ms
56,832 KB
testcase_07 AC 310 ms
56,776 KB
testcase_08 AC 308 ms
56,808 KB
testcase_09 AC 302 ms
57,032 KB
testcase_10 AC 307 ms
56,904 KB
testcase_11 AC 305 ms
56,756 KB
testcase_12 AC 1,219 ms
111,064 KB
testcase_13 AC 1,120 ms
108,972 KB
testcase_14 AC 1,209 ms
115,036 KB
testcase_15 AC 1,329 ms
119,340 KB
testcase_16 AC 1,204 ms
113,264 KB
testcase_17 AC 1,114 ms
111,236 KB
testcase_18 AC 753 ms
99,700 KB
testcase_19 AC 1,033 ms
108,044 KB
testcase_20 AC 801 ms
96,424 KB
testcase_21 AC 1,015 ms
112,332 KB
testcase_22 AC 1,452 ms
127,816 KB
testcase_23 AC 1,351 ms
123,500 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.PrintWriter
import java.util.*
import kotlin.math.*

fun floorSqrt(n: Long): Long {
    if (n <= 1) return n
    var lo = 1L
    var hi = n
    while (hi - lo > 1) {
        val mid = lo + (hi - lo) / 2
        val sq = mid.toBigInteger().let { it * it }
        if (sq <= n.toBigInteger()) lo = mid
        else hi = mid
    }
    return lo
}

fun primePi(n: Long): Long {
    if (n <= 3) return maxOf(0, n - 1)
    val v = floorSqrt(n).toInt()
    val smalls = (0..v).map { (it + 1) / 2 }.toMutableList()
    var s = (v + 1) / 2
    val roughs = (0 until s).map { 2 * it + 1 }.toMutableList()
    val larges = (0 until s).map { (n / (2 * it + 1) + 1) / 2 }.toMutableList()
    val skip = BooleanArray(v + 1) { false }
    var pc = 0
    for (p in 3..v step 2) {
        if (skip[p]) continue
        val q = p * p.toLong()
        pc++
        if (q * q > n) break
        skip[p] = true
        for (i in q..v step 2L * p) {
            skip[i.toInt()] = true
        }
        var ns = 0
        for (k in 0 until s) {
            val i = roughs[k]
            if (skip[i]) continue
            val d = i * p.toLong()
            val x = if (d <= v) larges[smalls[d.toInt()] - pc] else smalls[(n / d).toInt()].toLong()
            larges[ns] = larges[k] + pc - x
            roughs[ns] = i
            ns++
        }
        s = ns
        var i = v
        for (j in (p..v / p).reversed()) {
            val c = smalls[j] - pc
            val e = j * p.toLong()
            while (i >= e) {
                smalls[i] -= c
                i--
            }
        }
    }
    var res = larges[0] + (s + 2L * (pc - 1)) * (s - 1L) / 2 - larges.slice(1 until s).sum()

    for (l in 1 until s) {
        val q = roughs[l]
        val m = n / q
        val e = smalls[(m / q).toInt()] - pc
        if (e <= l) break
        val t = roughs.slice(l + 1..e).map { smalls[(m / it).toInt()].toLong() }.sum()
        res += t - (e - l) * (pc + l - 1L)
    }
    return res
}

fun PrintWriter.solve() {
    val l = nextLong()
    val r = nextLong()
    println(primePi(r) - primePi(l - 1) + primePi(2 * r) - primePi(2 * l))
}

fun main() {
    val writer = PrintWriter(System.out, false)
    writer.solve()
    writer.flush()
}

// region Scanner
private var st = StringTokenizer("")
private val br = System.`in`.bufferedReader()

fun next(): String {
    while (!st.hasMoreTokens()) st = StringTokenizer(br.readLine())
    return st.nextToken()
}

fun nextInt() = next().toInt()
fun nextLong() = next().toLong()
fun nextLine() = br.readLine()
fun nextDouble() = next().toDouble()
// endregion
0