結果

問題 No.2075 GCD Subsequence
ユーザー taiga0629kyoprotaiga0629kyopro
提出日時 2022-08-19 21:01:07
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,682 ms / 4,000 ms
コード長 3,214 bytes
コンパイル時間 218 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 144,640 KB
最終ジャッジ日時 2024-04-17 02:34:30
合計ジャッジ時間 42,437 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 490 ms
96,000 KB
testcase_01 AC 476 ms
95,616 KB
testcase_02 AC 479 ms
95,744 KB
testcase_03 AC 473 ms
95,872 KB
testcase_04 AC 471 ms
96,128 KB
testcase_05 AC 480 ms
96,000 KB
testcase_06 AC 489 ms
95,872 KB
testcase_07 AC 491 ms
95,744 KB
testcase_08 AC 877 ms
130,688 KB
testcase_09 AC 1,076 ms
144,640 KB
testcase_10 AC 873 ms
130,944 KB
testcase_11 AC 955 ms
137,088 KB
testcase_12 AC 926 ms
133,760 KB
testcase_13 AC 809 ms
126,976 KB
testcase_14 AC 952 ms
137,088 KB
testcase_15 AC 916 ms
128,640 KB
testcase_16 AC 851 ms
130,816 KB
testcase_17 AC 1,025 ms
144,128 KB
testcase_18 AC 2,170 ms
144,256 KB
testcase_19 AC 2,213 ms
144,384 KB
testcase_20 AC 2,682 ms
144,384 KB
testcase_21 AC 2,651 ms
144,512 KB
testcase_22 AC 2,168 ms
144,384 KB
testcase_23 AC 2,188 ms
144,384 KB
testcase_24 AC 2,208 ms
144,512 KB
testcase_25 AC 2,178 ms
144,384 KB
testcase_26 AC 2,174 ms
144,384 KB
testcase_27 AC 2,162 ms
144,384 KB
testcase_28 AC 2,269 ms
140,288 KB
testcase_29 AC 549 ms
95,872 KB
testcase_30 AC 532 ms
95,744 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 


from random import shuffle,randrange
rd=randrange

class primes():
    def __init__(self, n):
        self.prime_num = n
        self.min_prime = [-1] * (self.prime_num + 1)  # 2以上の自然数に対して最小の素因数を表す
        self.min_prime[0] = 0
        self.min_prime[1] = 1
        i = 2
        self.prime = []
        self.memo_prifac = {}
        while i <= self.prime_num:
            if self.min_prime[i] == -1:
                self.min_prime[i] = i
                self.prime.append(i)
            for j in self.prime:
                if i * j > self.prime_num or j > self.min_prime[i]: break
                self.min_prime[j * i] = j
            i += 1

    def prifac(self, n):
        # 素因数分解した結果を返す
        if n in self.memo_prifac:
            return self.memo_prifac[n]
        res = {}
        x = n
        while x > 1:
            p = self.min_prime[x]
            if p in res:
                res[p] += 1
            else:
                res[p] = 1
            x //= p

        # self.memo_prifac[n] = res  #場合によってはこの行を消すと高速化

        return res

    def divisors(self, n):
        # 約数列挙 メモした方がいいかも
        if n== 1: return [1]
        prf = self.prifac(n)
        keys = [key for key in prf]

        def divsearch(i):
            if i == len(keys) - 1:
                return [keys[i] ** j for j in range(prf[keys[i]] + 1)]
            else:
                res = []
                subres = divsearch(i + 1)
                p = keys[i]
                for j in range(prf[p] + 1):
                    for node in subres:
                        res.append(node * p ** j)
                return res

        return divsearch(0)

pri=primes(10**6+100)
u=[0]*(10**6+100)
u[1]=1
for x in range(2,10**6+10):
    u[x]=1
    d=pri.prifac(x)
    for p in d:
        if d[p]>=2:u[x]=0
        u[x]*=-1


from math import gcd
mod=998244353
def naive(n,p):
    ans=0
    for bit in range(1,2**n):
        x=[]
        for i in range(n):
            if (bit>>i)&1:x.append(p[i])
        k=len(x)
        flag=1
        for i in range(k-1):
            if gcd(x[i],x[i+1])==1:flag=0
        ans+=flag
    return ans%mod

def sol1(n,P):
    p=[0]+P[:]
    ma=max(p)
    dp=[0]*(ma+1)
    for i in range(1,n+1):
        res=0
        for j in range(1,ma+1):
            if gcd(p[i],j)==1:res+=dp[j]
        dp[p[i]]+=sum(dp)+1-res
        dp[p[i]]%=mod
    return sum(dp)%mod

def sol2(n,P):
    p=[0]+P[:]
    ma=max(p)
    dp=[0]*(ma+1)
    g=[0]*(ma+1)
    sumdp=0
    for i in range(1,n+1):
        f1=0
        div=[]
        for m in pri.divisors(p[i]):div.append(m)
        for m in div:
            f1+=u[m]*g[m]
        delta=sumdp+1-f1
        dp[p[i]]+=delta
        dp[p[i]]%=mod
        sumdp+=delta
        sumdp%=mod
        for m in div:
            g[m]+=delta
            g[m]%=mod
    return sum(dp)%mod

n=int(input())
a=list(map(int,input().split()))
print(sol2(n,a))

cnt=0
while 0:
    cnt+=1
    print(cnt)
    n=randrange(1,100)
    p=[rd(1,100) for i in range(n)]
    ans1=sol1(n,p)
    ans2=sol2(n,p)
    if ans1!=ans2:
        print(n)
        print(*p)
        exit()
0