青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

pku 1266 ural 1043 Cover an Arc. 解析幾何,注意浮點數下取整要用floor函數

題意:
給出一段圓弧的起點、終點以及第三點(三點不共線),求4個點在整點的最小矩形將其圍住。

測試數據:
/Files/yzhw/cover.rar

解法:
首先,肯定要求圓心的
大家數學都很牛,我就給一個圖,一個式子,你懂的

CE2=DB2+DE2
然后用解析法做吧- -,要討論下A、B是否為水平線或者豎直線

下面就是確定上下左右整點坐標的問題了
還記得math庫里有atan2?嘻嘻,那就好辦了
要分兩種情況討論,

A、B是弧的端點。
也許這兩種情況從幾何上看MS一樣的,但是對于求atan2就不同了
大家知道,atan2返回值是-PI~PI
也就是說,(-1,0)向量返回PI,(1,0)返回0,然后(-1,0)向逆時針方向轉一點點就是返回接近-PI的數了,-PI和PI是同一點,但atan2處理的時候PI位置是閉區(qū)間,而-PI位置是開區(qū)間。
廢話了一堆,大家應該明白了,上圖第一種情況C的atan2值是大于a、b的最大值或者小于a、b的最小值的,而第二種情況d的atan2值是介于e、f之間的。
分類討論,然后分別測試-PI向量、PI/2向量、0向量、-PI/2向量是否在圓弧內就可以了。
一個值得注意的地方,下取整應該使用floor函數,應為如果直接int取整,當浮點值小于0的時候就變成上取整了- -
我感覺自己的代碼寫的很到位的,再有不懂得大家看我代碼吧,不過話說java的效率真蛋疼。。一個O(1)的算法竟然能跑5秒。。。

代碼:
  1Source Code
  2
  3Problem: 1266  User: yzhw 
  4Memory: 3024K  Time: 5422MS 
  5Language: Java  Result: Accepted 
  6
  7Source Code 
  8import java.util.*;
  9public class Main {
 10
 11    /**
 12     * @param args
 13     */

 14    static int x1,x2,x3,y1,y2,y3,left,right,up,down;
 15    static double x,y,x4,y4,d,r;
 16    static double dis(double x1,double y1,double x2,double y2)
 17    {
 18        return Math.sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
 19    }

 20    static double cross(double x1,double y1,double x2,double y2)
 21    {
 22        return x1*y2-x2*y1;
 23    }

 24    static boolean gt(double num,double pos)
 25    {
 26        return Math.abs(num-pos)<1e-6||num>pos;
 27    }

 28    static boolean le(double num,double pos)
 29    {
 30        return Math.abs(num-pos)<1e-6||num<pos;
 31    }

 32    static int floor(double num)
 33    {
 34        num+=1e-8;
 35        int res=(int)Math.floor(num);
 36        if(Math.abs(res-num)<1e-6return res;
 37        else return res+1;
 38    }

 39    public static void main(String[] args) {
 40        Scanner in=new Scanner (System.in);
 41        x1=in.nextInt();
 42        y1=in.nextInt();
 43        x2=in.nextInt();
 44        y2=in.nextInt();
 45        x3=in.nextInt();
 46        y3=in.nextInt();
 47        x4=(x1+x2)*0.5;
 48        y4=(y1+y2)*0.5;
 49        d=dis(x1,y1,x4,y4);
 50        if(x1==x2)
 51        {
 52            y=(y1+y2)*0.5;
 53            x=(d*d+x4*x4+y4*y4-x3*x3-y3*y3-y*2.0*(y4-y3))*0.5/(x4-x3);
 54        }

 55        else if(y1==y2)
 56        {
 57            x=(x1+x2)*0.5;
 58            y=(d*d+x4*x4+y4*y4-x3*x3-y3*y3-x*2.0*(x4-x3))*0.5/(y4-y3);
 59        }

 60        else
 61        {
 62            double k=(x1-x2)/(double)(y2-y1);
 63            double t=d*d+x4*x4+y4*y4-x3*x3-y3*y3-2.0*(y4-y3)*(y4-k*x4);
 64            x=t/(2.0*(x4-x3)+2.0*(y4-y3)*k);
 65            y=k*x+y4-k*x4;
 66        }

 67        r=dis(x3,y3,x,y);
 68        //System.out.println(x+" "+y+" "+r);
 69        double start=Math.atan2(y1-y, x1-x),end=Math.atan2(y2-y, x2-x),nxt=Math.atan2(y3-y,x3-x);
 70        if(le(nxt,Math.min(start,end))||gt(nxt,Math.max(start,end)))
 71        {
 72            if(le(Math.PI,Math.min(start,end))||gt(Math.PI,Math.max(start,end)))
 73               left=(int)Math.floor(x-r+1e-8);
 74            else
 75                left=Math.min(x1, x2);
 76            if(le(0,Math.min(start,end))||gt(0,Math.max(start,end)))
 77                   right=floor(x+r);
 78                else
 79                   right=Math.max(x1, x2);
 80            if(le(-Math.PI*0.5,Math.min(start,end))||gt(-Math.PI*0.5,Math.max(start,end)))
 81                   down=(int)Math.floor(y-r+1e-8);
 82                else
 83                   down=Math.min(y1, y2);
 84            if(le(Math.PI*0.5,Math.min(start,end))||gt(Math.PI*0.5,Math.max(start,end)))
 85                   up=floor(y+r);
 86                else
 87                   up=Math.max(y1, y2);
 88        }

 89        else
 90        {
 91            if(le(Math.PI,Math.max(start,end))&&gt(Math.PI,Math.min(start,end)))
 92                   left=(int)Math.floor(x-r+1e-8);
 93                else
 94                    left=Math.min(x1, x2);
 95            if(le(0,Math.max(start,end))&&gt(0,Math.min(start,end)))
 96                       right=floor(x+r);
 97                    else
 98                       right=Math.max(x1, x2);
 99            if(le(-Math.PI*0.5,Math.max(start,end))&&gt(-Math.PI*0.5,Math.min(start,end)))
100                       down=(int)Math.floor(y-r+1e-8);
101                    else
102                       down=Math.min(y1, y2);
103            if(le(Math.PI*0.5,Math.max(start,end))&&gt(Math.PI*0.5,Math.min(start,end)))
104                       up=floor(y+r);
105                    else
106                       up=Math.max(y1, y2);
107        }

108        System.out.println((right-left)*(up-down));
109    }

110
111}

112
113

posted on 2011-01-18 10:44 yzhw 閱讀(634) 評論(0)  編輯 收藏 引用 所屬分類: geometry&phycise

<2011年1月>
2627282930311
2345678
9101112131415
16171819202122
23242526272829
303112345

導航

統(tǒng)計

公告

統(tǒng)計系統(tǒng)

留言簿(1)

隨筆分類(227)

文章分類(2)

OJ

最新隨筆

搜索

積分與排名

最新評論

閱讀排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲欧美区自拍先锋| 国产精品高潮粉嫩av| 亚洲区一区二| 玖玖精品视频| 国产视频一区在线| 亚洲免费在线精品一区| 欧美在线观看一区二区三区| 国产免费观看久久黄| 久久久久国产精品一区二区| 免费一级欧美片在线观看| 亚洲精品欧美日韩专区| 国产精品v欧美精品v日韩精品| 亚洲一区二区视频| 欧美国产免费| 欧美亚洲一级片| 亚洲第一页在线| 国产精品久久久久久久第一福利| 中国成人亚色综合网站| 免费成人黄色av| 亚洲天堂av在线免费| 国产婷婷色一区二区三区四区 | 亚洲国产一成人久久精品| 欧美精品性视频| 久久国产精品一区二区| 日韩一区二区免费高清| 久久人人97超碰国产公开结果| 一区二区三区精品视频| 在线欧美三区| 国内自拍一区| 国产人妖伪娘一区91| 欧美日韩午夜视频在线观看| 久久精品国内一区二区三区| 亚洲精品国产精品国自产观看| 亚洲尤物在线视频观看| 影音先锋成人资源站| 国产精品日韩欧美一区二区| 欧美国产精品| 久久综合给合久久狠狠色| 欧美中文字幕在线| 午夜精品久久久久久久99热浪潮| 亚洲美女av电影| 亚洲激情欧美| 欧美激情一区三区| 久久综合五月天婷婷伊人| 亚洲欧美日韩在线| 亚洲一区三区视频在线观看| 亚洲毛片一区| 亚洲国产精品悠悠久久琪琪| 一区二区三区在线高清| 精品动漫一区二区| 在线精品国精品国产尤物884a| 国产一区二区毛片| 国产亚洲欧美日韩美女| 国产精品免费一区二区三区观看| 91久久精品久久国产性色也91| 美女主播精品视频一二三四| 久久久国产精品亚洲一区| 欧美亚洲日本一区| 欧美在线看片| 久久久精品国产99久久精品芒果| 欧美在线视频不卡| 久久婷婷蜜乳一本欲蜜臀| 噜噜噜91成人网| 美女视频黄 久久| 亚洲福利精品| 亚洲精选中文字幕| 中文在线不卡视频| 亚洲免费中文| 久久精品国产v日韩v亚洲| 久久天天躁狠狠躁夜夜爽蜜月| 牛夜精品久久久久久久99黑人 | 一区二区激情视频| 亚洲综合不卡| 久久九九免费视频| 欧美国产日本韩| 一区二区日韩免费看| 亚洲综合国产| 牛牛精品成人免费视频| 欧美日韩一区二区三区在线视频| 国产精品美腿一区在线看| 国产又爽又黄的激情精品视频 | 久久全国免费视频| 欧美另类在线观看| 国产精品一区三区| 在线精品视频在线观看高清| 日韩天天综合| 久久精品女人天堂| 亚洲激情一区二区| 午夜精品久久久久久久蜜桃app | 欧美成人久久| 一区二区三区欧美激情| 久久精品72免费观看| 欧美日韩成人在线观看| 黄色成人av在线| 亚洲视频一区| 欧美激情第三页| 午夜视频在线观看一区| 欧美国产成人精品| 国产一区二区三区日韩欧美| 99视频有精品| 欧美成人精品激情在线观看 | 亚洲视频精品在线| 91久久精品日日躁夜夜躁国产| 亚洲一区二区三区四区五区午夜 | 欧美视频一区二区三区| 在线观看一区| 久久久久久一区二区| 一本色道久久88综合日韩精品| 久久综合99re88久久爱| 国产日韩欧美亚洲| 午夜精品视频一区| 亚洲精品视频在线看| 久久久综合网| 黄色资源网久久资源365| 久久精品91| 午夜精品福利一区二区蜜股av| 欧美日韩国产成人| 亚洲精品四区| 亚洲成人中文| 一本久道久久综合婷婷鲸鱼| 久久亚洲视频| 一本色道久久综合狠狠躁篇的优点| 亚洲欧美99| 久久在线免费观看视频| 一本色道88久久加勒比精品 | 久久精品成人一区二区三区| 国产精品五月天| 亚洲午夜精品久久久久久浪潮 | 欧美国产精品劲爆| 麻豆精品网站| 在线看成人片| 亚洲国产你懂的| 国产日韩免费| 午夜在线观看欧美| 亚洲一区高清| 国产精品视频网站| 欧美一区二区三区免费在线看| 99精品视频免费观看视频| 欧美日本精品在线| 亚洲婷婷综合久久一本伊一区| 亚洲人成网站影音先锋播放| 欧美成人精品在线视频| 日韩视频―中文字幕| 一区二区三区日韩在线观看| 国产精品视频一区二区三区| 久久精品首页| 久久亚洲二区| 一本大道av伊人久久综合| 日韩亚洲不卡在线| 国产欧美亚洲精品| 国产精品福利网站| 欧美一区二区性| 久久久久一区| 日韩午夜剧场| 午夜在线观看欧美| 亚洲国产成人在线| 亚洲精品自在久久| 国产亚洲精品资源在线26u| 久久久精品999| 欧美日韩国产一区二区| 欧美一区二区视频在线观看| 久久久91精品国产一区二区三区| 亚洲精品美女免费| 亚洲男同1069视频| 亚洲国产精品嫩草影院| 亚洲毛片一区| 伊人夜夜躁av伊人久久| 99re66热这里只有精品4| 国产亚洲欧美日韩精品| 亚洲人成网站影音先锋播放| 国产小视频国产精品| 91久久久在线| 黄色在线成人| 亚洲欧美日本伦理| 99精品久久久| 久久人人看视频| 欧美一区二区三区日韩视频| 欧美国产日韩亚洲一区| 久久久久久久精| 国产精品日韩在线| 99精品国产热久久91蜜凸| 亚洲国产日韩一级| 欧美一区二区三区四区在线观看地址| 亚洲精品一区二区三区不| 午夜国产精品视频| 亚洲一区二区三区四区五区黄| 免费成人高清视频| 久久综合电影| 国产日韩亚洲欧美| 亚洲午夜极品| 亚洲视屏在线播放| 欧美区一区二区三区| 欧美黄网免费在线观看| 国产亚洲综合性久久久影院| 这里只有精品丝袜| 亚洲欧美成人| 国产精品国产三级国产aⅴ9色| 亚洲丰满少妇videoshd| 亚洲福利免费| 麻豆精品视频在线观看视频|