« Award on Rails | main | ネット復活 »

Genbankファイルから塩基配列だけ抜き出すスクリプト

最近Perl書いてない上に眠れないのでやってみた。

while(<>){s/[^a-z]//g,print if($a=$a?1:/^O/)}

46バイト。
validationだけど、G-languageのパッケージからmgen.gbkを取ってきて
$ perl hoge.pl mgen.gbk | wc
0 1 580074
$ head -1 mgen.gbk
LOCUS NC_000908 580074 bp DNA circular BCT 15-OCT-2004
$

合ってるっぽい。簡単すぎるな。。
コードゴルフってこう言うのなのかなあ。
定石を知りたい。

トラックバック

このエントリーのトラックバックURL:
http://polog.org/mt-tb.cgi/387

コメント (2)

今日nzmさんとこれの話になって,ちゃんと読んでみたら以外と時間かかりました..すげー.

if ($a=$a?1:/^O/)

if ($a=$a?1:/^ORIGIN/)
でフラグ設定だったんですね.

そうそう
このスクリプトで頭ひねったのはそこだけ

コメントを投稿

Powered by
Movable Type 3.34