結果

問題 No.2211 Frequency Table of GCD
ユーザー 箱星箱星
提出日時 2023-02-10 21:58:30
言語 Kotlin
(2.1.0)
結果
AC  
実行時間 1,479 ms / 2,000 ms
コード長 3,735 bytes
コンパイル時間 17,782 ms
コンパイル使用メモリ 450,728 KB
実行使用メモリ 101,868 KB
最終ジャッジ日時 2024-07-07 18:00:24
合計ジャッジ時間 38,540 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

fun PrintWriter.solve() {
    val n = nextInt()
    val m = nextInt()
    val a = Array(n) { nextInt() }
    val freq = Array(m + 1) { 0L }
    val counter = Array(m + 1) { 0L }
    for (v in a) {
        freq[v]++
    }
    for (v in 1..m) {
        if (freq[v] == 0L) continue
        // region divisors
        val number = v.toLong()
        val divs = mutableListOf<Long>()
        var divisor = 1L
        while (divisor * divisor <= number) {
            if (number % divisor == 0L) {
                divs.add(divisor)
                if (divisor * divisor != number) {
                    divs.add(number / divisor)
                }
            }
            divisor++
        }
        // endregion
        for (d in divs) {
            counter[d.toInt()] += freq[v]
        }
    }
    val ans = Array(m + 1) { ModInt(0) }
    for (k in m downTo 1) {
        ans[k] = ModInt(2).pow(counter[k]) - ModInt(1)
        for (l in 2 * k..m step k) {
            ans[k] -= ans[l]
        }
    }
    println(ans.drop(1).joinToString("\n"))
}

// region ModInt
class ModInt(x: Long) {
    companion object {
        const val MOD = 998244353L
        //const val MOD = 1000000007L
    }

    constructor(y: Int) : this(y.toLong())

    val x = (x % MOD + MOD) % MOD

    operator fun plus(other: ModInt): ModInt {
        return ModInt(x + other.x)
    }

    operator fun minus(other: ModInt): ModInt {
        return ModInt(x - other.x)
    }

    operator fun times(other: ModInt): ModInt {
        return ModInt(x * other.x)
    }

    operator fun div(other: ModInt): ModInt {
        return this * other.inv()
    }

    fun pow(exp: Long): ModInt {
        if (exp == 0L) return ModInt(1L)
        var a = pow(exp shr 1)
        a *= a
        if (exp and 1L == 0L) return a
        return this * a
    }

    fun inv(): ModInt {
        return this.pow(MOD - 2)
    }

    override fun equals(other: Any?): Boolean {
        if (this === other) return true
        if (javaClass != other?.javaClass) return false

        other as ModInt

        if (x != other.x) return false

        return true
    }

    override fun hashCode(): Int {
        return x.hashCode()
    }

    override fun toString(): String {
        return "$x"
    }
}

class ModIntArray(n: Int) {
    private val a = LongArray(n) { 0L }

    constructor(n: Int, init: (Int) -> ModInt) : this(n) {
        for (i in 0 until n) {
            a[i] = init(i).x
        }
    }

    operator fun set(i: Int, x: ModInt) {
        a[i] = x.x
    }

    operator fun get(i: Int): ModInt {
        return ModInt(a[i])
    }

    fun joinToString(s: String = ", "): String {
        return a.joinToString(s)
    }
}

fun Long.toModInt() = ModInt(this)

fun Int.toModInt() = ModInt(this)

fun binomSimple(n: Long, k: Long): ModInt {
    var numer = ModInt(1)
    var denom = ModInt(1)
    for (i in 0 until k) {
        numer *= ModInt(n - i)
        denom *= ModInt(k - i)
    }
    return numer / denom
}
// endregion

fun main() {
    Thread(null, {
        val writer = PrintWriter(System.out, false)
        writer.solve()
        writer.flush()
    }, "solve", abs(1L.shl(26)))
        .apply { setUncaughtExceptionHandler { _, e -> e.printStackTrace(); kotlin.system.exitProcess(1) } }
        .apply { start() }.join()
}

// 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