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

The Fourth Dimension Space

枯葉北風寒,忽然年以殘,念往昔,語默心酸。二十光陰無一物,韶光賤,寐難安; 不畏形影單,道途阻且慢,哪曲折,如渡飛湍。斬浪劈波酬壯志,同把酒,共言歡! -如夢令

#

discovery(探索發現)下載地址

     摘要: 美國探索節目  閱讀全文

posted @ 2009-08-26 11:18 abilitytao 閱讀(1515) | 評論 (0)編輯 收藏

差分約束系統(System Of Difference Constraints)

(本文假設讀者已經有以下知識:最短路徑的基本性質、Bellman-Ford算法。)
    比如有這樣一組不等式:
   
X1 - X2 <= 0
X1 - X5 <= -1
X2 - X5 <= 1
X3 - X1 <= 5
X4 - X1 <= 4
X4 - X3 <= -1
X5 - X3 <= -3
X5 - X4 <= -3
不等式組(1)

    全都是兩個未知數的差小于等于某個常數(大于等于也可以,因為左右乘以-1就可以化成小于等于)。這樣的不等式組就稱作差分約束系統。
    這個不等式組要么無解,要么就有無數組解。因為如果有一組解{X1, X2, ..., Xn}的話,那么對于任何一個常數k,{X1 + k, X2 + k, ..., Xn + k}肯定也是一組解,因為任何兩個數同時加一個數之后,它們的差是不變的,那么這個差分約束系統中的所有不等式都不會被破壞。
   
    差分約束系統的解法利用到了單源最短路徑問題中的三角形不等式。即對于任何一條邊u -> v,都有:

d(v) <= d(u) + w(u, v)

    其中d(u)和d(v)是從源點分別到點u和點v的最短路徑的權值,w(u, v)是邊u -> v的權值。
    顯然以上不等式就是d(v) - d(u) <= w(u, v)。這個形式正好和差分約束系統中的不等式形式相同。于是我們就可以把一個差分約束系統轉化成一張圖,每個未知數Xi對應圖中的一個頂點Vi,把所有不等式都化成圖中的一條邊。對于不等式Xi - Xj <= c,把它化成三角形不等式:Xi <= Xj + c,就可以化成邊Vj -> Vi,權值為c。最后,我們在這張圖上求一次單源最短路徑,這些三角形不等式就會全部都滿足了,因為它是最短路徑問題的基本性質嘛。
    話說回來,所謂單源最短路徑,當然要有一個源點,然后再求這個源點到其他所有點的最短路徑。那么源點在哪呢?我們不妨自已造一個。以上面的不等式組為例,我們就再新加一個未知數X0。然后對原來的每個未知數都對X0隨便加一個不等式(這個不等式當然也要和其它不等式形式相同,即兩個未知數的差小于等于某個常數)。我們索性就全都寫成Xn - X0 <= 0,于是這個差分約束系統中就多出了下列不等式:
   
X1 - X0 <= 0
X2 - X0 <= 0
X3 - X0 <= 0
X4 - X0 <= 0
X5 - X0 <= 0

不等式組(2)

    對于這5個不等式,也在圖中建出相應的邊。最后形成的圖如下:


圖1

    圖中的每一條邊都代表差分約束系統中的一個不等式?,F在以V0為源點,求單源最短路徑。最終得到的V0到Vn的最短路徑長度就是Xn的一個解啦。從圖1中可以看到,這組解是{-5, -3, 0, -1, -4}。當然把每個數都加上10也是一組解:{5, 7, 10, 9, 6}。但是這組解只滿足不等式組(1),也就是原先的差分約束系統;而不滿足不等式組(2),也就是我們后來加上去的那些不等式。當然這是無關緊要的,因為X0本來就是個局外人,是我們后來加上去的,滿不滿足與X0有關的不等式我們并不在乎。
    也有可能出現無解的情況,也就是從源點到某一個頂點不存在最短路徑。也說是圖中存在負權的圈。這一點我就不展開了,請自已參看最短路徑問題的一些基本定理。

    其實,對于圖1來說,它代表的一組解其實是{0, -5, -3, 0, -1, -4},也就是說X0的值也在這組解當中。但是X0的值是無可爭議的,既然是以它作為源點求的最短路徑,那么源點到它的最短路徑長度當然是0了。因此,實際上我們解的這個差分約束系統無形中又存在一個條件:

X0 = 0

    也就是說在不等式組(1)、(2)組成的差分約束系統的前提下,再把其中的一個未知數的值定死。這樣的情況在實際問題中是很常見的。比如一個問題表面上給出了一些不等式,但還隱藏著一些不等式,比如所有未知數都大于等于0或者都不能超過某個上限之類的。比如上面的不等式組(2)就規定了所有未知數都小于等于0。
   
    對于這種有一個未知數定死的差分約束系統,還有一個有趣的性質,那就是通過最短路徑算法求出來的一組解當中,所有未知數都達到最大值。下面我來粗略地證明一下,這個證明過程要結合Bellman-Ford算法的過程來說明。
    假設X0是定死的;X1到Xn在滿足所有約束的情況下可以取到的最大值分別為M1、M2、……、Mn(當然我們不知道它們的值是多少);解出的源點到每個點的最短路徑長度為D1、D2、……、Dn。
    基本的Bellman-Ford算法是一開始初始化D1到Dn都是無窮大。然后檢查所有的邊對應的三角形不等式,一但發現有不滿足三角形不等式的情況,則更新對應的D值。最后求出來的D1到Dn就是源點到每個點的最短路徑長度。
    如果我們一開始初始化D1、D2、……、Dn的值分別為M1、M2、……、Mn,則由于它們全都滿足三角形不等式(我們剛才已經假設M1到Mn是一組合法的解),則Bellman-Ford算法不會再更新任合D值,則最后得出的解就是M1、M2、……、Mn。
    好了,現在知道了,初始值無窮大時,算出來的是D1、D2、……、Dn;初始值比較小的時候算出來的則是M1、M2、……、Mn。大家用的是同樣的算法,同樣的計算過程,總不可能初始值大的算出來的結果反而小吧。所以D1、D2、……、Dn就是M1、M2、……、Mn。
   
    那么如果在一個未知數定死的情況下,要求其它所有未知數的最小值怎么辦?只要反過來求最長路徑就可以了。最長路徑中的三角不等式與最短路徑中相反:

d(v) >= d(u) + w(u, v)
也就是 d(v) - d(u) >= w(u, v)


    所以建圖的時候要先把所有不等式化成大于等于號的。其它各種過程,包括證明為什么解出的是最小值的證法,都完全類似。
    
    用到差分約束系統的題目有ZJU 2770,祝好運。

轉自:http://imlazy.ycool.com/post.1702305.html

posted @ 2009-08-21 23:34 abilitytao 閱讀(2423) | 評論 (6)編輯 收藏

POJ 3322 Bloxorz I(BFS, 改編自同名游戲,很不錯的一題)

     摘要: Bloxorz I Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 2624 Acce...  閱讀全文

posted @ 2009-08-16 22:25 abilitytao 閱讀(1884) | 評論 (4)編輯 收藏

POJ 1631 Bridging signals——最長上升子序列 DP(nlogn)

題目說了很多,其實就是一個最長上升子序列的問題,傳統的n^2算法在這里會超時,所以改用二分的算法,1000MS的時限用掉110MS,差強人意,但不知還有沒有更快的算法,希望各位大牛能指點一二。
代碼如下:(閱讀前建議先閱讀我轉的上一篇關于算法介紹的文章)
#include<iostream>
#include
<cstdio>
using namespace std;


const int N =  40000;

int a[N], C[N], f[N]; // f[i]用于記錄a[0i]的最大長度

int bsearch(const int *C, int size, const int &a) 

{
    
    
int l=0, r=size-1;
    
    
while( l <= r ) 
        
    
{
        
        
int mid = (l+r)/2;
        
        
if( a > C[mid-1&& a <= C[mid] ) return mid; // >&&<= 換為: >= && <
        
        
else if( a < C[mid] ) r = mid-1;
        
        
else l = mid+1;
        
    }

    
}


int LIS(const int *a, const int &n){
    
    
int i, j, size = 1;
    
    C[
0= a[0]; f[0= 1;
    
    
for( i=1; i < n; ++i ){
        
        
if( a[i] <= C[0] ) j = 0;                 // <= 換為: <
        
        
else if( a[i] >C[size-1] ) j = size++;   // > 換為: >= 
        
        
else j = bsearch(C, size, a[i]);
        
        C[j] 
= a[i]; f[i] = j+1;
        
    }

    
    
return size;
    
}




int main()
{
    
int testcase;
    
int n;
    scanf(
"%d",&testcase);
    
int i,j;
    
for(i=1;i<=testcase;i++)
    
{
        scanf(
"%d",&n);
        
for(j=0;j<n;j++)
            scanf(
"%d",&a[j]);
        printf(
"%d\n",LIS(a,n));
    }

    
return 0;
}

posted @ 2009-08-12 19:14 abilitytao 閱讀(1860) | 評論 (1)編輯 收藏

最長不降子序列的nlogn算法 (轉)

所謂“最長不降子序列問題”,就是在一個給定的序列中尋找一個子序列{ai}滿足:

                       a1<=a2<=...<an

這個問題在一般教材上,往往會作為動態規劃的引例。

即使用如下的狀態轉移方程進行計算:

                   F[i]=max{F[j]}+1     aj<=ai

但是它的復雜度是O(n^2)的,對于稍大的規模便無法承受。

那么有沒有改進的方法呢?答案是肯定的。

----------------------------------------------------------------分割線------------------------------------------------------------------------------------------

我們維護一個數組C[i],這里C[i]表示F值為i的最小數。

不難發現   C[1]<=C[2]<=...<=C[n]

因此我們可以利用C[]通過二分查找確定F[j]的值。

----------------------------------------------------------------分割線------------------------------------------------------------------------------------------

實現如下:

const int N = 1001;

int a[N], C[N], f[N]; // f[i]用于記錄a[0...i]的最大長度

int bsearch(const int *C, int size, const int &a)

{

    int l=0, r=size-1;

    while( l <= r )

    {

        int mid = (l+r)/2;

        if( a > C[mid-1] && a <= C[mid] ) return mid; // >&&<= 換為: >= && <

        else if( a < C[mid] ) r = mid-1;

        else l = mid+1;

    }

}

int LIS(const int *a, const int &n){

     int i, j, size = 1;

     C[0] = a[0]; f[0] = 1;

     for( i=1; i < n; ++i ){

          if( a[i] <= C[0] ) j = 0;                 // <= 換為: <

         else if( a[i] >C[size-1] ) j = size++;   // > 換為: >=

         else j = bsearch(C, size, a[i]);

         C[j] = a[i]; f[i] = j+1;

     }

     return size;

}

------------------------------------------------------------------分割線------------------------------------------------------------------------------------------

至此,我們了解了O(nlogn)的算法,它主要是利用了二分查找的方法對樸素的動態規劃進行加速、優化,從而達到理想的效率。



轉自:http://fqq11679.blog.hexun.com/21632261_d.html

posted @ 2009-08-12 18:27 abilitytao 閱讀(421) | 評論 (0)編輯 收藏

POJ 計算幾何入門題目推薦(轉)

計算幾何題的特點與做題要領:
1.大部分不會很難,少部分題目思路很巧妙
2.做計算幾何題目,模板很重要,模板必須高度可靠。
3.要注意代碼的組織,因為計算幾何的題目很容易上兩百行代碼,里面大部分是模板。如果代碼一片混亂,那么會嚴重影響做題正確率。
4.注意精度控制。
5.能用整數的地方盡量用整數,要想到擴大數據的方法(擴大一倍,或擴大sqrt2)。因為整數不用考慮浮點誤差,而且運算比浮點快。

一。點,線,面,形基本關系,點積叉積的理解

POJ 2318 TOYS(推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=2318
POJ 2398 Toy Storage(推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=2398
一個矩形,有被若干直線分成N個格子,給出一個點的坐標,問你該點位于哪個點中。
知識點:其實就是點在凸四邊形內的判斷,若利用叉積的性質,可以二分求解。

POJ 3304 Segments
http://acm.pku.edu.cn/JudgeOnline/problem?id=3304
知識點:線段與直線相交,注意枚舉時重合點的處理

POJ 1269 Intersecting Lines
http://acm.pku.edu.cn/JudgeOnline/problem?id=1269
知識點:直線相交判斷,求相交交點

POJ 1556 The Doors (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1556
知識點:簡單圖論+簡單計算幾何,先求線段相交,然后再用Dij求最短路。

POJ 2653 Pick-up sticks
http://acm.pku.edu.cn/JudgeOnline/problem?id=2653
知識點:還是線段相交判斷

POJ 1066 Treasure Hunt
http://acm.pku.edu.cn/JudgeOnline/problem?id=1066
知識點:線段相交判斷,不過必須先理解“走最少的門”是怎么一回事。

POJ 1410 Intersection
http://acm.pku.edu.cn/JudgeOnline/problem?id=1410
知識點:線段與矩形相交。正確理解題意中相交的定義。
詳見:
http://hi.baidu.com/novosbirsk/blog/item/68c682c67e8d1f1d9d163df0.html

POJ 1696 Space Ant (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1696
德黑蘭賽區的好題目。需要理解點積叉積的性質

POJ 3347 Kadj Squares
http://acm.pku.edu.cn/JudgeOnline/problem?id=3347
本人的方法極度猥瑣。復雜的線段相交問題。這個題目是計算幾何的擴大數據運算的典型應用,擴大根號2倍之后就避免了小數。

POJ 2826 An Easy Problem?! (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=2826
問:兩條直線組成一個圖形,能容納多少雨水。很不簡單的Easy Problem,要考慮所有情況。你不看discuss看看能否AC。(本人基本不能)提示一下,水是從天空垂直落下的。

POJ 1039 Pipe
http://acm.pku.edu.cn/JudgeOnline/problem?id=1039
又是線段與直線相交的判斷,再加上枚舉的思想即可。

POJ 3449 Geometric Shapes
http://acm.pku.edu.cn/JudgeOnline/problem?id=3449
判斷幾何體是否相交,不過輸入輸出很惡心。
此外,還有一個知識點,就是給出一個正方形(邊不與軸平行)的兩個對角線上的頂點,需要你求出另外兩個點。必須掌握其方法。

POJ 1584 A Round Peg in a Ground Hole
http://acm.pku.edu.cn/JudgeOnline/problem?id=1584
知識點:點到直線距離,圓與多邊形相交,多邊形是否為凸

POJ 2074 Line of Sight (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=2074
與視線問題的解法,關鍵是求過兩點的直線方程,以及直線與線段的交點。數據有一個trick,要小心。

二。凸包問題

POJ 1113 Wall
http://acm.pku.edu.cn/JudgeOnline/problem?id=1113
知識點:赤裸裸的凸包問題,凸包周長加上圓周。

POJ 2007 Scrambled Polygon
http://acm.pku.edu.cn/JudgeOnline/problem?id=2007
知識點:凸包,按極角序輸出方案

POJ 1873 The Fortified Forest (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1873
World Final的水題,先求凸包,然后再搜索。由于規模不大,可以使用位運算枚舉。
詳見:
http://hi.baidu.com/novosbirsk/blog/item/333abd54c7f22c52574e0067.html

POJ 1228 Grandpa's Estate (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1228
求凸包頂點數目,很多人求凸包的模板是會多出點的,雖然求面積時能得到正確答案,但是在這個題目就會出問題。此外,還要正確理解凸包的性質。

POJ 3348 Cows
http://acm.pku.edu.cn/JudgeOnline/problem?id=3348
凸包面積計算

三。面積問題,公式問題

POJ 1654 Area
http://acm.pku.edu.cn/JudgeOnline/problem?id=1654
知識點:利用有向面積(叉積)計算多邊形面積

POJ 1265 Area
http://acm.pku.edu.cn/JudgeOnline/problem?id=1265
POJ 2954 Triangle
http://acm.pku.edu.cn/JudgeOnline/problem?id=2954
Pick公式的應用,多邊形與整點的關系。(存在一個GCD的關系)

四。半平面交

半平面交的主要應用是判斷多邊形是否存在核,還可以解決一些與線性方程組可行區域相關的問題(就是高中時的那些)。

POJ 3335 Rotating Scoreboard
http://acm.pku.edu.cn/JudgeOnline/problem?id=3335
POJ 3130 How I Mathematician Wonder What You Are!
http://acm.pku.edu.cn/JudgeOnline/problem?id=3130
POJ 1474 Video Surveillance
http://acm.pku.edu.cn/JudgeOnline/problem?id=1474
知識點:半平面交求多邊形的核,存在性判斷

POJ 1279 Art Gallery
http://acm.pku.edu.cn/JudgeOnline/problem?id=1279
半平面交求多邊形的核,求核的面積

POJ 3525 Most Distant Point from the Sea (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=3525
給出一個多邊形,求里面的一個點,其距離離多邊形的邊界最遠,也就是多邊形中最大半徑圓。
可以使用半平面交+二分法解。二分這個距離,邊向內逼近,直到達到精度。

POJ 3384 Feng Shui (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=3384
半平面交實際應用,用兩個圓覆蓋一個多邊形,問最多能覆蓋多邊形的面積。
解法:用半平面交將多邊形的每條邊一起向“內”推進R,得到新的多邊形,然后求多邊形的最遠兩點。

POJ 1755 Triathlon (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1755
半平面交判斷不等式是否有解。注意不等式在轉化時正負號的選擇,這直接影響到半平面交的方向。

POJ 2540 Hotter Colder
http://acm.pku.edu.cn/JudgeOnline/problem?id=2540
半平面交求線性規劃可行區域的面積。

POJ 2451 Uyuw's Concert
http://acm.pku.edu.cn/JudgeOnline/problem?id=2451
Zzy專為他那篇nlogn算法解決半平面交問題的論文而出的題目。

五。計算幾何背景,實際上解題的關鍵是其他問題(數據結構、組合數學,或者是枚舉思想)
若干道經典的離散化+掃描線的題目,ACM選手必做題目

POJ 1151 Atlantis (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1151
POJ 1389 Area of Simple Polygons
http://acm.pku.edu.cn/JudgeOnline/problem?id=1389
矩形離散化,線段樹處理,矩形面積求交

POJ 1177 Picture (推薦

http://acm.pku.edu.cn/JudgeOnline/problem?id=1177
矩形離散化,線段樹處理,矩形交的周長,這個題目的數據比較強。線段樹必須高效。

POJ 3565 Ants (推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=3565
計算幾何中的調整思想,有點像排序。要用到線段相交的判斷。
詳見:
http://hi.baidu.com/novosbirsk/blog/item/fb668cf0f362bec47931aae2.html


POJ 3695 Rectangles   
http://acm.pku.edu.cn/JudgeOnline/problem?id=3695
又是矩形交的面積,但是由于是多次查詢,而且矩形不多,使用組合數學中的容斥原理解決之最適合。線段樹是通法,但是除了線段樹,還有其他可行的方法。

POJ 2002 Squares   
http://acm.pku.edu.cn/JudgeOnline/problem?id=2002
枚舉思想,求平面上若干個點最多能組成多少個正方形,點的Hash

POJ 1434 Fill the Cisterns!(推薦)
http://acm.pku.edu.cn/JudgeOnline/problem?id=1434
一開始發昏了,準備弄個線段樹。其實只是個簡單的二分。

六。隨機算法
POJ 2420 A Star not a Tree?
http://acm.pku.edu.cn/JudgeOnline/problem?id=2420
多邊形的費馬點。所謂費馬點,就是多邊形中一個點P,該點到其他點的距離之和最短。四邊形以上的多邊形沒有公式求費馬點,因此可以使用隨機化變步長貪心法。
詳見:
http://hi.baidu.com/novosbirsk/blog/item/75983f138499f825dd54019b.html

七。解析幾何
這種題目本人不擅長,所以做得不多,模板很重要。當然,熟練運用叉積、點積的性質還是很有用的。
POJ 1375 Intervals
http://acm.pku.edu.cn/JudgeOnline/problem?id=1375
知識點:過圓外一點求與圓的切線

POJ 1329 Circle Through Three Points   
http://acm.pku.edu.cn/JudgeOnline/problem?id=1329
求三角形外接圓

POJ 2354 Titanic
http://acm.pku.edu.cn/JudgeOnline/problem?id=2354
求球面上兩個點的距離,而且給的是地理經緯坐標。

POJ 1106 Transmitters
http://acm.pku.edu.cn/JudgeOnline/problem?id=1106
角度排序,知道斜率求角度,使用atan函數。

POJ 1673 EXOCENTER OF A TRIANGLE
http://acm.pku.edu.cn/JudgeOnline/problem?id=1673
可以轉化為三角形的垂心問題。

八。旋轉卡殼

POJ 2187 Beauty Contest
http://acm.pku.edu.cn/JudgeOnline/problem?id=2187
凸包求最遠點對??梢员┝γ杜e,也可以使用旋轉卡殼。

POJ 3608 Bridge Across Islands(難)
http://acm.pku.edu.cn/JudgeOnline/problem?id=3608
兩個凸包的最近距離。本人的卡殼始終WA。郁悶。

九。其他問題
POJ 1981 Circle and Points
http://acm.pku.edu.cn/JudgeOnline/problem?id=1981
求單位圓最多能覆蓋平面上多少個點

posted @ 2009-08-07 01:15 abilitytao 閱讀(1785) | 評論 (0)編輯 收藏

POJ 1066——Treasure Hunt解題報告

我剛開始被這道題目的名字吸引了,因為它和寶藏有關,呵呵^_^不過把題目讀完以后才發現這道題是個無聊的計算幾何題 ,說實話有點失望。。。

題目的大意是這樣的:尋寶者在一個被分割成很多房間的正方形迷宮里尋寶,這個迷宮是100*100的正方形而且四個頂點坐標一定。尋寶者具有把墻鑿穿通過的能力,若尋寶者可以從正方形的任意一個邊界進入,問到達藏寶地點最少要穿過幾道墻?

這個題的解法是:枚舉每一個入口。然后在所有的情況中取穿墻數最少的輸出即可。
考察每一個入口的時候,枚舉每條邊,如果起點和終點在這條直線的兩側,那么尋寶者一定要穿過一道墻。于是此題轉化成了判斷2點是否在一條直線的異側的問題。模板解決~

由于自己寫的有點冗長,于是參考了下網上的代碼,發現將所有邊界上的點按照角度排序的確是個很巧妙的方法,學習了^_^

//coded by abilitytao
//Time:2009年8月5日17:50:19

#include
<iostream>
#include
<cmath>
#include
<algorithm>
using namespace std;
double const EPS = 1e-8;
const int INF = 0xf777777;
#define zero(x) (((x)>0?(x):-(x))<eps) 



struct Point
{
    
double x,y;
    Point()
{}
       Point(
double a, double b):x(a), y(b){}
       
bool operator<(Point a){return atan2(y - 50, x - 50< atan2(a.y - 50, a.x - 50); }
}


struct Line// 定義一條線段,用起點和終點來表示 
{               
    Point a, b; 
    Line() 
{} 
    Line(Point p10, Point p20): a(p10), b(p20) 
{} //Line a(p1,p2);
}



Point mypoint[
64], s, t;
Line myline[
30];
int n, countnum, minnum;

double xmult(Point p1, Point p2 , Point p0)
{
    
return (p1.x - p0.x)*(p2.y - p0.y)-(p2.x - p0.x)*(p1.y - p0.y);
}

int same_side(Point p1,Point p2,Line l)

    
return xmult(l.a,p1,l.b)*xmult(l.a,p2,l.b)>EPS; 
}


int main()
{
    
int i, j, ans;
    minnum 
= INF; countnum = 0;
    mypoint[countnum
++= Point(00);
    mypoint[countnum
++= Point(1000);
    mypoint[countnum
++= Point(0100);
    mypoint[countnum
++= Point(100100);
    
    cin 
>> n;
    
for(i = 0; i < n; i++)
    
{
        scanf(
"%lf%lf%lf%lf",&myline[i].a.x ,&myline[i].a.y ,&myline[i].b.x ,&myline[i].b.y);
        mypoint[countnum
++= myline[i].a;
        mypoint[countnum
++= myline[i].b;
    }

    scanf(
"%lf%lf",&s.x,&s.y);
    
       sort(mypoint, mypoint
+countnum);
       
       
for(i=0;i<countnum;i++ )
       
{
           ans 
= 0;
           t 
= Point( (mypoint[i].x + mypoint[(i+1)%countnum].x)/2, (mypoint[i].y + mypoint[(i+1)%countnum].y)/2 );
           
for(j = 0; j < n; j++)          
               
if(!same_side(s, t, myline[j]))
                   ans
++;
            
if(ans <minnum) minnum = ans;
       }

       
       printf(
"Number of doors = %d\n", minnum+1);
       
return 0;
}

posted @ 2009-08-05 23:58 abilitytao 閱讀(1578) | 評論 (0)編輯 收藏

2009年ACM-ICPC亞洲區預選賽共設十五個賽區如下(按現場賽日期排序):

1、Harbin 哈爾濱賽區(哈爾濱工業大學)
網絡選拔賽日期:2009年9月13日 14:00-17:00
現場賽日期:2009年9月26日~27日
http://acm.hit.edu.cn/

2、Dhaka 達卡賽區(孟加拉國 Northsouth University)
現場賽日期:2009年10月3日
http://www.northsouth.edu/acm/

3、Gwalior-Kanpur 瓜廖爾-坎普爾賽區(印度 IIITM Gwalior and Indian Institute of Technology, India)
現場賽日期:2009年10月3日~4日
http://www.cse.iitk.ac.in/users/acm/

4、Hefei 合肥賽區(中國科學技術大學)
網絡選拔賽日期:2009年9月6日 14:00-17:00
現場賽日期:2009年10月10日~11日
http://acm.ustc.edu.cn/

5、Ningbo 寧波賽區(浙江大學寧波理工學院)
網絡選拔賽日期:2009年9月12日
現場賽日期:2009年10月17日~18日
http://acmasia09.nit.net.cn/

6、Jakarta 雅加達賽區(印尼 Binus University)
現場賽日期:2009年10月21日
http://icpc.ewubd.edu/

7、Manila 馬尼拉賽區(菲律賓 Ateneo de Manila University)
現場賽日期:2009年10月22日~23日
http://www.math.admu.edu.ph/acm/

8、Shanghai 上海賽區(東華大學)
網絡選拔賽日期:2009年9月20日(12:00-17:00)
現場賽日期:2009年10月24日~25日
http://acm.dhu.edu.cn

9、Hsinchu 新竹賽區(交通大學)
報名截止日期:2009年8月19日
現場賽日期:2009年10月30日~31日
http://icpc2009.nctu.edu.tw/

10、Wuhan 武漢賽區(武漢大學)
網絡選拔賽日期:2009年10月3
現場賽日期:2009年10月31日~11月1日
http://acm.whu.edu.cn/

11、Amritapuri 阿姆里塔普里賽區(印度 Amrita Vishwa Vidyapeetham, Amritapuri Campus)
現場賽日期:2009年11月1日
http://icpc.amrita.ac.in/

12、Phuket 普吉島賽區(泰國 Prince of Songkla University, Phuket Campus)
報名截止日期:2009年9月30日
現場賽日期:2009年11月3日~4日
http://www.acmicpc-thailand.org/

13、Seoul 首爾賽區(韓國 Korea Advanced Institute of Science and Technology)
報名截止日期:2009年9月11日
現場賽日期:2009年11月5日~6日
http://acm.kaist.ac.kr/

14、Tehran 德黑蘭賽區(伊朗 Sharif University of Technology)
現場賽日期:2009年11月6日
http://sharif.edu/~acmicpc

15、Tokyo 東京賽區(日本早稻田大學)
現場賽日期:2009年11月7日~8日
http://www.waseda.jp/assoc-icpc2009/en/index.html

參賽報名網址

http://cm.baylor.edu/welcome.icpc

亞洲高??山M隊參加全部十五個賽區的預選賽,但每位參賽選手最多只能注冊參加兩個賽區的預選賽。

posted @ 2009-08-04 17:08 abilitytao 閱讀(1739) | 評論 (0)編輯 收藏

計算幾何模板整理

#include<iostream>
#include
<cmath>
#include
<algorithm>
using namespace std;

struct Point // 定義一個點
{              
    
double x, y;  
    Point() 
{} 
    Point(
double x0, double y0): x(x0), y(y0) {} //Point a(1.0,2.0);
}


struct Line// 定義一條線段,用起點和終點來表示 
{               
    Point p1, p2; 
    Line() 
{} 
    Line(Point p10, Point p20): p1(p10), p2(p20) 
{} //Line a(p1,p2);
}


//定義叉積:
//1.如果返回值為正數,表明sp在op->ep(op指向ep)這條射線的順時針方向;
//2.如果返回值為負數,表明sp在op->ep(op指向ep)這條射線的逆時針方向;
//3.如果返回值為0,表明三點共線
double multiply(Point sp,Point ep,Point op)
{
    
return((sp.x-op.x)*(ep.y-op.y) - (ep.x-op.x)*(sp.y-op.y));
}

posted @ 2009-08-04 16:35 abilitytao 閱讀(234) | 評論 (0)編輯 收藏

POJ 2318 toys 叉積的簡單運用

題目意思很簡單,把一個盒子分成很多部分,往盒子里扔小球,小球的落點當然要告訴你拉,讓你統計每個盒子里最后擁有的小球數。
我的做法是叉積+二分。
#include<iostream>
#include
<cmath>
#include
<algorithm>
using namespace std;

struct Point 
{              // 二維點或矢量 
    double x, y;  
    Point() 
{} 
    Point(
double x0, double y0): x(x0), y(y0) {} 
}


struct Line
 
{               // 二維的直線或線段 
    Point p1, p2; 
    Line() 
{} 
    Line(Point p10, Point p20): p1(p10), p2(p20) 
{} 
}



double multiply(Point sp,Point ep,Point op)
{
    
return((sp.x-op.x)*(ep.y-op.y) - (ep.x-op.x)*(sp.y-op.y));
}




Line myline[
5100];
int record[5100];


int main()
{

    
int n,m,i;
    Point left;
    Point right;
    
while(scanf("%d",&n))
    
{

        
if(n==0)
            
break;
        memset(record,
0,sizeof(record));
        scanf(
"%d%lf%lf%lf%lf",&m,&left.x,&left.y,&right.x,&right.y);
        myline[
0].p1.x=left.x;
        myline[
0].p1.y=right.y;
        myline[
0].p2.x=left.x;
        myline[
0].p2.y=left.y;
        
for(i=1;i<=n;i++)
        
{
            scanf(
"%lf%lf",&myline[i].p2.x,&myline[i].p1.x);
            myline[i].p1.y
=right.y;
            myline[i].p2.y
=left.y;

        }

        myline[n
+1].p1.x=right.x;
        myline[n
+1].p1.y=right.y;
        myline[n
+1].p2.x=right.x;
        myline[n
+1].p2.y=left.y;
        
        Point toy;
        
for(i=1;i<=m;i++)
        
{
            scanf(
"%lf%lf",&toy.x,&toy.y);
            
int front=0;
            
int rear=n+1;
            
while(front<=rear)
            
{

                
int mid=(front+rear)>>1;
                
if(multiply(toy,myline[front].p2,myline[front].p1)>0&&multiply(toy,myline[mid].p2,myline[mid].p1)<0)
                
{

                    
if(mid==front+1)
                    
{
                        record[front]
++;
                        
break;
                    }

                    rear
=mid;
                    
continue;
                }

                
else
                
{

                    
if(mid+1==rear)
                    
{
                        record[mid]
++;
                        
break;
                    }

                    front
=mid;
                    
continue;
                }


            }

        }


        
for(i=0;i<=n;i++)
        
{
            printf(
"%d: %d\n",i,record[i]);
        }

        printf(
"\n");
    }

    
return 0;
}


恭喜此題成為我收集計算幾何模板的第一題 呵呵~

posted @ 2009-08-04 16:26 abilitytao 閱讀(956) | 評論 (1)編輯 收藏

僅列出標題
共42頁: First 26 27 28 29 30 31 32 33 34 Last 
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            aa日韩免费精品视频一| 亚洲图片欧洲图片av| 亚洲成人在线视频播放| 国产亚洲精品aa午夜观看| 国产精品一二三四区| 国产精品美女久久久久久免费| 欧美日韩一区在线播放| 欧美性猛交xxxx乱大交蜜桃 | 亚洲精品之草原avav久久| 欧美激情一区二区三区成人| 欧美激情一区二区三级高清视频| 亚洲国产成人tv| 久久精品91| 亚洲大片一区二区三区| 亚洲日本国产| 午夜精品视频在线观看一区二区| 欧美中文字幕| 欧美精品一区二区视频 | 亚洲高清视频在线| 亚洲另类自拍| 欧美在线一区二区三区| 欧美aa在线视频| 9色精品在线| 久久久久一区二区三区| 欧美日韩国产综合网| 国产日韩欧美麻豆| 日韩一级不卡| 久久亚洲一区二区| 亚洲色图综合久久| 欧美成人a∨高清免费观看| 国产精品视频yy9299一区| 亚洲国产经典视频| 欧美在线免费视频| 日韩图片一区| 免费成人黄色av| 国产一区二区三区黄| 在线亚洲精品| 亚洲成色www8888| 欧美一区日本一区韩国一区| 欧美日本国产视频| 亚洲激情黄色| 久久久www成人免费毛片麻豆| 亚洲国产精品久久久久秋霞不卡 | 99亚洲视频| 欧美丰满高潮xxxx喷水动漫| 国产麻豆视频精品| 在线一区二区日韩| 亚洲激情网站免费观看| 久久久999精品免费| 国产精品电影在线观看| 亚洲欧洲一区| 欧美成人午夜| 久久婷婷丁香| 激情一区二区| 美女黄网久久| 久久艳片www.17c.com| 国产一区在线免费观看| 久久精品一区二区国产| 亚洲欧美视频在线| 国产精品一区在线观看| 亚洲欧美日韩精品在线| 在线亚洲观看| 国产精品久久久91| 亚洲一区欧美激情| 亚洲一区二区三区777| 国产精品国产三级国产aⅴ9色 | 欧美激情视频一区二区三区免费| 激情综合色综合久久| 久久看片网站| 久久国产手机看片| 国产亚洲va综合人人澡精品| 久久av二区| 久久人人爽人人爽| 亚洲激情自拍| 日韩视频不卡| 欧美性猛交一区二区三区精品| 一区二区欧美视频| 亚洲男女自偷自拍| 黄色在线一区| 亚洲人成网在线播放| 欧美日本在线| 欧美在线视屏| 久久综合久久综合这里只有精品 | 久久久久久九九九九| 亚洲精品免费网站| 在线亚洲伦理| 国产综合18久久久久久| 欧美ab在线视频| 国产精品久久久久7777婷婷| 久久―日本道色综合久久| 免费欧美日韩| 欧美一级电影久久| 欧美成人亚洲成人日韩成人| 亚洲一二三区在线观看| 久久精品视频在线播放| 夜夜嗨一区二区| 久久精品国产久精国产思思| 一本色道**综合亚洲精品蜜桃冫| 午夜精品久久久久久久 | 久久一日本道色综合久久| 蜜臀av性久久久久蜜臀aⅴ四虎 | 激情婷婷亚洲| 日韩视频在线免费观看| 国产香蕉97碰碰久久人人| 亚洲国产精品一区制服丝袜| 国产日韩欧美精品| 一本综合精品| 亚洲毛片在线看| 久久精品国内一区二区三区| 亚洲视频日本| 欧美jizz19性欧美| 久久国产精品久久精品国产| 欧美激情综合| 欧美成人按摩| 国产一区二区三区高清播放| 夜色激情一区二区| 亚洲七七久久综合桃花剧情介绍| 午夜精品视频一区| 午夜国产精品视频| 欧美三级电影大全| 91久久久在线| 亚洲欧洲精品天堂一级| 久久成人人人人精品欧| 亚洲欧美日韩电影| 欧美日韩综合精品| 亚洲欧洲一区二区在线播放| 伊人色综合久久天天五月婷| 亚洲——在线| 亚洲在线观看免费| 欧美三级视频| 999亚洲国产精| 99精品欧美一区二区蜜桃免费| 美女视频黄免费的久久| 欧美1区3d| 亚洲日韩欧美一区二区在线| 老司机精品视频网站| 欧美激情a∨在线视频播放| 在线观看一区欧美| 久久高清国产| 嫩草成人www欧美| 亚洲国产精品久久91精品| 久久只精品国产| 老司机午夜免费精品视频| 国产在线视频欧美一区二区三区| 亚洲欧美在线免费| 久久精品人人做人人综合| 国产午夜精品一区理论片飘花 | 欧美激情精品久久久久久| 欧美激情一区二区三区高清视频| 亚洲人成久久| 牛牛国产精品| 亚洲精品乱码久久久久久按摩观| 一区二区激情小说| 国产精品久久久久久av福利软件| 一区二区三区 在线观看视| 午夜精彩视频在线观看不卡| 国产精品专区第二| 久久成人免费| 欧美成人黄色小视频| 亚洲精品字幕| 国产精品久久9| 久久成人免费电影| 亚洲精品在线一区二区| 亚洲午夜激情网站| 国产日韩一区| 欧美高清在线一区二区| 日韩亚洲欧美一区二区三区| 一区二区冒白浆视频| 欧美日韩亚洲天堂| 久久精品日韩| 中文亚洲字幕| 蜜月aⅴ免费一区二区三区 | 欧美在线观看一区二区| 国产日韩av高清| 免费亚洲婷婷| 亚洲欧美成人综合| 亚洲激情在线激情| 久久久99免费视频| 亚洲欧美变态国产另类| 亚洲二区三区四区| 国产精品一区久久| 欧美精品一区在线播放| 久久狠狠婷婷| 亚洲图色在线| 亚洲国产天堂网精品网站| 久久www免费人成看片高清 | 亚洲自拍偷拍麻豆| 亚洲夫妻自拍| 国产区在线观看成人精品| 欧美精品在线视频| 久久久久久久尹人综合网亚洲 | 亚洲一本视频| 最近中文字幕日韩精品| 久色婷婷小香蕉久久| 亚洲欧美另类中文字幕| 亚洲美女色禁图| 亚洲激情成人在线| 黑人巨大精品欧美黑白配亚洲| 国产精品国产三级国产普通话蜜臀 | 久久精品72免费观看|