結果

問題 No.886 Direct
ユーザー 草苺奶昔
提出日時 2023-02-24 18:50:14
言語 Go
(1.23.4)
結果
AC  
実行時間 597 ms / 4,000 ms
コード長 1,331 bytes
コンパイル時間 16,367 ms
コンパイル使用メモリ 219,512 KB
実行使用メモリ 143,140 KB
最終ジャッジ日時 2024-09-13 03:13:25
合計ジャッジ時間 18,449 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

package main
import (
    "bufio"
    "fmt"
    "os"
)
func main() {
    in := bufio.NewReader(os.Stdin)
    out := bufio.NewWriter(os.Stdout)
    defer out.Flush()
    var ROW, COL int
    fmt.Fscan(in, &ROW, &COL)
    res := ROW*(COL-1) + COL*(ROW-1) // 原来相邻的线段数
    n := max(ROW, COL)
    A, B := make([]int, n+2), make([]int, n+2)
    for i := 0; i < ROW; i++ {
        A[i] = ROW - i - 1
    }
    for i := 0; i < COL; i++ {
        B[i] = COL - i - 1
    }
    C := GcdConvolution(A, B)
    fmt.Fprintln(out, (res+C[0]*2)%MOD) // pair of gcd(row,col)=1
}
func max(a, b int) int {
    if a > b {
        return a
    }
    return b
}
const MOD int = 1e9 + 7
// c[k] = ∑a[i]*b[j] mod MOD, gcd(i,j)=k
func GcdConvolution(nums1, nums2 []int) []int {
    n := len(nums1)
    pf := make([]int, n+1)
    copy1, copy2 := make([]int, n+1), make([]int, n+1)
    for i := 0; i < n; i++ {
        copy1[i+1] = nums1[i]
        copy2[i+1] = nums2[i]
    }
    for i := 2; i < n+1; i++ {
        if pf[i] == 0 {
            for j := n / i; j > 0; j-- {
                pf[j*i] = 1
                copy1[j] = (copy1[j] + copy1[j*i]) % MOD
                copy2[j] = (copy2[j] + copy2[j*i]) % MOD
            }
            pf[i] = 0
        }
    }
    res := make([]int, n+1)
    for i := 0; i < n+1; i++ {
        res[i] = copy1[i] * copy2[i] % MOD
    }
    for i := 2; i < n+1; i++ {
        if pf[i] == 0 {
            for j := 1; j < n/i+1; j++ {
                res[j] = (res[j] - res[j*i] + MOD) % MOD
            }
        }
    }
    return res[1:]
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0