米田氏が企画した「競プロ典型90問」を解いていきます。
ここでは私の試行錯誤やら解説の理解やら気づきやらをメモしていこうと思います。
競プロ初心者なのでどうか温かい目で見守ってください。
なお、本記事中のPythonコードはブラウザ上で実行可能です。
よろしければ遊んでいってください。
問題
006 - Smallest Subsequence(★5)
![](https://github.com/E869120/kyopro_educational_90/blob/main/problem/006.jpg?raw=true)
kyopro_educational_90/006.jpg at main · E869120/kyopro_educational_90
やってみよう!
テストデータ・テスト関数定義 ↓
# 縮小表示
test_data = [
{
"in":[7, 3, 'atcoder'],
"out": "acd"
},{
"in":[14, 5, 'kittyonyourlap'],
"out": "inlap"
}
]
def test_all(f):
for i, data in enumerate(test_data):
exp = data["out"]
ans = f(*data["in"])
result = "AC" if exp == ans else "WA"
print(f"{i+1} {result}: expected: {exp}, output: {ans}")
下のmain0
関数を完成させて「Run」ボタンをクリックしよう!
def main0(N, K, S):
pass
test_all(main0)
自由欄
try(1st.)
愚直に解いてみます。
def main(N, K, S):
# index:1~Kにansを作る文字のNでのインデックスを格納する
idx = [-1]
for k in range(1, K+1):
char = "~"
for i in range(idx[k-1] + 1, N - (K-k)):
if S[i] < char:
char = S[i]
w_i = i
idx.append(w_i)
return "".join([S[i] for i in idx[1:]])
test_all(main);
結果
一個だけTLE…!
提出 #32635326 - 競プロ典型 90 問
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAALgAAAB/CAYAAABc3gjLAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAAFiUAABYlAUlSJPAAABP5SURBVHhe7Z17jBVFvsdrcAB5DYLIQ3RAwYSHSnjoBnC5WYYLRAcwKmFjIBhAJCJgQBICyEtFIsYIXkgEcSUQ1CDmDhCW6wrcFYVsFBAV3D8Qecj7IoLiiIPMnU9N/w59evrMeUyfOWd6fp+k0qer+lVd3/rVr6q76+SUlmEUJaTUcZaKEkpqlcA3btxohgwZYr7//nsnphzWiZ89e7b57bffnFglDPgKPJYQgLQ33njDLgmZ4KeffjJPPPGEvQ4vpE2ZMsXs3bvXCpZ1if/www/NwIEDzR133GHj4iHn6dWrV9zAdnIuJXvwFfh9991n5s+fb4YPH16h4Pr162cOHTpkbrvtNrNt2zbfSlAdvP766+bMmTNWXO6KdtNNN5mCggLz2Wefmfvvv9+sWbPGin3AgAE2fdSoUbZiIH4RL/k8deqU2bp1q3nggQdsnLvyDB482Hz66afmiy++MHPmzLFxK1assOsff/yxufvuu22ckn3E7WQijvHjx5vRo0ebSZMmObH+yLYUfo8ePZzY9IJL8eKLL5r9+/ebpUuX+lrngwcPmldeecX+poIuX77cTJ8+3QwaNMhWCCrp5MmTTbdu3azwb7zxRrstFfvZZ58133zzjV2vDEROpeN4SvYQZcERCwXs9kURKtYLa+nnErhp1qyZadOmjfnhhx+s2KX5dgf3sd1ghUlnP0F8Y/d55biyHWJE4LQ4iNSvRenSpYuZO3euOX/+vBU3VnjEiBEJi1EteM2lgosyduxYK1CaahGjiMhtwUVoCFOEKE39ggULrCVHDIiAQAsAjzzySMRCIlzx9fGNERJWWFyiCxcu2OO1a9fOrgMViONyfLfQqYibNm2KWHCprF4XS66DeNI5hp+L4q5oSg0GF8WP4uLi0lmzZpX27NmzdM+ePU5sOWXCKy0Tik1nOzdlArVpbCMUFRXZ45Dm5vDhw6WFhYWR43AetmN7kP285xdke7/rcO8rv70h1nEFyaffvt7gzbOSHcQUuCAidIsBoRJHmhfE5E4TcfmJENxClEqFWMosauR3POFwPe7jyzV7K5RUCMkL6SJQd4iVN3Bfr5L9xBW4Fylgln6IiEgXwcYSN4gYvdu4j5MMYnXdFcMteI7JkvRVq1ZVECtpIvBYFSBWYHslu0jqQQ8+M/41PvDQoUOd2Giko3n06NGI704Qv9sLPjO+s3sbfOYysdjOG6MeiYIvzzAguEc0yoRul7fccktkWZZ3c/bsWbseC/oc0oeQzqUb8rl+/frINvFGmZTqJyGBS4fyo48+sgUaS9xAoTPcxqiL32hJPESkjHgw8pHISId0KOks0ll95513ovaj03zzzTebFi1a2HWWOTk55uTJk3Y9HlKxy6x+ROhLliyx+eSc2iHNXmIKXESDsGVExD1KEQuscKtWrayokhG4eyQG4VCR4p1LKh4jH5yPITs/K0prUuZCmObNm9t1llQCefgjIzKE1atX2ziQayKOYUH3sRs0aGBbHUTP/vGGUJUMUe6pXEd8X0JlvnNlJNMRE/84UR9W/Hq5xkTOIUjeZB+/6xQffOfOnXbp7QP47SPXJL67kj3EfZKpKDWZpDqZilLTUIEroUYFroQaFbgSalTgSqhRgSuhRgWuhBoVuBJqoh70FBcXm4sXL5pff/3ViVGU7KZevXomNzfXNG3a1L4+4SVK4MePH7dv2tWpU/sMO++y8CF1WAhbfirj2rVr5ty5c+b22293Yq4TpeSSkpJaKW6lZoNm0a4fqmYl1KjAlVCjAldCjQpcCTUqcCXUqMCVUFNlgW/evNk8+eSTUbNHVRWOxTGZZi0R2I4Zs44cOeLE1Gz4lnXevHn23tZE3NdPmYwbNy5jZVPjLTjiXrt2rTl9+rR5/PHHTZ8+fWxIlzgoKCqTnGffvn1OSu1E7keixqi6CUzgfIHuLniC16qKZXaLjxuTys2RG8uX9OvWrTO7du0yW7ZsMV27drVTR8gX80HC9b/00ktmzJgx9nzLli0zRUVFSc0eEDa458xSwAwHmbLSlZGSwKUJQsQLFy40Bw4csAUtQmPZunVrc++999qluDEwbNgwuy1i4ThMvpOfn2/TvGCZOQfncosIq4m1xmoTmAKCdM4LTHkca6KhqsA1M19L27Zt7Xr37t0j10b+pGJzfVKZp06dauP4vX37dt9tqOAS79fyECfp2WQpuX7ywYStvBaQyDTT1U1KAkc8FCxiHjlypLWaIioKG8EBM7gSx7wlgJV3/0acX331VUQwXjg256BS9O/f34qBm4qwiJfKxDFI52YvXrw4bXN0S6WdOHFilBi//fZbe16uR6z6lStXbBrXhbDJ4yeffGJ/z5w5M2obYF+J54U3gTwRh6VkXwwCcdkAZQidO3e2c8owu5fbEGUDgfvgb7/9tp0Bi2a8ffv2Ng7BIUoKCsHL788//9zOOBVvgh8R9NNPP20efPDBmAVMS0J6ZdawKkjFlhaK1guL2rt3bxvHORH/iRMnIuJF2OzXsmVLG/hNnHsb9gemquPdNyqxsHv37ki+qCzcW/bNNAiZMqRsKF+uHUNDZc8mAhW4NFlYosLCQie2HAqIN9yw2ggVi09hYZ0TtbgidJbcYMSGq4JVxbqRRhBfPF1Qcfm/H/KJRWVWK/LNebHgVeHHH3+MEjhISybBe28zgbS+4kaKy0iFzCYCFThCXblypW8BIIq33nrLLhEnk2NiyVLpDNKZ4YZyg/n3huqC806bNi0iwGPHjlmrfMMNN0QsWSoFLBYZH5YWTVo+wLrTkZYOHMf3VoBMINdKpZaKR4WnortdrExTZYHTQaTjlehNx22gqYUZM2bYJjtZxILiFtC5oSXgmOKa0FrQrAeNzFIrbhDWG1cMEYol4336vLw8u12iSOcT927WrFn25X2BisM5ZAiUCsbfqWQSKhvX6m19pX/1/vvv22U2EPXBw+HDh5N+SR5rvGjRIisyL7gJ0unjpjCiQK2P1xGksvAnURQu7kws5NzgriyyPwWQaHOeiQ8E5DrJI3kNktr0wQOQ3zvvvNNZu06VBR4WVOA1m1gCD3wURUkc6bMELW7lOipwJdSowJVQowJXQo0KXAk1KnAl1KjAlVCjAldCTYUHPYpSU0noSabfRrWBsOW9tpVlrPyqi6KEGhW4EmpU4EqoUYEroUYFroQaFbgSatIucF7qZ/qIjRs3OjHpJxPnVLKTpATOPBhDhgyx3+MhIObCcAfSZK6MVOATtNmzZ9tZkgCBes8xevRoex5Zr+o5MwH540v8sJJN+UtY4FjF+fPnm27dutmPbmHOnDl2shfCihUrbFwsFixYECVUvxtw8OBBu5Qv5YcOHWo/sEXUfNbFOVavXm02bdoU2DmDRFoOzsWSdT++/PLLqG9FpWKzXzZX2FTzJ4axusohCp5kCt99953zqyJLly61AS5cuFBaJrrSnj17RoWyTJUePnzYbiPItkVFRU5MbDi+e7s9e/ZUOIc3varnFCrLeyIUFxeXvvzyy5FrYck68W64tvfee89ZK4d8kxcgfdasWXZZFaqaHy9VyR95S6YsUiFWfhO24JMmTbKB2jdlyhQ7VUQsC04TJbWVeU+YQ8NrTcUNEbAGZ86cMf369XNiynnsscfMc889Z8/FMhapnDNITp06ZZdt2rSJWkq8wPRtd911l7N2Pd/ymJnvNFu1apV17wWlmj/I5AfQSfngCESmhyBjbgGNHz/exrmnS0P0UgH8KoIbCpSC9U4nwVRg8se0r776qv0tzTnn9JLMOYOkzHKZRo0aRaauYMk68QKuyNGjR6PemSC9cePGUfPDtGvXzooim0g1f0AcZUWZVebapIOEBY4fNXfuXDuvyZIlS6IEhJ88ePBg6ytj5ZOFG8NEPn379nVirsOfezLpDf44VlyQ1gN/PN7chtkCBoCJKiubE6YmEyt/aEK0MnnyZGuoKPPqIGGBY5mlOeICxXITmNFo69attgPoHpqTWusOflaX4/7888++b4NxXI7PvrQY8Uj0nEHTrFkzc/ny5UjBsWSdeGHHjh2mQ4cOzlo5pP/yyy9RBY7Fy7Y5TVLNnxcpY3d+00lSPrhYUCaexwcTK+q24FjaHj16RGqsN7AP+7pvzNdff20KCgp8LRvH5fiyrxf6BIRkzxk0Xp9UlhJPs4wgZF3w+txenzxbSDV/XmK5oukiKR9c4OYzVzeT1vgNaeGrMyzEpJGInuE//GaESCcSN2fVqlW2FrP/hg0bzD333OPsHU2iFjyZc6YDfNIJEybYoVSulSXr4rPS+cK3dvvaAvcSF439Ro0aZSeUzzY3JtX8IXwZWiSQT/JbbTijKZZ4Q0sM9TD8VtaxsIHf7777rh2qkyE8ho9YZ+iLISSGvGQYjG1kyEh+Vwbp7C9DUexH4HgybEgoc5+qfM6gh9UyTdjyE49Y+dUvehzClvfaVpax8puSi6IoNQUVuBJqVOBKqFGBK6FGBa6EGhW4EmpU4EqoUYEroabCgx5Fqan4PeipIHD9l7VwoP+yVo66KEqoUYEroUYFroQaFbgSalTgSqhRgSuhRocJHSobVtt5aqdZ/916c7nkshOTGo3qNjLDOww3f27zZycmfcQbJqyJeaqMjA8TLl++3AY3fKnP93kyTwbfS86bN8/06dMnZiC9ur7IFoIQAnAMjpUNhDFPfqRd4IgXEZ89e9Z07949SqwLFy40Bw4csHMd8qU+0yXAyJEjza5du8yyZctM69atzbp16+w68ZkgCCEIqR5L7iP3jXt15MgRJyU1siFP4DVq+/btc1LKjaLEe41jolSLBR84cKCdEWvatGlm5syZVqwEfnft2tVs2bLFfm3NV9mwdu1am6mJEyea06dP24l/WCe+tkIlZwJS7ttrr71mPvjgg2pvydIBsyDwtT35QgdFRUW2MovQRSvgFn+ipF3gTH/AdA6I/KmnnrKTwzBxEIJ1W3Cs08WLF+0+2WbBY9GleRez9D+Wmhm9ZjgxxvS9ta9ZWbDSrBm0xgZ3WizEOsu0d17XjeVDDz1kW0DgnnonCwoKvzw92vFR87f//FskT+O6jnNSKkess1hfBOpufUjv2LFjZCZappxgGjvJdxCkXeBkBguMBX/zzTfN7t27TadOnSI1UwJzrDRt2tTuU1MseJ82fUyD3AYmv0m+FTbiGHHXCPOv0/8yo/5nlA37/2+/s3VsECz5B/J57Ngxuy5zo7CU1g2YS6Zly5aR9CDx5olQcHuB2fz95kieTlw+4WxdOQgWgffu3dvmC+uMsWrfvn0kHS0IlDVg1KjMsh+B31LBkyHtAiczuB+4I1jx7du3m7Zt20YuXAI3AliK6L0WnEC6d3KZTNHxpo7m+C/HzdVrV03nZp1N5+adTUlpidl1qrxJhb8f+bvzKz7cF6AwY4HBoKJT6dNBhTyVhUu/XzIbDm1wtkguT8CMYpQjrkisssNqM0kTswmzDdYeYyjlzm9p4ZIh7QKXZgphM//gokWLzKFDhyK+NxeP+GH//v1RovdacAmp+GJBg2XLq5dnDpw/YE5ePmmF0ax+M1M3p66zRXK4C5SlXx4RNyJ4/vnn02K9/fLUokELJzU1uGb6Cxgp8BMp4n7hhRfs9IBi3bkH7orOb1q2ZKkWgSPsW2+91TRp0sT6WOfOnXNSo+HfIyhghE8FcCPWn/RUmqqgwbI1qdfEPNzhYdP15q6mVcNWpnG9xia3Tq615MlAAVOgY8aMset0JjECbl803eIGvzyRH0SP+JOFsqcsx40bZ60yvnbDhg1tXgQ/cUN+fr69J8KJEydsXLKkXeBkgAkZxYd85plnbG2UziUWmc6mQIeE+GHDhkVcFP4TCLDk7puTKSjsP7X+k/nnD/+M+KVnfj1jWjdsbbYd32YK7yiMdMgW9V3k7BUbBIuo3c338OHDI0Imz1OnTrUFLvcsiKFCN7Hy1LhuY7P37F4z4Z4JkTwl0nEG8kO+3BWyf//+ESGjjenTp9t8uVtpWi/peEoc1lvikqHanmTSNDHT64wZM+yk9gh58eLFNvOSxgSZ/HME24gfRu1mWIybIgXNPzhw44Kksid/Y3eMdX4Fw6q/rHJ+pY94TzJrYp4qI9aTTH1U76ACrxrZKvC0uyhhgPctgiLIY1WFMObJDxV4AvAyURCFyDE4VjYQxjz5oS6KQ7wmvaYRtvzEI5aLogJ3qE0Cv3psm7n677Wm9PdfnJiqk8MQaaeRJje/wImpXmIJXF2UWkjQ4gaOx3GzDRV4LSRocQtBHZchZJ5+86CoqqjAlaxBXuvgz3554h0EKnDF1M3/i8kb+p5p+khRJDQZtMLkFa6xaX7bNBnwXzY+UUS8sV6dBR4A8gSbF66CQgWumJJjO8yljX81vx/5h7lWfN5c3vm8KTn+v06qi2sl5sq/3zcXPxxmfv74GScyMXgyjcB5TYNH795XZ4F3jIJ+z0gFrlQribw6GyQqcCVx6tQ19TuNME0f3mBu7JL8++i4I/FenQ0aFbiSOOKi/Pej5reD5SJNlERenU0HKnAlJjm5DU3DXs/aTmW99gOuW/AUOpnxXp0FXp/lW1ReneUrMNL9PvxIBn2S6VCbnmQWb07fuyMNCjMzR4o+yVRqJSrwWgjvjaSDdB23KqjAayG8FBW0GOVlq2xDfXCH2uSDhxH1wZVaiQpcCTUqcCXURAm8Th3Vu1LzyMnJcX5VpILAr1275qwpSs3gjz/+MLm5uc5aNFGjKJcuXbKPS69everEKEr2g7h5BSAvL8+JuU6UwOHKlSumpKTEeKIVJSvBPalbt66pX7++ExNNBYErSpjQXqUSYoz5fznXT2IesKawAAAAAElFTkSuQmCC)
try(2nd.)
速度改善に資するかは微妙ですが、内側のループの代わりにmin
とindex
で次の文字を探してみます。
def main(N, K, S):
# index:1~Kにansを作る文字のNでのインデックスを格納する
idx = [-1]
chars = list()
for k in range(1, K+1):
begin = idx[k-1] + 1
end = N - (K-k)
char = min(S[begin:end])
i = S.index(char, begin, end)
chars.append(char)
idx.append(i)
return "".join(chars)
test_all(main);
結果
お、AC!やったー
提出 #32651671 - 競プロ典型 90 問
解説・想定コード
一応、解説も見てみます。
すると先程の達成感もつかの間、私の実装では本来TLEとのこと…。とほほ。
kyopro_educational_90/006.jpg at main · E869120/kyopro_educational_90
想定コード
kyopro_educational_90/006.cpp at main · E869120/kyopro_educational_90
ちょっと解説が何言ってるのか分からないのでコードを見ます。
Pythonで書き換えたのが以下。
def main(N, K, S):
'''想定コード(下記URL)参照
https://github.com/E869120/kyopro_educational_90/blob/main/sol/006.cpp
'''
# 前処理
nex = [[N]*26]
for i in range(N-1, -1, -1):
a = nex[-1].copy()
a[ord(S[i])-ord("a")] = i
nex.append(a)
nex = nex[::-1]
# 一文字ずつ貪欲に決める
ans = []
curr_pos = 0
for i in range(K):
for j in range(26):
nex_pos = nex[curr_pos][j]
if N - nex_pos >= K-i:
ans.append(j)
curr_pos = nex_pos + 1
break
return "".join([chr(ord("a")+c) for c in ans])
test_all(main)
結果
提出 #32641757 - 競プロ典型 90 問
は、早い!
上の自力コードでは最大1005msだったものが、120msにまで短縮されました!!
解説の理解
コード全体を見るとようやく理解できました。これ思いつく人、天才かよ。
というより、どうやら頻出のパターンのようです。
この配列は極めて汎用性が高いので、実にさまざまな問題で活用できます!!!
いっそライブラリ化しても良いように思います。
競プロ典型 90 問 006 - Smallest Subsequence(★5) - けんちょんの競プロ精進記録
とは言え、初見ではなかなか分かりづらいので、以下順を追ってみていきます。
各index以降で各文字が最初に登場するindexを格納する配列
やはり言葉だとなんのこっちゃ、となりますので例のkittyonyourlap
で確認します。
N, K, S = test_data[1]["in"]
N, K, S
nex = [[N]*26]
for i in range(N-1, -1, -1):
a = nex[-1].copy()
a[ord(S[i])-ord("a")] = i
nex.append(a)
nex = nex[::-1]
nex
配列を見てみましょう。
# 縮小表示
import string
import pandas as pd
pd.options.display.max_columns = 30
import jinja2 # DataFrame.styleを扱うのに必要。pyodideでは明示的にimportしなければならない
df = pd.DataFrame(
nex, index=[f'S[{i:}]: {c}' for i, c in enumerate(S)]+[""],
columns=list(string.ascii_lowercase))
df.style.set_table_attributes('class="dataframe"').set_td_classes(pd.DataFrame(
[[None]*i+["updated"]+[None]*(25-i) for i in [ord(c)-ord("a") for c in S]]+[[None]*26], index=df.index, columns=df.columns)
).set_table_styles(
[{
'selector': '.updated',
'props': [
('color', 'red'),
('font-weight', 'bold')
]
}]
)
最終的にnex
のサイズはN+1
になりますが、nex[N]
は作成時(の最初の1回)のみ参照されます。
S
の後ろから見て、nex
の下の行から埋めていきます。
S[13]: p
なので、nex[13]
のp
の位置(13)に13を入れます(偶然S
のインデックスと文字のインデックスが一致して分かりづらいです)。他の値はnex[14]
と同じです。
続けてS[12]
を見るとa
なので、nex[12]
のa
の位置(0)に12を入れます。
・・・ということを繰り返すと、『各index以降で各文字が最初に登場するindexを格納する配列』が出来上がります。
貪欲法
今回の場合では辞書順に(アルファベットの"a"から)探索し、条件を満たせば採用する点が貪欲法なんでしょう。
そして条件というのが、「その文字($S_i$)を採用した場合、$S_{i+1}$文字目以降で最終的に長さ$K$の文字列を作れるか」ということになります。
コードに沿って考えてみます。
ans = []
curr_pos = 0
まずnex[0]
を見ます。この時点ではまだ文字を選んでいないので、採用できる文字の$S$におけるインデックスは9以下となります。
nex[0]
を右から(辞書順に)探索して初めて9以下になるのは"i"なので採用します。
# for i in range(K):
i = 0
for j in range(26):
nex_pos = nex[curr_pos][j]
if N - nex_pos >= K-i:
ans.append(j)
curr_pos = nex_pos + 1
break
"".join([chr(ord("a")+c) for c in ans])
次に、"i"のインデックスを一つ進めたnex[2]
を探索し、10以下になるものを探します。
$S[6]$で"n"を見つけたので採用します。
# for i in range(K):
i = 1
for j in range(26):
nex_pos = nex[curr_pos][j]
if N - nex_pos >= K-i:
ans.append(j)
curr_pos = nex_pos + 1
break
"".join([chr(ord("a")+c) for c in ans])
今度はnex[7]
で11以下になるものを探索し・・・ということをK回繰り返せば答えが得られるというわけです。
for i in range(2, K):
for j in range(26):
nex_pos = nex[curr_pos][j]
if N - nex_pos >= K-i:
ans.append(j)
curr_pos = nex_pos + 1
break
"".join([chr(ord("a")+c) for c in ans])
まとめ
今更ながら、アルゴリズム、奥が深い。
今回の解法って一般化して活用することはできるんだろうか?🤔
自由欄