project euler 66を解いていて力任せなコードでunsigned long型のxについてx2を計算したらオーバーフローしてしまった。 オーバーフロー検出としては以前利用したlimits.hをインクルードして、 #include<limits.h> if (x > ULONG_MAX / x) return 1; とするか、gccなら</limits.h>…

C++の練習がてらProject Eulerの続きを解いている。66番はx2-Dy2=1を満たす自然数x,yが存在するD(ただし1<=D<=1000、平方数でない)のうちxが最大になるものを求めるというもの。 ペル方程式についてはペル方程式とチェビシェフ多項式（その１）に性質の解説…

簡易DICOM Viewerをつくろうを参考にしつつ、DICOM sample image setsのデータを試しに開いてみることにする。