AtCoder Beginner Contest 132 F / Small Products
AtCoder Beginner Contest 132 F / Small Products を解答ACしたので、
行間の証明をメモしておく。
問題はこれ
F - Small Products。
解説
以下、""はガウス記号、" " は集合 の要素の個数とする。
自然数 の積が 以下になるようなものをどうやって求めればいいかを考える必要がある。
一つ Remark をしておく。
各 に対して、 となる の数は、 である。
さて、
何か、 を境目に別れるらしい。
これは、「 かつ 」または、
「 かつ 」が常に満たされるので、そんな気がする。
次のようなことを考えてみる。
を適当な自然数として、
いつ が満たされるかを調べてみる。
であればよい。
これを変形すると、
となる。
つまり、
であれば、 が満たされる。
従って、 であることが分かる。
さらに、正の実数 として、
という対称性から、
となることが分かる。
というわけで、この区間ごとに調べれば良い。
あとは、dp とか累積和とか使っていい感じにやればよい。
解答
解答
Q = 10**9+7 from collections import deque from itertools import accumulate def main(): N, K = map( int, input().split()) ans = 0 d = deque() e = deque() M = int(N**(1/2))*2 for i in range(1, int(N**(1/2))+1): if i**2 == N: d.append(i) M -= 1 continue d.append(i) e.appendleft(N//i) T = list(d)+list(e) CT = [1]*M for i in range(1,M): CT[i] = T[i] - T[i-1] accT = list( map(lambda x:x%Q, accumulate(CT) )) dp = [[0]*M for _ in range(K)] for i in range(M): dp[1][i] = accT[-i-1]*CT[i]%Q for k in range(1, K-1): accT[0] = dp[k][0] for i in range(M-1): accT[i+1] = (accT[i]+dp[k][i+1])%Q for i in range(M): dp[k+1][i] = accT[-i-1]*CT[i]%Q print(sum(dp[K-1])%Q) if __name__ == '__main__': main()