結果
問題 | No.1241 Eternal Tours |
ユーザー |
![]() |
提出日時 | 2020-09-06 03:37:33 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 1,805 ms / 6,000 ms |
コード長 | 1,497 bytes |
コンパイル時間 | 111 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 21,248 KB |
最終ジャッジ日時 | 2024-11-29 10:56:10 |
合計ジャッジ時間 | 25,832 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 40 |
ソースコード
#!/usr/bin/env python3 X, Y, T, a, b, c, d = list(map(int, input().split())) md = 998244353 T = (T - 1) % (md - 1) + 1 rx, ry = pow(3, (md - 1) // (1 << (X + 1)), mod=md), pow(3, (md - 1) // (1 << (Y + 1)), mod=md) rxi, ryi = pow(rx, -1, mod=md), pow(ry, -1, mod=md) rxa, ryb = pow(rx, a, mod=md), pow(ry, b, mod=md) rxai, rybi = pow(rxa, -1, mod=md), pow(ryb, -1, mod=md) rxc, ryd = pow(rx, c, mod=md), pow(ry, d, mod=md) rxci, rydi = pow(rxc, -1, mod=md), pow(ryd, -1, mod=md) rxpows, rxs2 = [], [] rxpow, rxpowi = 1, 1 rxapow, rxapowi = 1, 1 rxcpow, rxcpowi = 1, 1 for k in range(1 << X): rxpows.append(rxpow + rxpowi) rxs2.append((rxapow - rxapowi) * (rxcpow - rxcpowi) % md) rxpow = rxpow * rx % md rxpowi = rxpowi * rxi % md rxapow = rxapow * rxa % md rxapowi = rxapowi * rxai % md rxcpow = rxcpow * rxc % md rxcpowi = rxcpowi * rxci % md rypows, rys2 = [], [] rypow, rypowi = 1, 1 rybpow, rybpowi = 1, 1 rydpow, rydpowi = 1, 1 for l in range(1 << Y): rypows.append(rypow + rypowi) rys2.append((rybpow - rybpowi) * (rydpow - rydpowi) % md) rypow = rypow * ry % md rypowi = rypowi * ryi % md rybpow = rybpow * ryb % md rybpowi = rybpowi * rybi % md rydpow = rydpow * ryd % md rydpowi = rydpowi * rydi % md ret = 0 for k in range(1 << X): for l in range(1 << Y): fkl = rxpows[k] + rypows[l] + 1 ret += pow(fkl, T, mod=md) * rxs2[k] * rys2[l] % md print(ret * pow(1 << (X + Y + 2), -1, mod=md) % md)