2019-07-13

ICPC2019　国内予選　参加記

noy先輩とHao君と一緒に「UnRated」というチームで出場しました。
3完54位でした。
f:id:wk1080id:20190713022124p:plain

〜前日

チーム練は3回くらいしたと思います。
HCPC 2019 Virtual Vol.20 - 足跡-sokuseki-
ICPC2019模擬国内予選 - 足跡-sokuseki-
four-t practice 2019 Vol.12 - 足跡-sokuseki-

国内突破には4完が必須なので、ABCを早くACしてD以降の考察に時間を割きたいと思っていました。
何度かチーム練して、ABCを自分とnoy先輩で考察&実装、D以降をHao君に考察してもらう戦略を立てました。
自分は今年でICPC引退なので、悔いの残らないようにしたいなあと思って個人でAOJ-ICPCをいっぱいときました。

コンテスト前

15時ごろに大学に着きそうだったので2人にpracticeを任せる。
ライブラリを印刷したやつと蟻本とティッピーを持って部屋に行く。

へへーん pic.twitter.com/cTVjaupDVo
— そすうさ（素数うさぎ） (@wk1080id) July 12, 2019

実装練習にICPC2017国内Cを解く。
チーム練では自分のPC（JIS）を使っていたのですが、突然外付けUSキーボードを使うことになり焦る。
さらに、プリンタを接続して印刷テストしようとしたらPC側でプリンタを認識してくれなくて、これどうするの？使うPC変える？とかいう話になる。（3人はみんなMacBookだけどエディタ等の環境がバラバラすぎるのでア）
いろいろググったりして試行錯誤したところ、開始20分前くらいにうまく認識してくれるようになったので、一安心。
ぞいぞい。

コンテスト

ACを自分が、Bをnoy先輩が、DをHao君が担当することにして、Aから読む。やるだけなので書く。AC。Bの方針を聞いて良さそうなので書いてもらっている間にCを読む。
Cがちょっと難しそうなのでHao君と一緒に考えると、右に載せる分銅と左に載せる分銅をビット全探索すれば0かどうかの判定はできることがわかる。1つ追加して〜のところがよくわかんないなーとか言ってるうちにBがかけたみたいなので見る。answerとprogramをまちがえそうになりながらもAC（複数人で提出チェックしててよかった）。
Cは右と左を全探索すると、そのとき必要になる分銅が一意に決まるので、その重さの分銅がすべてのAiで必要であればそれが答えとなりうるみたいな考察をして、実装する。Hao君に3^Mって言われたけど、計算量的に4^Mでも大丈夫そうだったので、4^Mを書く。

ICPCにはTimeLimitがないので、素早く実装できる方針を立てましょうということです(3^Mはバグらせる自信があって4^Mはバグらせない自信があったので、4^Mをしました
— そすうさ（素数うさぎ） (@wk1080id) July 12, 2019

↑いかにもすんなり通りましたみたいな書き方をしているんですが、実はデバッグ出力と答えを見間違えていて「は？答えが出力されてないやが」とか言いつつ、解法が出てたDに実装を変わったりして、20分ほど溶かした（バカか？）。
DはCを実装している間にDP解が浮かんだらしくて、C終わった後に聞いたけどよくわからなかった（は？）。E以降を読んで解けそうなのを探したりしてると、Dのサンプルが合わないらしいので、Dを一緒に考えることにする。Dは区間add区間minの遅延セグ木でできるんじゃない？とかいう嘘解法をはやしつつ順位表を見てみるとEが0ACでFにACがでているのでFを考えてみる。FはAtCoderぽいな〜（Dもだけど）とか言いつつわからないのでDに戻ってくる。順位表をみるとEもACがでてきたので、noy先輩にEを任せてDの遅延セグ木解のコーナーを探す作業をしてた。コーナーみつからないのでサンプルあうか試そうとおもってコードを書いていると、Eは実装激重やるだけゲーらしいのでPCをnoy先輩に譲る。ここで残り45分くらい。Dは諦めてHao君とFを考えてた。
ラスト1分くらいでEの実装が終わったけど、サンプルが合わず終了。

コンテスト後

4完がアジアボーダーだと思っていたのですが、4完以上のチームが少なく、3完でもアジア行けそうな感じがする。
焼肉を食べた。

放送見ながら打ち上げ〜 pic.twitter.com/mD0WecN36U
— そすうさ（素数うさぎ） (@wk1080id) July 12, 2019

感想とか

Cのバカがなかったら10個くらい順位上がってたので泣いてる。
ICPC引退が（たぶん）半年延びたので嬉しいです。もっと精進したいです。

2019-07-11

four-t practice 2019 Vol.12

競技プログラミング

https://onlinejudge.u-aizu.ac.jp/beta/room.html#fourt_2019_12/info
チームでときました

ABCを自分とnoy先輩、D以降をHao君が読むことにする。
Aを読む。やるだけなので実装する。AC。
Bを書いてもらう。入力形式がアでつらそう。その間Cを読む。なんだかんだBが通る。開始20分くらい。
Cは実は解いたことがあったけど実装することになる。実装つらいんだよね〜とか言いながら書く。その間noy先輩はHao君と一緒にDを考えてもらう。
Cの実装が終わってサンプルを試しても合わないので、Dを書いてもらうことにする（といいつつ、全探索しかわかんないっぽい状況だった）。
デバッグをするとlower_boundとかupper_boundのあたりがまずそうなのでかわってもらう、サンプルがあったので投げる、RE。intじゃだめなことに気づいてもう一回投げる、WA。完全にわからなくなる。ここで1時間20分くらい。
Cがわからなくなったので、Cを捨ててDを一緒に考えることにする。Dは答えが8より大きくならない（1文字ずつ区切った場合、最大が9で最小が1になる）ので、A桁で分割したあとにありうる分割の桁数はA-1, A, A+1 くらいしかないんじゃないかって話になってたんだけど、「0」は使われないのでA>2のときはAでしか分割できないのでは？となって、A>2の場合はNの約数だけ考えることにする。A=2のときどうすればいいかよくわかってなかったんだけど、Hao君がいい感じの案を思いついたらしくて、全部実装してもらうことに。ばちゃ開始から2時間15分後、見事DをAC！！わーい！

反省

一度ACした問題は、バグらせずに解けるようになりましょう（終わった後昔のACコード見たけど同じことやってるんだよね、いまだに今日のコードのだめなところがわかりません、たすけて）

2019-06-30

ICPC2019模擬国内予選

競技プログラミング

問題文：
http://icpc2019.jag-icpc.org/icpc2019/contest/problems_guest_ja.php
順位表：
2019/Practice/模擬国内予選/順位 - ACM-ICPC Japanese Alumni Group

チームUnRatedで出てABC3完75位（有資格者46位）でした。
BCあたり時系列ちょっと怪しいですが、コンテスト中を書きます。

コンテスト

Aを自分が、Bをnoy先輩が、CをHao君が担当することになる。
Aは見たらわかるので提出、AC。Bは場所と沼地の個数をもってbfsらしいので書いてもらう、Cが読解がやばいらしいので一緒に読む。
Bは実装がつらそうで、Cはx軸とy軸独立に考えたら良くて、偶奇で良い感じにできるので、先にCを書く。WA。
Cのコーナーを考えつつBを書いてもらう。BつらそうなのでDを考える。
Dは回転だけ最初にしてそのあと編集距離を求めればできそうで、書きます。投げます。WAになります。ここで1時間経過。
BはそのあといろいろがんばってもらってAC。
Cは実験しようということになって、Cのコードをかいてもらう。20くらいまで実験すると規則性が見えます。Cが通ります。
Dのコーナーがわかって、削除してから回転したほうが得な場合があって、でもそれをどうしたらいいかわからなくなったので、先を読みます。ここでコンテスト半分くらい。

E：Mが小さいので半分全列挙とかするのかな、いや何もわからねえ
F：構文解析マジ？と思ったら文字列数え上げじゃーん、むり
G：一番長いやつから1つ内側にはいると2番目になりませんか？どうなんですか？みたいなことを考えていた
H：グラフだし（部屋、コイン）で拡張ダイクストラしたらできませんか？とか考えてた

なんだかんだ自分にはHが一番できそうだったのでHを考えると、まずコインボタンは途中で押す必要がなく、脱出直前にどの部屋のボタンを押すかを考えればよくて、これは昨日のABC-D。よって、（今いる部屋、通った部屋の数、ワープで得たコインの数）を持ってDPをして、 $d_i$ にきたら、「 $e_i$ - ワープで得たコインの数」を「通った部屋の数」に分割するcombinationを求めて答えに足せばいい、みたいになる。しかし、状態が2000*2000*6000になって配列で持てないので、 $a_i < b_i$ を使ってクエリ先読みして、各部屋でその都度求めることにする。あとdp[2][2000][6000]もつらいのでpriority_queueでがんばる。ここまで考えた時点で残り1時間で、とりあえずコードを書き始めたんだけど、途中でcombinationの引数が1e9までなことに気づいてつらくなる。combination以外を書いて、どうするかかんがえることにする。残り30分くらい。
FができそうらしいのでHao君にバトンタッチする。
combinationの第二引数は高々Nなので、階乗の逆元だけ前計算することにして、分子はO(N)かけて求めて計算機の力を信じることにする。ラスト10分になってFつらそうなので変わってHを書く。サンプルがいくつか合わなくて終了。お疲れ様でした。

反省など

Cですぐ実験できなかったところ
多少時間がかかってもいいから全探索するという選択肢が頭になかった
ABCあたりを高速に通したい

2019-06-20

HCPC 2019 Virtual Vol.20

競技プログラミング

チームで解きました。
onlinejudge.u-aizu.ac.jp

Aを自分が、Bをnoy先輩が、CをHao君が担当することになったので、まずAを読む。
やるだけなので書く。ACしたのでBを書いてもらう。その間Dを読む。
Bが通ったのでCを書いてもらう。Dは拡張ダイクストラをすればできそう（最近拡張ダイクストラいっぱい書いてるので嬉しい）。
Cバグらせてつらそうだったけど、なんとかなってACしたのでDを書く。ライブラリ写経つらい。
Dも難なくACしてこの時点で開始40分くらい。
Dを書いてるうちに、2人がEの考察をしてくれていて、Eはdfsで全探索すればいいらしくて、noy先輩に実装を任せる。
FGHを読んで、Fを考えてみるけど、3^40はつらいなあとか言ってた。HはHao君がフローでできそうと言ってたんだけど、3000*3000頂点がつらいってなって、もうちょっとかんがえると、DPでできそうって言われたので、次元を増やしたらできるんじゃない？とか適当なことを言って考察を任せることにした（最悪）。Gを考えることにする。
Eが実装つらそうだったけどAC。ここで1時間40分。
Hao君にHの実装を任せて、FとGの考察をnoy先輩と2人で考察してわかんないって言うやつをした。一応2時間ばちゃということだったから2時間すぎたあたりで若干感想戦みたいなことを周囲としていた（最悪）。
Hは開始から2時間20分くらいで通りました。遷移がつらそうだったけど、コード長さはそこまで長くない（でも実装できるのつよい）。

初見のセットでのチーム練ははじめてだった*1けど、いい感じかなーっておもいました（小並感）

*1:前回の第一回目は去年の国内予選をした

2019-06-01

AOJ2333 僕の友達は小さい

競技プログラミング AOJ

AOJ-ICPC 500 自力AC！！！！

My friends are small | Aizu Online Judge

問題概要

$N$ 人の友達がいて、 $i$ 番目の友達の重さは $w_i$ である。
この友達を好きな順で選択し、選択した友達の重さの総和が $W$ を超えないようにしたとき、選択した友達の組み合わせは何通りあるか。ただし、選択できる友達が残っている場合はその友達は選択しなければならない。

制約
$1 \le N \le 200$
$1 \le W \le 10^4$
$1 \le w_i \le 10^4$

解法

「選択できる友達が残っている場合はその友達は選択しなければならない」の制約がない場合は普通のナップサック問題となり、

dp[i][j] := iまで見て重さの総和がjのときの場合の数

とすればよい。

「選択できる友達が残っている場合はその友達は選択しなければならない」を言い換えると、「選択した友達の重さの総和 + 選択してない友達の最小の重さ $> W$ を満たさなければならない」となる。
よって、選択していない友達のうち、重さが最小の情報を保持してあげればよい。

dp[i][j][k] := iまで見て重さの総和がjで、選択していない重さ最小の友達のindexがkのときの場合の数

とすると、DP後で $j + w_k > W$ を満たす dp[n][j][k]のみを答えに加算してあげればよい。
計算量は $O(N^2W)$ となって間に合うが、dp[200][10000][200]では配列が足りなくなるので、1次元目を使いまわしてdp[2][10000][200]とすればよい。

AIZU ONLINE JUDGE: Code Review

2019-05-15

AOJ 2642 Dinner

競技プログラミング AOJ

Dinner | Aizu Online Judge
AOJ-ICPC 550点…？

問題概要

$N$ 日間の夕食を食堂か自炊のどちらかで済ませる。
$i$ 日目に食堂で夕食をとると幸福度が $C_i$ 得られる。
各日に自炊をすると得られる幸福度は、 (自炊をした時点での自炊パワー) $\times P$ である。
自炊パワーは、初期値が $Q$ であり、食堂に行くたびに 1 減り、自炊をするたびに 1 増える（食事の終了時に変動する）。
$N$ 日間の幸福度の総和の最大値を求めよ。

制約

$1 \le N \le 500,000$
$0 \le P \le 500,000$
$|Q| , |C_i| \le 500,000$

解法

気づくべき重要な点は、ある日に自炊によって得られる幸福度は、「その日以前に何日自炊し、何日食堂に行ったか」のみに依存するということである。
自炊を $x$ 日行い、食堂に $y$ 日行ったとき、 $x+y+1$ 日目に自炊をすると幸福度は $(Q+x-y)P$ 得られる。

ここで、 $N$ 日間で $K$ 日自炊したときの幸福度を式にしてみる。
$sum := C_0 + C_1 + … + C_{N-1}$ とし、
自炊した日を $a_0$ 日目, $a_1$ 日目, … , $a_{K-1}$ 日目（0-indexed）とする。

$a_i$ 日目以前に自炊した日数が $i$ 日で食堂に行った日数が $a_i - i$ 日なので、 $a_i$ 日目の自炊パワーは、 $Q+i-(a_i-i) = Q - a_i + 2i$ となる。
よって、 $a_i$ 日目に自炊したときの幸福度は $(Q - a_i + 2i)P$ であり、 $K$ 日の自炊で得られる幸福度は $\displaystyle \sum _{i=0}^{K-1} (Q-a_i+2i)P$ である。
一方、食堂で得られる幸福度の総和は $sum - \displaystyle \sum_{i = 0}^{K-1} C_{a_i}$ である。

以上から、自炊した日を $a_0$ 日目, $a_1$ 日目, … , $a_{K-1}$ 日目としたときの幸福度の総和は
$sum - \displaystyle \sum_{i = 0}^{K-1} C_{a_i} + \sum _{i=0}^{K-1} (Q-a_i+2i)P$
となる。
これを変形すると、
$sum - \displaystyle \sum_{i = 0}^{K-1} C_{a_i} + \sum _{i=0}^{K-1} (Q - a_i + 2i)P$
$= sum + \displaystyle \sum _{i=0}^{K-1} \{(Q - a_i + 2i)P-C_{a_i} \}$
$= sum + \displaystyle \sum _{i=0}^{K-1} (QP - a_iP + 2iP-C_{a_i} )$
$= sum + KQP + 2 \frac{K(K-1)}{2} P - \displaystyle \sum _{i=0}^{K-1} ( a_iP + C_{a_i} )$
$= sum + KQP + K(K-1) P - \displaystyle \sum _{i=0}^{K-1} ( a_iP + C_{a_i} )$
となり、これを最大化したいので、 $a_0, a_1, … , a_{K-1}$ は $iP + C_i$ が小さいものから $K$ 個選んでくるのが最適となる。

よって、 $K$ を全探索してあげると良い。

2019-04-28

GitHubで過去のcommit logのメールアドレスを変更する

GitHubのメールアドレス変更が終わった後の話
基本はこれ：https://help.github.com/en/articles/changing-author-info#

変更したい repository の名前を repo とする

ローカルのその repository の外で
git clone --bare https://github.com/user/repo.git
cd repo.git
をして repo.git に入り、その中で

#!/bin/sh

git filter-branch --env-filter '

OLD_EMAIL="your-old-email@example.com"
CORRECT_NAME="Your Correct Name"
CORRECT_EMAIL="your-correct-email@example.com"

if [ "$GIT_COMMITTER_EMAIL" = "$OLD_EMAIL" ]
then
    export GIT_COMMITTER_NAME="$CORRECT_NAME"
    export GIT_COMMITTER_EMAIL="$CORRECT_EMAIL"
fi
if [ "$GIT_AUTHOR_EMAIL" = "$OLD_EMAIL" ]
then
    export GIT_AUTHOR_NAME="$CORRECT_NAME"
    export GIT_AUTHOR_EMAIL="$CORRECT_EMAIL"
fi
' --tag-name-filter cat -- --branches --tags

を書いて、hoge.sh として保存する。
OLD_EMAIL（変更前のアドレス）、CORRECT_NAME（名前）、CORRECT_EMAIL（変更後のアドレス）は書き換える。

sh hoge.shでスクリプトを走らせる。
これでローカルでは変更ができたはず。
（git logで変更されているか確認する。）

git push --force --tags origin 'refs/heads/*'
これでpushできる。あとは
cd ..
rm -rf repo.git
で repo.git を消せば良い。

その後、ローカルの repo を書き換えたいので、 repoのフォルダに入り
git pull --rebaseし、
git logで書き換わっていれば完了。