参加記です。

予選

• A: C早解きでぎりぎり通過。ただし、剰余のあたりバグらせてしまい若干危うかった。
• B: 本番の時間では参加しませんでしたが、自力で解けたのはBまででした。
• C: (同上)、自力で解けたのはAだけという悲惨さ。
• 予選Aで通って本当に良かった…。

11/8

• 博士予備審査を通す。まあ、退学は予定していましたが、予備審査は通したかったので。
• ちなみに通っていなければ 本選辞退案件 でした…

水曜日〜金曜日

• この時期にバス旅行とか色々とあれ。
• 本選ページを見て、パーカーボーダーに愕然とする。ABCD+E部分点が必要。
• 実は一回スーツケース(含PC)をホテルに忘れかけそうになった(これも 本選辞退案件 )

1日目

朝

• 9:50起床、 Code Festivalの伝統的正装とされるパーカー を来て出発。汐留のサブウェイで食べようと思ったけど中止。結局昼食は歌舞伎揚げで済ませた。
• スタンプラリー、今年はないみたい。スケジュールを巻かなくてよいので楽。
• 席、学年順なんですか？老人を固めるのはやめてくだされ。

本選

• D問題、嘘解法だと思っていたのが想定解法だったらしく、びっくり。今年も運が良かった。
  • UTPCのラミーの経験が役に立った？
• 帰宅してから思ったけどDは一般マッチングであって二部グラフじゃないですよね…

夕方

• tourist氏トークライブ。
  • everything except gifted brain
  • FAR Manager (そこか)
• 体力測定。11ptだった。だめすぎ。
• ペアプロライブ。二人ゲームは難しいです…。でもコードは短いという…
  • http://cf16-pair.contest.atcoder.jp/ 、公開になりませんか?
• エキシビション、さらに難易度上がってませんか。
  • sevenkplus氏が、tourist氏のPCのFAR Managerの使い方がわからず苦戦していた模様(そこか)
• リレー顔合わせ。自己紹介(英語)がぐだってしまった。
  • ショートコードランキングについて。丸写しはまれですが、 写経はそれなりにしてます。すみません。

夜

• 起床フェーズがやばいことは明白なのに今年は宿泊なしですよ？酷いよ？
• 特に理由もなくセガ神楽坂に寄ってから帰宅した。

2日目

朝

• 7:50起床。 コンテストはACだけど開会式はTLE 。

Elimination Tournament

Round1 (group19)

• A問題、後で確認したら、long long案件だった。これは酷い。

Round2 (group10)

Round3 (group5)

まあ、下以外のステッカー、計3枚もらえてうれしかったので、その場で貼ってしまいました。

昼

• 知床鮨。
• LT。きゅうりさんの擬人化の よさみ 。
• 休憩時間にちょろっと退出しUFOキャッチャー。まあ、1回のトライではタオルはおろか何も取れるはずはない。ご苦労様でした。
• りんごさんの挑戦状。おもしろかった。

早解きリレー

• W4yneb0t氏のチームでした。
• Tシャツの色、不幸にも去年と同じ色…。あの色を2着というのは厳しい…。
• ゼッケン。アルファベットにしてくれと去年要望出したのにどうしたー。
• 順位順でGを担当。思いついた解法がどう考えてもTLEだったが、状態を配列に取っておけばいいことを教えてもらった。実装は一発ACだった。
• 9完でFinish。
  • 役に立てなくてすみません。 W4yneb0tさん、Grand Final優勝おめでとうございます。
• Kは、最長パスの接点2の数+葉の数が答えらしい。
  • 実際には少し違った。

帰路

• 他の人の実装見たりしていた。

まとめ

• ゆっくり楽しめました。ただ、ゆっくりすぎて太鼓やらなかった(でもまあ太鼓はメイン機種じゃないのでいいか…)
• お土産少なくないですか。パーカーあるからいいですが…

お土産

• パーカー
  • 白だと汚れるので家宝として封印かなぁ。
• 本選Tシャツ、リレーTシャツ(カーネーション?)
• ステッカー(本選、トーナメント3枚)
• トートバッグ
  • 今回も手提げショルダー両用
• ボールペン、ノート

• 無矛盾な単位系
  • http://www.utpc.jp/2010/
  • http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=2207
• マス目と整数
  • http://code-festival-2016-quala.contest.atcoder.jp/tasks/codefestival_2016_qualA_d

無矛盾な単位系

概要

• 無向グラフが与えられるので、誤った閉路がないかどうか判定せよ。なお点は文字列として与えられる。

方法

• まだ訪れていない点を調べる。
  • 解 説 には深さ優先探索とありますが、実装が適切であれば幅優先探索でも大丈夫っぽいです。ただ、2倍ほど遅くなります。
  • memoに加えるのは、深さ優先探索では入れる時出す時のいずれでも良いですが、幅優先探索だと入れる時に加える必要があります。
  • 訪れていたら、最初に訪れた時の値と同じかどうか確認します。

• なお、独立な単位系が複数存在するかもしれない(連結成分が1個ではないかもしれない)ので、グラフが空になるまで取り出し続ける実装にする必要がある。

• コード(深さ優先探索と幅優先探索をベンチマークするためにdequeにしていますが、深さ優先探索だけで良いならvectorの方が速いです)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40

#include <string>
#include <cstdio>
#include <unordered_map>
#include <deque>
#include <algorithm>
using namespace std;
typedef string key;

char b1[99],b2[99];
void main2(int n){
  unordered_map<key,deque<pair<key,long long>>>m;
  for(;n;n--){
      int d;
      scanf(" 1 %s = 10^%d %s",b1,&d,b2);
      string s1=b1,s2=b2;
      m[s1].emplace_back(s2,d);
      m[s2].emplace_back(s1,-d);
  }
  for(;!m.empty();){
      auto s=m.begin()->first;
      deque<pair<key,long long>>st={ {s,0} };
      unordered_map<key,long long>memo={ {s,0} };
      for(;!st.empty();){
          auto p=*st.rbegin();st.pop_back();
          auto cur=p.first;long long d=p.second;
          for(auto &e:m[cur]){
              if(memo.find(e.first)==memo.end()){
                  st.emplace_back(e.first,d+e.second);
                  memo[e.first]=d+e.second;
              }else if(memo[e.first]!=d+e.second){
                  puts("No");
                  return;
              }
          }
      }
      for(auto &e:memo)m.erase(m.find(e.first));
  }
  puts("Yes");
}
int main(){int n;for(;~scanf("%d",&n)&&n;)main2(n);}

マス目と整数

• 入力を変更し、連結成分重み判定を加えるだけで、マス目と整数の解答に変更することができる。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54

--- utpc2010_10125(aizu2207)_stack.cpp    2016-10-13 04:19:40.000000000 +0900
+++ atcodercodefestival2016qualaD.cpp 2016-10-13 04:42:46.000000000 +0900
@@ -1,22 +1,22 @@
-#include <string>
 #include <cstdio>
 #include <unordered_map>
 #include <deque>
 #include <algorithm>
 using namespace std;
-typedef string key;
+typedef int key;

-char b1[99],b2[99];
-void main2(int n){
+int main(){
+  int r,c,n;
+  scanf("%d%d%d",&r,&c,&n);
  unordered_map<key,deque<pair<key,long long>>>m;
  for(;n;n--){
-      int d;
-      scanf(" 1 %s = 10^%d %s",b1,&d,b2);
-      string s1=b1,s2=b2;
-      m[s1].emplace_back(s2,d);
-      m[s2].emplace_back(s1,-d);
+      int s1,s2,d;
+      scanf("%d%d%d",&s1,&s2,&d);
+      m[s1].emplace_back(s2+r,d);
+      m[s2+r].emplace_back(s1,-d);
  }
  for(;!m.empty();){
+      long long a=-1LL<<61,b=1LL<<61;
      auto s=m.begin()->first;
      deque<pair<key,long long>>st={ {s,0} };
      unordered_map<key,long long>memo={ {s,0} };
@@ -29,12 +29,17 @@
                  memo[e.first]=d+e.second;
              }else if(memo[e.first]!=d+e.second){
                  puts("No");
-                  return;
+                  return 0;
              }
          }
+          if(cur<=r)a=max(a,d);
+          else b=min(b,d);
      }
+      if(a>b){
+          puts("No");
+          return 0;
+      }
      for(auto &e:memo)m.erase(m.find(e.first));
  }
  puts("Yes");
 }
-int main(){int n;for(;~scanf("%d",&n)&&n;)main2(n);}

• https://github.com/cielavenir/procon/commit/073d4d0c991750b056836de69cc7b556864b6353

あとがき

• Code Festival 予選中、無矛盾な単位系が類似問題として頭に浮かんだのですが、当時書いていたコードが 一行入力するごとにグラフを全走査する くそコードだったため、どうしようもできませんでした。
• 結果的には復習になってよかったです。
• ただ、連結成分重み判定はやはり本番中に思いつかなかったと思います…

URL

• https://www.hackerrank.com/contests/tozangezan-birthday-contest/challenges

結果

• ---ox△△◎oo-oo△-ooo△△--△- 1752.66点で2位でした。

A:BonusMondai

• 解いたらあかんやつやこれ
• てか負の配点って設定できちゃだめやろ

B:NA RI KI RI knapsack

• よくある問題だけど得点が謎。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

#include <vector>
#include <algorithm>
#include <cstdio>
using namespace std;
int main(){
  int N,W;
  scanf("%d%d",&N,&W);
  vector<pair<int,long long>>v(N);
  for(int i=0;i<N;i++)scanf("%d%lld",&v[i].first,&v[i].second);
  vector<long long>bag(W+1);bag[0]=1;
  for(auto &e:v){
      for(int j=W;j>=e.first;j--){
          if(bag[j-e.first]&&bag[j]<bag[j-e.first]+e.second)bag[j]=bag[j-e.first]+e.second;
      }
  }
  printf("%lld\n",*max_element(bag.begin(),bag.end())-1);
}

C:Number Guessing Master

• 8桁の数値の当該桁がヒットならACになるかと思ったけどそうでもなかった
• 最初の4ケースがhit、後ろの4ケースがblowだそうです…

D:ポッポカウント(AC)

• 素直に解いて良い。なおNの最大は100であるようだ。

1
2
3
4
5
6
7

#include <stdio.h>
int main(){
  long long N,i,j,R=0;
  scanf("%lld",&N);
  for(i=1;i<N;i++)for(j=i+1;j<=N;j++)R+=__builtin_popcountll(i*j);
  printf("%lld\n",R);
}

E:Hedgehog(TLE)

• 入力拾うだけでTLEになる のでどうしたもんだろう。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

#include <algorithm>
#include <cstdio>
using namespace std;
int main(){
  long long v[1000];
  long long R=0;
  int N,M;
  scanf("%d",&N);
  for(int i=0;i<N;i++){
      scanf("%d",&M);
      for(int j=0;j<M;j++)scanf("%lld",&v[j]);
      std::sort(v,v+M);
      for(int j=1;j<M;j++)R^=v[j]*j;
  }
  printf("%lld\n",R);
}

F:ライツアウト(部分点)

• 5x5を反転させるのでは満点にならないらしい。サイズ、いくらなんだろう。
• 解説読んだけどあれはひどい–;

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

15
0 0
0 1
1 0
1 1
1 3
1 4
2 2
2 3
2 4
3 1
3 2
3 3
4 1
4 2
4 4

G:Surprising Language(部分点)

• Perlで解いたものをBashに埋め込んで210/300点。300点は難読言語専用かな？
• Whitespaceで解くと満点なの?テストケース埋め込みが必要やん…

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

#!/usr/bin/perl
use feature qw(say);
use strict;
my($x,$y)=split(/ /,<>);
my $M=10**9+7;
my @I=(1,1);
my @F=(1,1);
my @R=(1,1);
for(my $i=2;$i<=$x;$i++){
  my $z=($M-($M/$i^0))*$I[$M%$i]%$M;
  $I[$i]=$z;
  $F[$i]=$i*$F[-1]%$M;
  $R[$i]=$z*$R[-1]%$M;
}
say $F[$x]*$R[$x-$y]%$M*$R[$y]%$M;

H:OCR Killer(AC)

• FAです。やった＾ω＾
• http://kjjk.weblio.jp/ の手書きに頑張って書く。途中から疲れてiPadになった(twitter 自分宛DMでPCに転送)。タッチパネル最強＾＾
• 問題文(一部おかしいけど、まあ。これでもGoogle翻訳結果見ながら手直しした方ですよ？)

1
2
3
4
5
6
7
8
9
10
11

달의 구면 상의 러리를 구하심시모
달의 구면 상의 점은 두 실수
(a,b) 로 주어진다 이건 지구의
"동겸a 도, 북왼 b도", 와 같다
서경, 남위는 음수를 사용 찬다
단, 이 몬세에서는 달을 위키백라
기사에 쓰이는 직경의 구라고
가정해야 찬다
(出力)
주머신 두 점의 거리를 츨력한다
오차는 적담히

• 要約

• 東経aと北緯bが2組入力される。月における2点間の距離を求めよ。月の半径はWikiにあるものを使え。

• 答案

• 北緯と東経が(慣習と)逆に入力されることにさえ注意すれば、難しくはない。

1
2
3
4
5
6
7
8
9

#!/usr/bin/python
import sys,math
def calcdist(a):
  lat1,lon1,lat2,lon2=map(lambda e:e*math.pi/180,a)
  return 3474.3/2*math.acos( (math.sin(lat1)*math.sin(lat2)+math.cos(lat1)*math.cos(lat2)*math.cos(lon1-lon2)) )

b1,a1=map(float,sys.stdin.readline().split())
b2,a2=map(float,sys.stdin.readline().split())
print(calcdist([a1,b1,a2,b2]))

I:雑な問題(AC)

• piyo

J:雑な問題2(AC)

• 問題IDが「the-answer-is-to2ange2an」である。
• to2ange2an

K:Classical Cipher

• もういいや

• 鍵がT,Z,Nであることからエニグマだと思ったんですけどね、 ローターの選択とプラグボード配線もエニグマの鍵に含まれますよね？？

• 問題文

1
2
3
4
5
6
7
8

YOUAREGIVENASTRING
WRITEAPROGRAMTODECODETHESTRING
THESTRINGISENCODEDINFOLLOWINGWAY
FIRSTCONVERTEACHCHARACTERBYTHETABLE
THETABLEISGIVENRIGHTBELOW
KSYWQJDRZUIPGOLMXEBTVAFNHC
NEXTCHANGEKTHCHARACTERTOITSKTHNEXTALPHABET
FOREXAMPLEIFAFOURTHCHARACTERISYTHENCHANGEITTOC

• 答案

1
2

#!/usr/bin/ruby
puts gets.chomp.bytes.map.with_index{|b,i|((b-65-i)%26+65).chr}.join.tr('KSYWQJDRZUIPGOLMXEBTVAFNHC','A-Z')

L:洞窟(AC)

• Googleフォームを開いたらソースを開くと答えを見つけることができる。
• SNUKEGORM_BY_DARUDE

M:Matryoshka(AC)

(1)

• CPUに 10時間ほど 頑張って頂きました。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

#!/usr/bin/ruby
s=DATA.gets
314161.times{|i|p i
  s=IO.popen('php','r+b'){|io|
      io.puts s
      io.close_write
      io.read
  }
  puts s if i%1000==0
}
puts s
# <?printf($a=($b=0)?'<?printf($a=($b=%d)?%c%s%c:"%s-%s",$b-1,39,$a,39,dechex($x=%d),dechex($y=%d),$x+$x+$y,-$x);':"bf-cbc",$b-1,39,$a,39,dechex($x=3642),dechex($y=-191),$x+$x+$y,-$x);
# bf-cbc
__END__
<?printf($a='<?printf($a=($b=%d)?%c%s%c:"%s-%s",$b-1,39,$a,39,dechex($x=%d),dechex($y=%d),$x+$x+$y,-$x);',314159,39,$a,39,'','',-1084159067,1084162518);

(2)

• こちらもCPUに10分ほど頑張って頂き、後半戦には入れたのですが、そこで死にました。
• funzip、こういう目的には割と良いツール。

1
2
3
4
5
6

#!/usr/bin/ruby
loop{
  system("<data.zip funzip -bf-cbc >data.bin")
  break if File.binread('data.bin',2)!='PK'
  File.rename('data.bin','data.zip')
}

• ヒントに助けられました。

1
2
3
4
5
6
7
8
9
10

#!/usr/bin/php
<?php
$s='';
for(;;){
  $stdin=trim(fgets(STDIN));
  if(!$stdin)break;
  $s.=$stdin;
}
for($i=0;$i<39;$i++)$s=openssl_decrypt($s,'bf-cbc','bf-cbc');
echo $s;

(3)

• これもヒントに助けられました。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

#!/usr/bin/ruby
require 'net/http'
s=''.dup
Net::HTTP.start('goo.gl'){|http|
  key='xFnLDH'
  location=''
  loop{
      resp=http.get('/'+key)
      h=Hash[*resp]
      break if !h.has_key?('location')
      location=h['location'][0]
      domain,key=location.split('/')[-2..-1]
      s<<domain.scan(/\d/).join
      sleep(1)
  }
  puts location
}
puts s.tr('10','ls')

• Herbert Online Judgeに流し込む。

• 答案

1
2
3

bf-cbc
xFnLDH
SHRIFT

N:Answer the Number(部分点)

• (1)だけです

1
2

def f(n) r=1.0;(1..1/0.0).find{|i|(r=r*(n-i)/n)<=0.5}+1 end
p f(10000)

O:THE EMPTY

• 問題文は透過pngになっているのですが、画像中の暗号を解読できず。
• 0,6,8,9の輪っかでπを表現するとか…うーむ。

1
2
3

#!/usr/bin/ruby
s=Math::PI.to_s.tr('.','')
puts s[gets.to_i-1]

P:THE EMPTY 2(AC)

• 解答言語としてWhitespaceを選択すると問題文が現れます。

1
2
3
4

#!/usr/bin/ruby
require 'prime'
a,b=`dd`.split.map(&:to_i)
p Prime.each(b).drop_while{|e|e<a}.size

Q:THE EMPTY 3(AC)

• Download sample test casesすると、output00.txtに問題文があります。

1
2
3
4
5
6

#!/usr/bin/ruby
require 'prime'
class Integer
  def divisornum() self.prime_division.reduce(1){|s,(n,p)|s*(p+1)} end
end
p (gets.to_i**2).divisornum

R:正しいダジャレギュレーション(AC)

• 落ち着いて問題文を読みましょう。sの種類数と同じ意味です(ただしsが1文字の場合は無効)。

1
2
3
4
5
6
7

#!/usr/bin/ruby
n,a=gets.split.map(&:to_i)
if n<2
  p 0
elsif
  p a**n%100000007
end

S:健康計算機(部分点)

• こういう問題の相場として、複雑な字句解析ってのがあるんですが、どうにもならなかった感じがします。ぐぬぬ。
• まじめ版の字句解析器が通らなかったのは「x'y」という入力のせいですかあーあーそこまで考えてなかったわ

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

def bmi(kg,m)
  kg/m/m
end
a,b,c,d=gets.downcase.scan(/([\d'"\.]+)|([a-zA-Z]+)/).map{|e|e.compact[0]}
if a.scan(/\d/).empty?
  a,b=b,a
end
if c.scan(/\d/).empty?
  c,d=d,c
end
if b.end_with?('g')
  a,c=c,a
  b,d=d,b
end
a=a.to_f
a*=0.3048/12 if b=='in'
a*=0.3048 if b=='ft'
a*=0.9144 if b=='yd'
a*=1000 if b=='km'
a/=100 if b=='cm'
a/=1000 if b=='mm'
c=c.to_f
c*=1000 if d=='t'
c/=1000 if d=='g'
c/=1000000 if d=='mg'
c*=0.45359237 if d=='lb' || d=='lbs'
p bmi(c,a)

T:明日使えない英単語(部分点)

• 知るかいな
• 本番中は2,4に正解。
• alopecoidという単語は「プログレッシブ英和・和英中辞典」「ウィズダム英和・和英辞典」「New Oxford American Dictionary」「Oxford Dictionary of English」のいずれにも載っていない単語ですよ
  • OSXの辞書ってこんなにくそでしたっけ(んなわけはない) ＾ω＾＃

1
2
3

#!/usr/bin/ruby
T=%w(eaglet camaraderie oblast thesaurus alopecoid)
puts T[gets.to_i-1]

U:うんちくビ〜ム

• 解説すら読めない。

V:Side Story of WolfHockeyTeamEasy/Hard

• 問題、読めず。私の国語力も地に落ちたもの。
• jpg埋め込みかよ
• HOUSTON COSMONAUTS

W:Shiritori(部分点)

• 知るかいな(2回目)
• 本番中は1,2,3,4,5,7,8,9に正解できました。
• ウルフサザ?なんですかそれGoogle検索結果0件なんですが…?「sothe」は不正解だしねぇ。

1
2
3

#!/usr/bin/ruby
T=%w(ke-ki ki kitune nezi zi ziyu-nomegami idou u uruhusaza huryoku kuku donou)
puts T[gets.to_i-1]

X:Please ignore.

• 解いたらあかんやつ(2回目)

• 先日iPad Air 9.3.2を脱獄してから久しぶりにいろいろ遊んでます。
• ech0chrome iPhone Repo on Cydia (昔私がビルドしたソフト群)も一部使えました。ncompressとか(tar.Zとか今は使われてるんでしょうか…)
• しかし、nkfやPuTTY、OpenHSPは「Illegal Instruction: 4」と出て動かない。
• 少し調べたところ、http://theiostream.tumblr.com/post/63923259800/patching-iphone-gcc-binaries-to-armv7s に行き当たりました。ldr r3, [r3]命令がarmv6/v7専用命令になっていて、armv7sで解釈できないのが原因とのこと。armv7sで解釈できるように手直しすればよいそうです。crt1.oだけがおかしく、原因は不明とのこと。

#!/bin/sh
for file in "$@";do
sed -i'' 's/\x00\x30\x93\xe4/\x00\x30\x93\xe5/g;s/\x00\x30\xd3\xe4/\x00\x30\xd3\xe5/g;' $file
ldid -S $file
done

• まあ、当時はiPhone 3GSやiPhone 4の時代だったので、v7sのことなんて頭になかったですからねぇ…
• ビルドスクリプト、もしあったら、パッケージ化しようと思います。バックアップは今何処＾＾；(開発機のiPod Touch 3rdはヒューズが飛んだ状態)

• てか、私の出してるiOS gccは間違いなく壊れたバイナリを出力するんですが、どうしたもんか。3.1.2用のtoolchainをセルフコンパイル用にポートしただけなので、新たにtoolchainを作る能力は私にはない＾＾；

• ちなみに、Illegal Instruction: 4はこのパッチで直せますが、mmap error 22は直せません。

• https://forums.developer.apple.com/thread/7524 によると、32bit dylibを読み込むにはMach-Oのヘッダがどうこうという話っぽいです。

• 話は変わりますが、アプリ領域にアクセスするのが不便なのでapplinkなるソフトを作ってましたが、これもiOS8で使えなくなってた模様。com.apple.mobile.installation.plistがなくなったため。
• 突貫工事で新たに作りました。Pythonが必要ではありますが。しかもリンク名がアプリ識別子になる制限付き。
• https://gist.github.com/cielavenir/089ae47301458db8e2e04a67f72132b3

http://abc042.contest.atcoder.jp/tasks/arc058_a http://arc058.contest.atcoder.jp/tasks/arc058_a

概要

• N以上の数値で、与えられた数字を使わずに作れる最小のものを答える。

高速解法

• 使って良い数字の集合(入力の補集合)Aを求めておく。
• まず、与えられた数値の文字列表記Sの前から後ろに向かって、使って良い数字以外が含まれる桁jを求める。
  • もし見つからないなら、Sを出力して終了する。
• 桁jから前側に向かって、Siよりも大きい数字がAに含まれるような桁kを求める。見つからなければ-1とする。
  • 大きい数字がAに含まれるとは、繰り上がりを受け切れるという意味である。
• kが-1ならばAの先頭(ただし先頭が0ならば次の要素)を出力。0以上ならばSの先頭からk文字とSkよりも大きい数字のうち最小のものを出力。
• Aの先頭要素を(Sの文字数-k-1)だけつなげて出力。
  • kを-1にしておくとこの部分での場合分けが不要になるというわけ。
• 2回Sを走査するので、計算量はO(S)あるいはO(logN)となる。

まとめ

• 土曜日、帰宅後に問題ページを開いて最初に思いついた解法がこれです。O(NlogN)解法の方が明らかに実装が楽なので、コンテストとしては負け。
• 復習コンテストのストックにしなかったのは、テストデータ作成が面倒だったからです。

http://abc042.contest.atcoder.jp/submissions/818591

• 1/3、何を思ったか研究室から貸与されたiPadにインストール。まあ、共有物ではないし…
• ULC、面白いシステムだと思う。(最後までプレイすることで完了させられるようになったのは最近らしいけど、それを差し置けば)難しい譜面も繰り返しプレイすることでクリアできるようになる点は良い。
• 節分キャンペーンを進める。所持していない楽曲はフレンド募集スレッドからプレイさせていただくことを学習する。シェアソングはソーシャル要素の答えの１つだと思った。
• 1/15、Break down獲得(RP600突破)。しかし、節分キャンペーンを進めた結果、同時期開催の個人戦報酬を落とす。まあ、ソシャゲを始めた時のイベント報酬が厳しいのと近いので、しかたがない。個人戦の課題曲が全て未所持であることも理由。
• ランダムセレクト祭。ランダム選曲すると過去の曲が常駐するかもしれないというイベント。とりあえずフルコンボ譜面を増やす。
• 1/24、Azitate獲得。条件は難しい難易度であるMASTER譜面をULTIMATEで5曲クリアということで、当時はかなりきつい条件だった。ランダムセレクト祭で引き当てたEternallyが5曲目だった。
• 期間限定RP曲、未確認XX生命体を購入すれば5曲とも所持することになるので、購入。結果、2/8、RP1000突破。NEXT FRONTIER獲得。簡単な難易度であるSTANDARD譜面の難易度値がかなり高いんですが、慣れれば自然と指が動く面白い譜面。
• 2/17、RP1100突破。期間限定RP曲+上位20曲のため、いろいろ埋めないと上がらない。あと、Lobi入れた。
• 2/18、バレンタインキャンペーン完走。きつかったのがMake It Fresh HARDフルゲージと最後の1ページ MASTER Sランク。まあ、クリア出来てよかった(節分キャンペーンリバイバルの最後の報酬を諦めた理由がMASTERフルゲージで、それよりは簡単ですし…)。
• 2/28、RP1200突破。Touch of Goldと雨の音が虹を呼ぶでMASTER ULTできたのが大きい。
• 3/1、期間限定RP曲のうち4曲しか所持していないため、RP1200を割る。
• この進行がハイペースかどうかを訊くのは愚問とのこと。落ち着いて考えればそうなのですが…。ごめんなさい。
• まあ、個人個人自由に楽しむのが第一ですよね。
• あとは、アンロックゲージが溜まったからといってすぐにアンロックチャレンジを開始すべきかという質問も愚問のよう。チケットの消費に慎重な人もいるかもしれないので。

• CROSS X BEATSを始めて2ヶ月、42曲となりました。(RP以外の)ミッション曲を全部回収したこともあり、今後曲がどんどん増えていくということはないでしょう。
• というわけで、時系列でまとめておきます。
• 緑鍵
• 黄鍵
• 祭曲をラ祭以外の方法で入手

• 月別獲得数
• 30,12,4,7,10,19,27,9,37,5,12,5
• 8,2,2,5,8,7,10,7,12,8,10,11

• 8,7,1

• 2年半ありがとうございました！

http://abc032.contest.atcoder.jp/tasks/abc032_b

概要

• sからk文字取り出して、その部分文字列を連想配列に突っ込んで、重複しない数を答える。

高速解法

• ローリングハッシュを生成する。ローリングハッシュには、hash[a:c] = (hash[a:b]*B**(c-b) + hash[b:c]) % Pという性質がある。
• 今回はhash[i:i+k]が得られれば良いので、この式をhash[b:c] = ...と変形すれば良い。
• 部分ハッシュを得るのは、B**(c-b)を事前に計算しておくことでO(1)とできる。
• また、連想配列がハッシュテーブルならば重複判定もO(1)なので、トータルでO(s)となる。
• http://abc032.contest.atcoder.jp/submissions/608838
  • ※C++98互換のためunordered_setではなくsetを使っているためこの実装はO(slogs)

まとめ

• 今回の制約では全くの無駄ですので、解法は賢く選びましょう。
• なぜyukicoder復習コンテストのストックにしなかったし。提出後に気づいたので後の祭り。

• P,Qがどちらも0の場合はX==0とY==0だけを数えます。
• 簡単のため、P,Qの最大公約数を単位とする座標系に変換します。X,Yが変換できなければ弾きます。
• ここで、各方向に進む回数をa,b,c,dと置くと、

1
2

(a+b)P+(c+d)Q==X - (1)
(c-d)P+(a-b)Q==Y - (2)

• となります。
• α=a+b、β=c+dとおきます。
• αP+βQ==1を満たす(α,β)を拡張ユークリッド互除法で求めます。
  • 代表値を(z1,z2)とすると、式(1)を満たすのは、α=z1X+Qn、β=z2X-Pnとなります。
• 式(2)を(β-2d)P+(α-2b)Q==Yと変形し、αとβを代入して、整理します。2*P*d+2*Q*b==(Q*Q-P*P)*n+z2*X*P+z1*X*Q-Yとなります。
• この不定方程式において、 整数解(d,b)を持つようなnが存在する ことが問題が求める条件となります。
  • ここで、上で媒介変数表示により変数を一個減らしたことが効いてきます。
• (Q*Q-P*P)*n+z2*X*P+z1*X*Q-Yがgcd(2*P,2*Q)の倍数であることが条件ですので、-(z2*X*P+z1*X*Q-Y)%gcd(Q*Q-P*P,gcd(2*P,2*Q))==0を判定します。
• なお、このnを用いてd,bを拡張ユークリッド互除法から求めることで、a,cも求まりますので、原点から何回で目標に到達するかを実際に計算することも可能です。
  • プログラムにおいてeとfが可変。abs(a)+abs(b)+abs©+abs(d)が最小となるようなe,fを求めるのもまた一興かと(やってませんが)。

http://yukicoder.me/submissions/65751

参加記、書く予定なかったけど、私のことがコンテストレポートに載る可能性が濃厚なので、書くことにする。

• CodeIQパーカー。Code Festivalでパーカー獲得が叶ったので、今日も願掛けです。
• orisano氏とjmc.miz-miz.bizの話とか(jmc、A問題であふあふなんですが…)。kcs/koj復活お願いします(開発頑張ってください)。
  • てか、yukicoderもkojも、Docker-Builder部分は公開してはどうでしょうか。
  • セキュリティの問題があるので公開はあれですが、問い合わせは可とするのは。どうしてもサイトを作りたい人には、みたく。

• 10位以内のパーカーはないんですね。願掛けはあくまでパーカーに対するものと考えれば前提が(略)
• 懇親会、時間中に3ブースかつ先着でトートバッグとかいろいろひどくないですか。
  • 私は特別に頂けましたが。kenkoooo氏からもらったようなものなので、使うことにします。。

• 食事はあまり摂れませんでした。昼食のサンドイッチを食べたのがコンテスト開始後2時間だったので、まあよかった(サンドイッチ食べたの遅すぎたかなと心配していましたが杞憂に終わった(よくない))。

• 帰り、レジャーランドでmusecaやった。となりにsnuke氏がいた。対戦とかはしてませんが。twitterで募集かけてればまた違ったかも。そういえばmusecaって現在は対戦機能ないんだった。

• そういえばインタビューでシエルについて訊かれたんだっけ。スカイガンナー、機会があったらやってみてください。キャラは可愛いのにゲーム内容がかなりガチだったりするので＾＾；