他人のHaskell日記 RSSフィード

Haskell初心者が、リハビリがてらに「ふつける」と「入門Haskell」片手に、試行錯誤するサイト。

2006-06-17

transposeは非汚染  transposeは非汚染 - 他人のHaskell日記 を含むブックマーク

transposeは縦横共に無限リストが扱える

*Main > take 5 $ map (take 5) $ transpose [[1..],[100..]]
[[1,100],[2,101],[3,102],[4,103],[5,104]]

やるなλ星人っ


IOモナドもprimitive  IOモナドもprimitive - 他人のHaskell日記 を含むブックマーク

http://www.sampou.org/haskell/report-revised-j/standard-prelude.html

うーん、自分が知りたいところが分からなくて歯がゆい


seqの定義はprimitive seqの定義はprimitive - 他人のHaskell日記 を含むブックマーク

http://www.sampou.org/haskell/report-revised-j/standard-prelude.html

seq :: a -> b -> b
seq = ...       -- Primitive

これはつまり、逐次的な評価は、

他の関数の組み合わせだけでは表現が出来ないという事だろうか。


関数も型定義できる? 関数も型定義できる? - 他人のHaskell日記 を含むブックマーク

Haskell ReportのPreludeより

 type  ReadS a  = String -> [(a,String)]

別名だからデータコンストラクタは宣言できない(ってデータじゃなくて関数だから当たり前か)

・・なのにどうやって使うんだろ

strHead ::(Read a) => ReadS a
strHead (x:xs) = [(read [x],xs)]
  where types =(x::Char,xs::String)

うーん・・確かに関数の型宣言に使える。

「ReadS a」における型変数aの存在理由が分からなかったので

その型変数を使うような関数を作ってみる

strHeadAs :: (Read a,Num a) => a -> ReadS a
strHeadAs n (x:xs) = [(read [x] `asTypeOf` n,xs)]
  where types =(x::Char,xs::String)

実行結果

*Main> strHeadAs (0::Int) "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Int, String)]
*Main> strHeadAs (0::Integer) "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Integer, String)]

確かに多相的になっている。

なんだか使いこなすのは大変そうだ

あっ・・部分適用して「関数を返す」ようにすれば便利なのでは

strHeadAsInt     = strHeadAs (0::Int)
strHeadAsInteger = strHeadAs (0::Integer)

実行結果

*Main> strHeadAsInt "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Int, String)]
*Main> strHeadAsInteger "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Integer, String)]

うん、動く。

あああ、つまりまり、もしかして

*Main> (strHead::ReadS Int) "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Int, String)]
*Main> (strHead::ReadS Integer) "1abc"
[(1,"abc")]
*Main> :t it
it :: [(Integer, String)]

なるほど、strHeadAsなんていらない!

型変数恐ろしい子


クラスメソッドがわからない クラスメソッドがわからない - 他人のHaskell日記 を含むブックマーク

クラスメソッドの中では引数を二つとるタイプ関数が定義されているけど

その引数がともに多相型で、しかも別の種類だった場合、

そのクラスメソッドはどちらの多相型クラスに書くべきなの?(というか書けるの?)


素数列をつくろう 素数列をつくろう - 他人のHaskell日記 を含むブックマーク

http://ll.jus.or.jp/2006/blog/doukaku1

面白そうなのでミラーテストバージョンを書いてみました

http://www.faireal.net/articles/7/03/#d30129

このサイトを参考にしました。

import List

main =  putStrLn $unwords $ map show $ primes 100
primes n = 2:filter (check.map head.group.modList)[3..n] where
  check (1:[]) =  True
  check (1:(-1):_) = True
  check _ = False
  modList n = map (g.(`mod` n).(2^)) (f (n-1)) where
    f x | even x = x:f (x`div`2)
        | otherwise = x:[]
    g s | s > 1 = s - n 
        | otherwise = s

実行結果

*Main> main
2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97

素数狩り 擬素数狩り - 他人のHaskell日記 を含むブックマーク

素数のフリしてお高くとまっている偽素数どもを発見し晒し上げる

正しい素数列としてid:nobsunのエラトステネスのふるいを

(http://haskell.g.hatena.ne.jp/nobsun/20060617)

フェルマーテストとしてd:id:nurseの書いたものをhaskellに書き換えたものを使います。

(http://d.hatena.ne.jp/nurse/20060616)

import List

--エラトステネスのふるい(id:nobsun)
sieve (p:ps) = p : sieve [q | q <- ps, q `mod` p /= 0]
primes = sieve [2..]

--フェルマーテスト(id:nurse)のHaskell版
primes' = 2:filter(\x -> ((2^(x-1)) `mod` x) == 1)[3..]

--ミラーテスト(id:taninsw)
primes'' = 2:filter millerTest [3..] 
millerTest = check.map head.group.modList where
  check (1:[]) =  True
  check (1:(-1):_) = True
  check _ = False
  modList n = map (g.(`mod` n).(2^)) (f (n-1)) where
    f x | even x = x:f (x`div`2)
        | otherwise = x:[]
    g s | s > 1 = s - n 
        | otherwise = s

--テスト用コード

--それぞれの素数列のn番目の要素を抜き出す
(!!!) n =  map ($ n)  [(primes !!),(primes' !!),(primes'' !!)]

--「正しい素数列」と「擬素数を含む素数列」から「(素数の振りをした)擬素数の配列」を作る
unmatch x@(xh:xt) y@(yh:yt) 
  | xh >  yh = yh:unmatch x yt  
  | xh <  yh = xh:unmatch xt y  
  | xh == yh = unmatch xt yt

実行結果

*Main> take 10 $ unmatch primes primes'                       --擬素数(最初の10個)
[341,561,645,1105,1387,1729,1905,2047,2465,2701]
*Main> take 10 $ unmatch primes primes''                      --強擬素数(最初の10個)
[2047,3277,4033,4681,8321,15841,29341,42799,49141,52633]

*Main> takeWhile (<10000) $ unmatch primes primes'            --擬素数(10000以下)
Loading package haskell98-1.0 ... linking ... done.
[341,561,645,1105,1387,1729,1905,2047,2465,2701,2821,3277,4033,4369,4371,4681,5461,6601,7957,8321,8481,8911]
*Main> takeWhile (<10000) $ unmatch primes primes''           --強擬素数(10000以下)
[2047,3277,4033,4681,8321]

ミラーテストで現れる擬素数は、フェルマーテストでも擬素数になっています。

逆にいえば、フェルマーテストを通過できなかったものは、ミラーテストも通過できない、ということです。

ちなみに、フェルマーテストミラーテストも中身はまったく理解していません。

コードの短さ コードの短さ - 他人のHaskell日記 を含むブックマーク

print "2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97";

より

main=print$takeWhile(<101)$s[2..]where s(n:p)=n:s[q|q<-p,q`mod`n/=0]

の方が短い。

フェルマーテストならもっと短い。

main=print$2:filter(\x->2^(x-1)`mod`x==1)[3..100]

どうして誰もベーシックを使わないんだー素数計算BASICどうして誰もベーシックを使わないんだー素数計算BASIC版 - 他人のHaskell日記 を含むブックマーク

Haskell日記じゃない・・。)

N88-Basic フェルマーテストバージョン

10 CLS
20 PRINT 2;
30 FOR I = 3 TO 100
40   IF (2^(I-1)) MOD I = 1 THEN PRIN I
50 NEXT I

う・・オーバーフローした・・。

とにかく自分より小さい数で割っちまえばいいんだ!

100 CLS
110 DIM I,J
120 PRINT 2;
130 FOR I = 3 TO 100
140 FOR J = 2  TO SQR(I)
150 IF I MOD J = 0 THEN GOTO 180
160 NEXT J
170 PRINT I;
180 NEXT I

f:id:taninsw:20060618004909j:image

昔はQuickBasicとVisualBasic人間だった。

NScripter素数計算 NScripter版素数計算 - 他人のHaskell日記 を含むブックマーク

(もはやLLでは無い・・・。)

使ったことがないのでリファレンスを読みながら作成。

最初はfor文があることに気付かなかった。

雰囲気はなんとなくアセンブラっぽいですね。

*define
game
*start
caption "100までの素数列を作ろう。"
“2” /
for %0=3 to 100
   mov %4,0
   mov %5,%0
   dec %5
   for %1=2 to %5
     mov %2,%0
     mod %2,%1
     if %2=0 mov %4,1
   next
   if %4=1 jumpf
  “ /
   %0 /
   ” /
   ~
next
click
end 

f:id:taninsw:20060618013237j:image

Excel素数計算。 Excelで素数計算。 - 他人のHaskell日記 を含むブックマーク

(・・・すでに言語ですらない。)

VBSは禁じ手にしてセルを利用して求める方針。

Excelも使ったことがないので四苦八苦。(Lotus 1-2-3なら弄ったことある)

アドイン入れないと最大公倍数関数GCDが使えないし

ヘルプみようとしてもイルカ邪魔するしで大変。

セルの相対指定絶対指定に戸惑ったもののなんとか出来た。

	2	3		5		7				11
		13				17		19		
		23						29		
31						37				
41		43				47						53						59		61		
				67				71		
73						79				
83						89				
				97			

f:id:taninsw:20060618042834j:image

R1 C1~R100 C1 は順番に1~100

R1 C1~R1 C100 は順番に1~100

R2 C2~R2 C100 は「OR(GCD(RC1,R1C)=1,GCD(RC1,R1C)=R1C)」

 つまり2との最大公約数が1かまたはその値になるときTrue

全体

R3 C2~R100 C100 は「AND(OR(OR(GCD(RC1,R1C)=1),RC1>=R1C),R[-1]C)」

 つまり「一つ上のセルがTrue」かつ、一番左のセルの値と一番上のセルの値の公倍数が、

 一番上のセルの値あるいは1になるときTrue

結果表示部

R101 C1~R101は「IF(R[-1]C,R1C,"")」

 つまり「一つ上のセルがTrue」ならば値を表示し、そうでなければ表示しない

極めてシンプル

素数を求めるプログラムBATファイル版 素数を求めるプログラムBATファイル版 - 他人のHaskell日記 を含むブックマーク

・・を作ろうと弄ってたら朝になってしまった・・。

for /L %%I in (2,1,100) do echo "ttt" > %%1.prime

ファイルを作って削除していくことでフルいにかけようと思って

ファイル生成してからDIRしたらどうしても順番が変わってしまう。

ファイル名で並び替えてもダメ、日時で並び替えてもダメ

うーん・・・巨大なファイルでも生成させるべきかな・・・

・・寝よう



独り言 独り言 - 他人のHaskell日記 を含むブックマーク

brainfuckやwhitespaceで作る人が出てこないかな.

正規表現で出来るみたい(http://blog.livedoor.jp/dankogai/archives/50534089.html)

だからオートマトンやパーサーでも出来そうだ。

専用のLL言語作って1バイトで出来るとか言ったらやっぱ怒られるんだろうな。

DavianDavian2012/01/06 16:14God help me, I put aside a whole afternoon to fuigre this out.

udveweudvewe2012/01/06 19:2621ObXs <a href="http://szbmaiurkppr.com/">szbmaiurkppr</a>

pjpdsaocpjpdsaoc2012/01/07 22:36jPUNSS , [url=http://dwepjvfdsppz.com/]dwepjvfdsppz[/url], [link=http://ckcdiyjoding.com/]ckcdiyjoding[/link], http://oxjlddkdkqxu.com/

ywxbdkqywxbdkq2012/01/09 23:22fRxEP3 <a href="http://fyiftgcnicxn.com/">fyiftgcnicxn</a>

qpjdaenmqpjdaenm2012/01/11 03:44hHRgBW , [url=http://yrdxnqbyzdzz.com/]yrdxnqbyzdzz[/url], [link=http://tulnkdnbqovm.com/]tulnkdnbqovm[/link], http://djxnfmrqmwqd.com/

HedHed2013/07/27 04:52Weeeee, what a quick and easy sotnilou.

OmarOmar2013/07/30 01:33Wow, this is in every recpset what I needed to know.

MikiMiki2013/08/02 22:23This was so helpful and easy! Do you have any artielcs on rehab? http://mtlhgz.com [url=http://hsyjlrwhamf.com]hsyjlrwhamf[/url] [link=http://zskglfaz.com]zskglfaz[/link]

DennisRourbDennisRourb2017/03/28 07:50metronidazol creme bijsluiter
<a href= http://jadalipfor.e-monsite.com >metronidazol hond bijsluiter</a>
<a href=http://lonsibounm.eklablog.net>propranolol</a>
strattera add inattentive
<a href= http://bisacodyl.tilda.ws/bisacodylhtml >bisacodyl 5 mg zetpil</a>
<a href=http://brandimroz.bloggsite.se>flagyl ovules kopen</a>
flagyl 500 mg bg

WilliamGlushWilliamGlush2017/03/30 08:33cipralex 10 mg nebenwirkungen
<a href= http://nouw.com/rematague >gabapentin 300 und alkohol</a>
<a href=http://nouw.com/rematague>gabapentin micro labs 400mg</a>
azithromycin 500 hec
<a href= https://agnusdinap.yolasite.com >nifurantin und pille</a>
<a href=http://lonnielenn.edublogs.org/2017/03/27/azithromycin/>azithromycin 500 1a pharma in der schwangerschaft</a>
gabapentin teva

EdwardBidgeEdwardBidge2017/03/30 09:10flagyl 400 metronidazole
<a href= http://antionehoy.e-monsite.com >dalacin c penicillin allergy</a>
<a href=http://kathyrnkes.blogg.no>ibux barn dosering</a>
pergotime 50 mg pris
<a href= http://devinnair.startbewijs.nl >proscar finasteride msd</a>
<a href=http://lincolnmaw.ek.la>doxylin mot lungebetennelse</a>
ibux barn 10 år

SamuelRexSamuelRex2017/04/01 09:08ibux gel barn
<a href= http://delmydavid.blogg.se >ibux 400 mg dose</a>
<a href=http://emmanuelje.tribalpages.com>proscar (finasteride) and avodart (dutasteride)</a>
pergotime bivirkninger
<a href= http://ellsworthut.page4.me >dalacin</a>
<a href=http://trinidadpi.bloggsite.se>zovirax 5 cream for genital herpes</a>
proscar apoteket

ReubencoarpReubencoarp2017/04/01 14:03comprar qlaira preço
<a href= http://www.onweb.fr?company_id=17318 >qlaira engorda muito</a>
<a href=http://jonsuchars.oneminutesite.it>ciprofloxacino posologia para que serve</a>
ciprofloxacino colirio preço
<a href= http://andreahatt.wifeo.com >excedrin pm ingredients</a>
<a href=http://janethorua.guildomatic.com>augmentin</a>
ciprofloxacino preço nissei

WilliamdenWilliamden2017/04/03 09:12qlaira bula efeitos colaterais
<a href= http://jorgemanwa.jimdo.com >diclofenaco pomada preço</a>
<a href=http://lauriceeag.freeblog.site>ciprofloxacino hidrocortisona otologico</a>
nexium 20 mg preço ultrafarma
<a href= http://lauriceeag.freeblog.site >ciprofloxacino preço drogasil</a>
<a href=http://chucklyshe.edublogs.org/2017/04/02/nexium/>remedio nexium 20 mg preço</a>
nexium 20 mg bula pdf

MichaeldowMichaeldow2017/04/06 07:34finasteride 1 mg precio
<a href= http://elwoodeast.soup.io >nolvadex 20 mg effets secondaires</a>
<a href=http://project200559.tilda.ws>amoxicilline acide clavulanique grippe</a>
avis clomid duphaston
<a href= http://leroypusey.startspot.nl >tolexine 50 prix</a>
<a href=http://earledelli.e-monsite.com/pages/clomid.html>clomid et ovitrelle opk</a>
amoxicilline 250 mg posologie

MatthewfubMatthewfub2017/04/06 09:51dapoxetine 30 mg price in india
<a href= http://salomemahe.page4.me >prix finasteride 1mg</a>
<a href=http://jolenemale.onlc.be>nolvadex dosage</a>
nolvadex 20 mg effets secondaires
<a href= http://salomemahe.page4.me >finasteride 5mg prix</a>
<a href=http://margorieza.spruz.com>doxyval 100 prix</a>
doxyval chien prix

RichardQuavaRichardQuava2017/04/10 09:36etoricoxib torrinomedica
<a href= http://cruzliedy.guildomatic.com >bijsluiter celecoxib sandoz</a>
<a href=http://markettaro-blog.logdown.com>amoxicillina sciroppo torrino</a>
amoxicillina sciroppo torrino
<a href= http://elishakoew.ayosport.com/article-15572-zovirax >zovirax pomata per herpes zoster</a>
<a href=http://jinnychoin.familybelle.com>janumet 50/1000 mg</a>
janumet 50 mg/500

GuestspuckGuestspuck2018/03/11 22:49guest test post
<a href=" http://temresults2018.com/ ">bbcode</a>
<a href="http://temresults2018.com/">html</a>
http://temresults2018.com/ simple