Integer/Summation.py

View this file on GitHub
Last update: 2023-03-18 02:55:12+09:00

Code

#約数のK乗和
def Divisor_Sigma(N, K=1):
    if N==1:
        return 1

    R=1
    p=2
    while p*p<=N:
        if N%p==0:
            e=0
            while N%p==0:
                N//=p
                e+=1

            if K:
                s=pow(p,K)
                R*=(pow(s,e+1)-1)//(s-1)
            else:
                R*=e+1
        p+=1

    if N>1:
        R*=pow(N,K)+1

    return R

Traceback (most recent call last):
  File "/opt/hostedtoolcache/Python/3.13.5/x64/lib/python3.13/site-packages/onlinejudge_verify/documentation/build.py", line 71, in _render_source_code_stat
    bundled_code = language.bundle(stat.path, basedir=basedir, options={'include_paths': [basedir]}).decode()
                   ~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/opt/hostedtoolcache/Python/3.13.5/x64/lib/python3.13/site-packages/onlinejudge_verify/languages/python.py", line 96, in bundle
    raise NotImplementedError
NotImplementedError