結果

問題 No.2187 三立法和 mod 333
ユーザー sepa38sepa38
提出日時 2023-01-13 22:23:48
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 824 bytes
コンパイル時間 267 ms
コンパイル使用メモリ 87,104 KB
実行使用メモリ 78,508 KB
最終ジャッジ日時 2023-08-26 04:07:47
合計ジャッジ時間 4,560 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 TLE -
testcase_02 TLE -
testcase_03 TLE -
testcase_04 TLE -
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 TLE -
testcase_10 TLE -
testcase_11 TLE -
testcase_12 TLE -
testcase_13 TLE -
testcase_14 TLE -
testcase_15 TLE -
testcase_16 TLE -
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
testcase_22 TLE -
testcase_23 TLE -
testcase_24 TLE -
testcase_25 TLE -
testcase_26 TLE -
testcase_27 TLE -
testcase_28 TLE -
testcase_29 TLE -
testcase_30 TLE -
testcase_31 TLE -
testcase_32 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import bisect
ls = [[] for i in range(333)]
for i in range(1, 4444):
  x = i ** 3
  ls[x%333].append(x*i)
lim = 4444 ** 4
ans = 0
a = int(input())
e3 = [i**3 for i in range(4445)]
e4 = [i**4 for i in range(4445)]
for i in range(1, 4444):
  s3 = e3[i] + e3[i]
  k = (a - s3) % 333
  s4 = e4[i] + e4[i]
  cnt = bisect.bisect_right(ls[k], lim-s4) - bisect.bisect_left(ls[k], e4[i])
  cnt = max(cnt, 0)
  ans += cnt * 3
  if (e3[i] * 3) % 333 == a and (e4[i] * 3) <= lim:
    ans -= 2
  for j in range(i+1, 4444):
    s3 = e3[i] + e3[j]
    k = (a - s3) % 333
    s4 = e4[i] + e4[j]
    if s4 > lim:
      break
    cnt = bisect.bisect_right(ls[k], lim-s4) - bisect.bisect_left(ls[k], e4[j])
    cnt = max(cnt, 0)
    ans += cnt * 6
    if (e3[i] + e3[j] * 2) % 333 == a and (e4[i] + e4[j] * 2) <= lim:
      ans -= 3
print(ans)
0