Codeforces Round #176 div 2
久しぶりに(一年ぐらい)Codeforcesに出ました。
正解したのはA問題のみというひどい結果でした。
B問題はコーナーケースで落ち、C問題以降は
手も足も出ませんでした。
もう一度初心に帰って勉強する必要性を感じました。
A. IQ Test
問題概要
4x4の升目が与えられます。各升には"."または"#"が書かれています。
16個の升のどれか一つを書き換えることによって同じ記号が書かれた
2x2の正方形を得ることができるなら"YES"、できないなら"NO"と出力しなさい。
ただし、何も書き換えない状態で条件を満たすなら"YES"と出力しなさい。
解法
16個の升のそれぞれについて書き換えて正方形が得られるかどうかを調べるだけです。
ソースコード
#include <string> #include <vector> #include <algorithm> #include <iostream> using namespace std; bool check(vector<string> &vs){ for(int i=0;i<3;i++) for(int j=0;j<3;j++) if(vs[i][j]==vs[i+1][j]&&vs[i][j]==vs[i][j+1]&&vs[i][j]==vs[i+1][j+1]) return true; return false; } char opposite(char c){ if(c=='#')return '.'; else return '#'; } int main(){ vector<string> vs; for(int i=0;i<4;i++){ string s; cin>>s; vs.push_back(s); } string ans="NO"; if(check(vs)) ans="YES"; for(int i=0;i<4;i++) for(int j=0;j<4;j++){ vs[i][j]=opposite(vs[i][j]); if(check(vs)) ans="YES"; vs[i][j]=opposite(vs[i][j]); } cout<<ans<<endl; return 0; }
B. Pipeline
問題概要
自然数nとkが与えられます。
k以下の互いに異なる自然数l_1,l_2,l_3,...,l_mが
l_1-1+l_2-1+l_3-1+...+l_m-1=n-1
を満たすときmの最小値を求めなさい。
ただし、n=1のときはm=0です。
(これは問題から主要な部分を抽出したものであり
実際の問題文の日本語訳とはかなり異なります。)
解法
l_1-1+l_2-1+l_3-1+...+l_m-1>=n-1
となるmを見つければ良いので(勘)
l_1=k,l_2=k-1,l-3=k-2,...,l-m=k-m+1とします。
このとき
l_1-1+l_2-1+l_3-1+...+l_m-1=(k+k-m+1)*(k-k+m-1+1)/2-m であるので
mの値を二分探索すれば良いです。
ソースコード
#include <iostream> #include <algorithm> using namespace std; typedef long long LL; int main(){ LL n,k; cin>>n>>k; if(n==1){ cout<<0<<endl; return 0; } if(n<=k){ cout<<1<<endl; return 0; } LL u=k-1,d=0; for(int i=0;i<100;i++){ LL m=(u+d)/2; LL a=(k-m)*(k+m-1)/2; if(n-1<=a) d=m; else u=m; } // cout<<"u="<<u<<" d="<<d<<endl; if(n-1<=(k-u)*(k+u-1)/2) cout<<k-u<<endl; else if(n-1<=(k-d)*(k+d-1)/2) cout<<k-d<<endl; else cout<<-1<<endl; return 0; }
Marathon Match 78 の記録
最終的なアルゴリズム
左上(0,0)を出発点として、盤面上に一本の線を引いていく。
1マス進むことを一手として10手探索してはもっとも高い点数が
得られるものを最終的な解に追加していく。このとき輪っかができなくなるものは除外する。
探索と解の追加を繰り返して出発点に戻ってきたら終了する。
経過
1:左上を始点として点数が最大になる輪っかを作るようなビーム探索を書く。
(当時は自分が書いているものがビーム探索であるとは知らなかった。)
※ビーム探索とは幅優先探索において評価の高い解のみを残すような探索のこと。
2:うーん、遅すぎる。
3:だんだん保持する解の個数を減らす。
4:あーもうこれは保持する解を一つにしてもいいんじゃない。
5:保持する解を一つにする。
6:10手づつ深さ優先探索して暫定解に付け足していって輪っかができたら終了するアルゴリズムに変更する。
結果
105人中57位。
初めてのマラソンマッチにしては健闘したと見るか、
修行が足りないと見るかは判断が分かれるところ。(僕の頭の中では)
反省点
大量のテストを簡単に実行できる環境を早めに構築するべきだった。
(最終的にはコマンド一つで100個のテストケースを試すようにした。)
更に、テストから得られた結果を簡単に判断できるようにするべきだった。
良かった点
公式のビジュアライザを改造して、サイズやfillRatioをコマンド引数で変更できるようにし、
更に正しくない(輪っかになっていない)解も目で見れるようにしたのは良かった。
最終提出コード
#include <iostream> #include <cstdio> #include <cstdlib> #include <cmath> #include <string> #include <vector> #include <queue> #include <map> #include <set> #include <algorithm> #include <functional> #include <cstring> #include <sstream> #define FOR(i,k,n) for (int i=(k); i<(int)(n); ++i) #define REP(i,n) FOR(i,0,n) #define FORIT(i,c) for(__typeof((c).begin())i=(c).begin();i!=(c).end();++i) #define SZ(i) ((int)i.size()) #define pb push_back #define mp make_pair #define ALL(X) (X).begin(),(X).end() typedef long long LL; // #define DFS_DEPTH 11 #define MAX_ANS_LENGTH 20200 #define INF 100000000 using namespace std; int dfsDepth=10; int mv[4][2]={{-1,0},{1,0},{0,-1},{0,1}}; int c2m(char c){if(c=='U')return 0;if(c=='D')return 1;if(c=='L')return 2;if(c=='R')return 3;return INF;} int iabs(int n){if(0<=n)return n;else return -n;} vector<string> dia; int H,W; int isVisited[110][110]; int nVisit[110][110]; int nTouch[110][110]; int bfs[110][110]; string udlr="UDLR"; string snake; int point; pair<pair<int,int>,string> ans; int ystart,xstart,yhead,xhead,ygoal,xgoal; class BIT{ public: int bit[MAX_ANS_LENGTH+10]; int sum(int i){ i++; int s=0; while(0<i){ s+=bit[i]; i-=i&-i; } return s; } void add(int i,int x){ i++; while(i<=MAX_ANS_LENGTH+5){ bit[i]+=x; i+=i&-i; } } }; BIT turnCount; int rorlturn(char s,char t){//left turn = 1 right turn = -1 if(s=='U'){if(t=='L')return 1;if(t=='R')return -1;} if(s=='L'){if(t=='D')return 1;if(t=='U')return -1;} if(s=='D'){if(t=='R')return 1;if(t=='L')return -1;} if(s=='R'){if(t=='U')return 1;if(t=='D')return -1;} return 0; } bool isInLoop(const int w,const int v,const char n,const int y,const int x,const char m,const int nv){ if(y==0||y==H||x==0||x==W)return false; //右回りのループを検出 int yy=y,xx=x; if(m=='U')xx--; if(m=='D')xx++; if(m=='L')yy++; if(m=='R')yy--; if(0<isVisited[yy][xx]&&turnCount.sum(nv)-turnCount.sum(nVisit[yy][xx]+1)<=-3) return true; //左回りのループを検出 yy=y,xx=x; if(m=='U')xx++; if(m=='D')xx--; if(m=='L')yy--; if(m=='R')yy++; if(0<isVisited[yy][xx]&&turnCount.sum(nv)-turnCount.sum(nVisit[yy][xx]+1)>=3) return true; int yf=w+mv[c2m(n)][0],xf=v+mv[c2m(n)][1]; if(yf==0||yf==H||xf==0||xf==W)return false; //右回りのループを検出 yy=yf,xx=xf; if(n=='U')xx--; if(n=='D')xx++; if(n=='L')yy++; if(n=='R')yy--; if(0<isVisited[yy][xx]&&turnCount.sum(nv-1)-turnCount.sum(nVisit[yy][xx]+1)<=-3&&rorlturn(n,m)==-1) return true; //左回りのループを検出 yy=yf,xx=xf; if(n=='U')xx++; if(n=='D')xx--; if(n=='L')yy--; if(n=='R')yy++; if(0<isVisited[yy][xx]&&turnCount.sum(nv-1)-turnCount.sum(nVisit[yy][xx]+1)>=3&&rorlturn(n,m)==1) return true; return false; } pair<pair<int,int>,string> DFS(const int y,const int x,const string s,const int p,const char l){ if(dfsDepth==SZ(s)){ return make_pair(make_pair(p,iabs(ygoal-y)+iabs(xgoal-x)),s); }else{ vector<pair<pair<int,int>,string> > res; for(int i=0;i<4;i++){ int yy=y+mv[i][0]; int xx=x+mv[i][1]; if(yy<0||H<yy||xx<0||W<xx) continue; if(ygoal==H&&xgoal==W&&yy<xx) continue; if(ygoal==H&&xgoal==W&&yy==xx&&i==0) continue; if(ygoal==H&&xgoal==W&&y==x+1&&i==2) continue; vector<int> tpos; string ss=s+udlr[i]; char ll=udlr[i]; if(i==0){ if(0<x) tpos.push_back((y-1)*1000+x-1); if(x<W) tpos.push_back((y-1)*1000+x); }else if(i==1){ if(0<x) tpos.push_back(y*1000+x-1); if(x<W) tpos.push_back(y*1000+x); }else if(i==2){ if(0<y) tpos.push_back((y-1)*1000+x-1); if(y<H) tpos.push_back(y*1000+x-1); }else{ if(0<y) tpos.push_back((y-1)*1000+x); if(y<H) tpos.push_back(y*1000+x); } int pp=p; for(int j=0;j<SZ(tpos);j++){ int yyy=tpos[j]/1000,xxx=tpos[j]%1000; if(nTouch[yyy][xxx]==dia[yyy][xxx]-'0') pp--; if(nTouch[yyy][xxx]+1==dia[yyy][xxx]-'0') pp++; nTouch[yyy][xxx]++; } isVisited[yy][xx]++; int oldnVisit=nVisit[yy][xx]; nVisit[yy][xx]=SZ(snake)+SZ(ss); turnCount.add(SZ(snake)+SZ(ss),rorlturn(l,ll)); if((x==0&&i==0)||(y==H&&i==2)||(x==W&&i==1)||(y==0&&i==3)){ ; }else if(ygoal==yy&&xgoal==xx){ if(2<snake.size()+ss.size()&&ans.first.first<point+pp) ans=make_pair(make_pair(point+pp,iabs(ygoal-yy)+iabs(xgoal-xx)),snake+ss); }else if(1<isVisited[yy][xx]){ ; }else{ if(!isInLoop(y,x,l,yy,xx,udlr[i],SZ(snake)+SZ(ss))){ res.push_back(DFS(yy,xx,ss,pp,ll)); } } for(int j=0;j<SZ(tpos);j++){ int yyy=tpos[j]/1000,xxx=tpos[j]%1000; nTouch[yyy][xxx]--; } isVisited[yy][xx]--; nVisit[yy][xx]=oldnVisit; turnCount.add(SZ(snake)+SZ(ss),-1*rorlturn(l,ll)); } if(SZ(res)==0) return make_pair(make_pair(-INF,iabs(ygoal-y)+iabs(xgoal-x)),""); else{ sort(res.begin(),res.end()); return res.back(); } } } string solve(void){ memset(turnCount.bit,0,sizeof(turnCount)); turnCount.add(100,1); if(H<=20){ dfsDepth=12; }else if(H<=30){ dfsDepth=11; }else if(H<=40){ dfsDepth=10; }else if(H<=50){ dfsDepth=10; }else if(H<=60){ dfsDepth=9; }else if(H<=70){ dfsDepth=9; }else if(H<=80){ dfsDepth=9; }else if(H<=90){ dfsDepth=8; }else if(H<=100){ dfsDepth=8; } point=0; FOR(i,0,H)FOR(j,0,W)if(dia[i][j]=='0')point++; memset(isVisited,0,sizeof(isVisited)); memset(nVisit,-1,sizeof(nVisit)); memset(nTouch,0,sizeof(nTouch)); ystart=0;xstart=0;isVisited[0][0]=1;nVisit[0][0]=0; yhead=1;xhead=0;isVisited[1][0]=1;nVisit[1][0]=1; ygoal=H;xgoal=W; snake="D"; ans=make_pair(make_pair(-INF,iabs(ygoal-yhead)+iabs(xgoal-xhead)),""); nTouch[0][0]=1; if(dia[0][0]=='0')point--; if(dia[0][0]=='1')point++; for(int time=0;time<MAX_ANS_LENGTH/dfsDepth/2;time++){ pair<pair<int,int>,string> res=DFS(yhead,xhead,"",0,snake[SZ(snake)-1]); if(res.first.first==-INF)break; point+=res.first.first; snake+=res.second; FOR(i,0,SZ(res.second)){ char ch=res.second[i]; if(ch=='U')yhead--; if(ch=='D')yhead++; if(ch=='L')xhead--; if(ch=='R')xhead++; isVisited[yhead][xhead]++; nVisit[yhead][xhead]=SZ(snake)-SZ(res.second)+i+1; turnCount.add(SZ(snake)-SZ(res.second)+i,rorlturn(snake[SZ(snake)-SZ(res.second)+i-1],snake[SZ(snake)-SZ(res.second)+i])); } } fprintf(stderr,"ans.first.first=%d ans.second=\"%s\"\n",ans.first.first,ans.second.c_str()); if(ans.first.first==-INF){ stringstream tss; tss<<0<<" "<<0<<" "<<snake<<endl; return tss.str(); } pair<pair<int,int>,string> fans=ans; memset(isVisited,0,sizeof(isVisited)); memset(nTouch,0,sizeof(nTouch)); memset(nVisit,-1,sizeof(nVisit)); memset(turnCount.bit,0,sizeof(turnCount.bit)); yhead=0,xhead=0; isVisited[yhead][xhead]++; nVisit[yhead][xhead]=0; FOR(i,0,SZ(fans.second)){ char ch=fans.second[i]; if(ch=='U')yhead--; if(ch=='D')yhead++; if(ch=='L')xhead--; if(ch=='R')xhead++; isVisited[yhead][xhead]++; nVisit[yhead][xhead]=i+1; if(0<i) turnCount.add(i,rorlturn(fans.second[i-1],fans.second[i])); } point=0; for(int y=0;y<H;y++) for(int x=0;x<W;x++){ nTouch[y][x]+=(isVisited[y][x]&&isVisited[y+1][x]); nTouch[y][x]+=(isVisited[y+1][x]&&isVisited[y+1][x+1]); nTouch[y][x]+=(isVisited[y+1][x+1]&&isVisited[y][x+1]); nTouch[y][x]+=(isVisited[y][x+1]&&isVisited[y+1][x]); point+=(nTouch[y][x]==dia[y][x]-'0'); } snake=fans.second+"U"; ystart=H;xstart=W;isVisited[H][W]=1; turnCount.add(SZ(snake)-1,rorlturn('R','U')); yhead=H-1;xhead=W;isVisited[H-1][W]=1;nVisit[H-1][W]=snake.size(); ygoal=0;xgoal=0; ans=make_pair(make_pair(-INF,iabs(ygoal-yhead)+iabs(xgoal-xhead)),""); for(int time=0;time<MAX_ANS_LENGTH/dfsDepth/2;time++){ pair<pair<int,int>,string> res=DFS(yhead,xhead,"",0,snake[SZ(snake)-1]); if(res.first.first==-INF)break; point+=res.first.first; snake+=res.second; FOR(i,0,SZ(res.second)){ char ch=res.second[i]; if(ch=='U')yhead--; if(ch=='D')yhead++; if(ch=='L')xhead--; if(ch=='R')xhead++; isVisited[yhead][xhead]++; nVisit[yhead][xhead]=SZ(snake)-SZ(res.second)+i+1; turnCount.add(SZ(snake)-SZ(res.second)+i,rorlturn(snake[SZ(snake)-SZ(res.second)+i-1],snake[SZ(snake)-SZ(res.second)+i])); } } fprintf(stderr,"ans.first.first=%d ans.second=\"%s\"\n",ans.first.first,ans.second.c_str()); fprintf(stderr,"snake=\"%s\"\n",snake.c_str()); if(ans.first.first==-INF){ stringstream tss; tss<<0<<" "<<0<<" "<<snake<<endl; return tss.str(); } stringstream ss; ss<<0<<" "<<0<<" "<<ans.second; return ss.str(); } class FixTheFence{ public: string findLoop(const vector<string>& __diagram){ dia=__diagram; H=SZ(dia); W=SZ(dia[0]); return solve(); } };
HaskellとOpenGLを使って四角形を描画しテクスチャを貼り付ける
こんな感じになります↓
結構苦労しました。誰かの役に立ってくれると嬉しいです。
参考サイト
いろいろなサンプルがある:Index of /GLUT/examples/RedBook
今回役に立ったのは Tex〜.hs
視野の設定:Lichu's_Base
テクスチャの読み込み・設定:https://github.com/fumieval/free-game/blob/master/Graphics/FreeGame/Backends/GLFW.hs
import Graphics.UI.GLUT hiding (Bitmap) import qualified Graphics.Rendering.OpenGL.GL as GL import System.Exit ( exitFailure, exitWith, ExitCode(ExitSuccess) ) import Data.IORef import Codec.Picture.Repa import Data.Array.Repa as R hiding (reshape) import qualified Data.Array.Repa.Repr.ForeignPtr as RF import Foreign.ForeignPtr import Data.Word import Control.Applicative import Control.Monad import Unsafe.Coerce -- テクスチャをファイルから読み込む -- 何をやっているのか未だによくわからない -- FreeGameのサンプルをちょっといじったもの loadTextureFromFile :: FilePath -> IO GL.TextureObject loadTextureFromFile path = do content <- delay <$> (flipVertically.imgData) <$> either error id <$> (readImageRGBA path) let (Z :. width :. height :. _) = R.extent content [tex] <- GL.genObjectNames 1 GL.textureBinding GL.Texture2D GL.$= Just tex GL.textureFilter Texture2D $= ((Nearest, Nothing), Nearest) fptr <- liftM RF.toForeignPtr $ R.computeP $ content withForeignPtr fptr $ GL.texImage2D Nothing GL.NoProxy 0 GL.RGBA8 (GL.TextureSize2D (gsizei width) (gsizei height)) 0 . GL.PixelData GL.RGBA GL.UnsignedInt8888 return tex gsizei :: Int -> GL.GLsizei {-# INLINE gsizei #-} gsizei x = unsafeCoerce x --タイマの間隔 timerInterval = 40 main = do --回転の角度を初期化 rot <- newIORef 0.0 arg <- newIORef 10.4 --GLUTの初期化 initialDisplayMode $= [RGBAMode, DoubleBuffered] initialWindowSize $= Size 640 480 initialize "" [] --ウィンドウを作る createWindow "DrawPictureWithOpneGL" --テクスチャを読み込む tex <- loadTextureFromFile "読み込みたい画像ファイルへのパス" texture Texture2D $= Enabled GL.blend $= GL.Enabled GL.blendFunc $= (GL.SrcAlpha, GL.OneMinusSrcAlpha) --表示に使うコールバック関数の指定 displayCallback $= display rot arg tex --ウィンドウのサイズが変更された時に呼ぶコールバック関数の指定 reshapeCallback $= Just reshape --キーボードやマウスのコールバック keyboardMouseCallback $= Just (keyboardProc arg) --タイマを作る addTimerCallback timerInterval $ timerProc (display rot arg tex) --GLUTのメインループに入る mainLoop display rot arg texname = do --回転させる w<-readIORef arg --w <- readIORef hoge modifyIORef rot (+w) r <- readIORef rot --背景を黒にする clearColor $= Color4 0.0 0.0 0.0 0.0 clear [ColorBuffer] --単位行列を読み込む loadIdentity -- texCoord,vetexは多くの型について定義されているので -- オーバーロードに関する問題を避けるために以下の関数が必要 let texCoord2f = texCoord :: TexCoord2 GLfloat -> IO () vertex3f = vertex :: Vertex3 GLfloat -> IO () -- 表示 -- 描写する直前に色を指定するようだ -- Color4 赤 緑 青 透明度(0で完全に透明、1で完全に不透明) currentColor $= Color4 1 1 1 1 -- テクスチャの設定 textureBinding Texture2D $= Just texname preservingMatrix $ do translate (Vector3 100 100 (0::Double)) -- rotate 度数法での角度 rotate r (Vector3 0.0 0.0 1.0 :: Vector3 GLdouble) -- 描画 renderPrimitive Quads $ do -- ここから texCoord2f (TexCoord2 0 0) vertex3f (Vertex3 (-50) (-50) 0.0) -- ここまでで一つの頂点の設定(テクスチャに関するもの、場所に関するもの) texCoord2f (TexCoord2 0 1) vertex3f (Vertex3 (-50) 50 0.0) texCoord2f (TexCoord2 1 1) vertex3f (Vertex3 50 50 0.0) texCoord2f (TexCoord2 1 0) vertex3f (Vertex3 50 (-50) 0.0) --バッファの入れ替え swapBuffers --タイマが呼ばれるたびにactを繰り返す timerProc act = do act addTimerCallback timerInterval $ timerProc act --ウィンドウのサイズが変更された時の処理 reshape (Size w h)=do viewport $= (Position 0 0, (Size w h)) --ウィンドウ全体を使う --ビューボリュームの設定- matrixMode $= Projection loadIdentity -- 視野のセッティング ortho 0.0 640.0 0.0 480.0 (-1000.0) 1000.0 -- これ大事(理由はよくわからない) matrixMode $= Modelview 0 --キー入力の処理 keyboardProc arg ch state _ _ -- Qが押されたら終了 | ch == Char 'q' = exitWith ExitSuccess -- Aが押されたら回転速度を二倍にする | ch == Char 'a' = modifyIORef arg (*(2)) -- Sが押されたら回転速度を二分の1にする | ch == Char 's' = modifyIORef arg (*(0.5)) -- それ以外なら回転の方向を変える | state == Down = modifyIORef arg (*(-1)) | otherwise = return ()
HaskellとOpenGLをつかって四角形を描画する
うまく動いたら雛形にでも使ってやってください。
参考サイト
GLUTによる「手抜き」OpenGL入門
プログラミング/Haskell/HSDL - Flightless wing
import Graphics.UI.GLUT import Graphics.Rendering.OpenGL.GLU import System.Exit import Data.IORef --タイマの間隔 timerInterval = 40 main = do --回転の角度を初期化 rot <- newIORef 0.0 arg <- newIORef 10.4 --GLUTの初期化 initialDisplayMode $= [RGBAMode, DoubleBuffered] initialWindowSize $= Size 640 480 initialize "" [] --ウィンドウを作る createWindow "GLUTpractice" --表示に使うコールバック関数の指定 displayCallback $= display rot arg --ウィンドウのサイズが変更された時に呼ぶコールバック関数の指定 reshapeCallback $= Just reshape --キーボードやマウスのコールバック keyboardMouseCallback $= Just (keyboardProc arg) --タイマを作る addTimerCallback timerInterval $ timerProc (display rot arg) --GLUTのメインループに入る mainLoop display rot arg= do --回転させる w<-readIORef arg --w <- readIORef hoge modifyIORef rot (+w) r <- readIORef rot --背景を黒にする clearColor $= Color4 0.0 0.0 0.0 0.0 clear [ColorBuffer] --単位行列を読み込む loadIdentity --表示 currentColor $= Color4 0 0 1 0 -- 描写する直前に色を指定するようだ preservingMatrix $ do translate (Vector3 100 100 (0::Double)) rotate r (Vector3 0.0 0.0 1.0 :: Vector3 GLdouble) -- rotate 度数法での角度 {-- renderPrimitive Quads $ mapM_ vertex [ Vertex3 0.10 0.10 0.0, Vertex3 (-0.10) 0.10 0.0, Vertex3 (-0.10) (-0.10) 0.0, Vertex3 0.10 (-0.10) 0.0 :: Vertex3 GLfloat] --} renderPrimitive Quads $ mapM_ vertex [ Vertex3 50 50 0.0, Vertex3 (-50) 50 0.0, Vertex3 (-50) (-50) 0.0, Vertex3 50 (-50) 0.0 :: Vertex3 GLfloat] --バッファの入れ替え swapBuffers --タイマが呼ばれるたびにactを繰り返す timerProc act = do act addTimerCallback timerInterval $ timerProc act --ウィンドウのサイズが変更された時の処理 reshape (Size w h)=do viewport $= (Position 0 0, (Size w h)) --ウィンドウ全体を使う --ビューボリュームの設定- matrixMode $= Projection loadIdentity ortho 0.0 640.0 0.0 480.0 (-1000.0) 1000.0 -- ortho2D 0.0 640.0 0.0 480.0 matrixMode $= Modelview 0 -- これ大事 {- --少し後ろから撮影 lookAt (Vertex3 0.0 0.0 (600.0)) (Vertex3 0.0 0.0 0.0) (Vector3 0.0 1.0 0.0) matrixMode $= Modelview 0 -} --キー入力の処理 keyboardProc arg ch state _ _ | ch == Char 'q' = exitWith ExitSuccess --qが押されたら終了 | ch == Char 'a' = modifyIORef arg (*(2)) | ch == Char 's' = modifyIORef arg (*(0.5)) | state == Down = modifyIORef arg (*(-1)) --それ以外なら回転の方向を変える | otherwise = return ()
cabal install を使う前に
apt-cache search パッケージ名
とするとよいかも
AOJ 0090 Overlaps of Seals
問題の内容
10*10の正方形上に半径1の円をn個(n<=100)配置する。
最も多くの円の内部に含まれる正方形上の点はいくつの円に含まれるか。
(円の内部とは円周上も含む。)
解法
想定解法では任意の二つの円についてその交点を求め
いくつの円に含まれるかを数えるらしい。
AOJ : 0090 - Overlaps of Seals - Respect2Dの日記
AOJ Problem 0090 : Overlaps of Seals - kyuridenamidaのチラ裏
AOJ 0090: Overlaps of Seals:Snowing day:So-net blog
嘘解法
正方形上の点をできるだけ多く列挙して
その一つ一つについていくつの円に含まれるかを調べる。
円に含まれるかどうかの判定を甘くすると通った。
ソース
#include <stdio.h> #include <string.h> #include <algorithm> #include <vector> #include <iostream> #include <map> #include <set> using namespace std; typedef long long LL; double x[100],y[100]; int main(){ int n; while(1){ scanf("%d",&n); if(n==0) break; for(int i=0;i<n;i++) scanf("%lf,%lf\n",&x[i],&y[i]); int ans=0; for(double a=0;a<=10.1;a+=0.01) for(double b=0;b<=10.1;b+=0.01){ int t=0; for(int i=0;i<n;i++) if((a-x[i])*(a-x[i])+(b-y[i])*(b-y[i])<=1.001) t++; ans=max(ans,t); } printf("%d\n",ans); } return 0; }
cabal install monadius したら結構ハマった
結局こんなことをやった
cabal install OpenGLRaw --reinstall
cabal install StateVar --reinstall
cabal install Tensor --reinstall
cabal install ObjectName --reinstall
cabal install OpenGL --reinstall
cabal install GLURaw --reinstall
cabal install array-0.3.0.1 --reinstall
cabal install containers --reinstall
cabal install GLUT --reinstall
cabal install directory-1.1.0.0 --reinstall
cabal install monadius --reinstall
- history からインストールに必要だと思われるコマンドを抜粋した。
- よく分からないのでとりあえず全部 --reinstall フラグをつけているけど要らないと思う。
ハマりポイント
cabal install array-0.3.0.1 --reinstall
cabal install directory-1.1.0.0 --reinstall
- "cabal install array"だけでバージョン指定をしなかったり、"cabal install directory"だけだと最新版を入れようとしてくれるが、入らない。