結果
問題 | No.906 Y字グラフ |
ユーザー |
![]() |
提出日時 | 2020-01-06 11:44:30 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 37 ms / 2,000 ms |
コード長 | 789 bytes |
コンパイル時間 | 146 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-11-22 23:55:56 |
合計ジャッジ時間 | 2,548 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 28 |
ソースコード
import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines # a + b + c = N - 1 を数えればよい # どうせ mod 6 ごとに 2 次関数 N = int(read()) MOD = 10 ** 9 + 7 f = [0] * 50 for a in range(1,51): for b in range(1,a+1): for c in range(1,b+1): x = 1 + a + b + c if x >= 50: continue f[x] += 1 def lagrange(a,b,c,x): # f(0) = a, f(1) = b, f(2) = c のとき、f(x) を求める f = a * (x-1) * (x-2) // 2 f += b * x * (2 - x) f += c * x * (x-1) // 2 return f def F(N): if N < 50: return f[N] i = N % 6 + 18 a,b,c = f[i::6][:3] x = (N - i) // 6 return lagrange(a,b,c,x) answer = F(N) % MOD print(answer)