結果

問題 No.1728 [Cherry 3rd Tune] Bullet
ユーザー 👑 KazunKazun
提出日時 2021-08-17 04:00:25
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 314 ms / 2,000 ms
コード長 810 bytes
コンパイル時間 316 ms
コンパイル使用メモリ 82,044 KB
実行使用メモリ 76,556 KB
最終ジャッジ日時 2024-10-07 08:49:18
合計ジャッジ時間 2,940 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
52,440 KB
testcase_01 AC 38 ms
51,968 KB
testcase_02 AC 38 ms
51,712 KB
testcase_03 AC 38 ms
51,968 KB
testcase_04 AC 38 ms
52,608 KB
testcase_05 AC 42 ms
57,344 KB
testcase_06 AC 45 ms
59,776 KB
testcase_07 AC 46 ms
57,856 KB
testcase_08 AC 48 ms
59,648 KB
testcase_09 AC 45 ms
57,728 KB
testcase_10 AC 51 ms
61,312 KB
testcase_11 AC 49 ms
59,648 KB
testcase_12 AC 51 ms
60,032 KB
testcase_13 AC 50 ms
60,032 KB
testcase_14 AC 48 ms
58,752 KB
testcase_15 AC 50 ms
59,648 KB
testcase_16 AC 52 ms
61,568 KB
testcase_17 AC 49 ms
59,776 KB
testcase_18 AC 50 ms
60,544 KB
testcase_19 AC 50 ms
60,416 KB
testcase_20 AC 50 ms
60,672 KB
testcase_21 AC 48 ms
59,776 KB
testcase_22 AC 50 ms
59,904 KB
testcase_23 AC 314 ms
76,160 KB
testcase_24 AC 227 ms
76,556 KB
testcase_25 AC 47 ms
57,472 KB
testcase_26 AC 49 ms
59,904 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def Divisors(N):
    N=abs(N)
    L,U=[],[]
    k=1
    while k*k <=N:
        if N%k== 0:
            L.append(k)
            if k*k!=N:
                U.append(N//k)
        k+=1
    return L+U[::-1]

#==================================================
T=int(input())
Mod=10**9+7

Ans=[]
for _ in range(T):
    N,C=map(int,input().split())
    # gcd(N,x) の分布を調べる.
    D=Divisors(N)[::-1]; M=len(D)
    H={}

    for i in range(M):
        d=D[i]
        H[d]=N//d
        for j in range(i):
            e=D[j]
            if e%d==0:
                H[d]-=H[e]

    # σ^k 型の計算
    A=0
    for g in H:
        m=2*g
        A+=H[g]*pow(C,m,Mod); A%=Mod

    # τ σ^k 型の計算
    B=N*pow(C,N,Mod)

    X=(A+B)*pow(2*N,Mod-2,Mod)
    X%=Mod

    Ans.append(X)

print(*Ans,sep="\n")
0