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

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>
            国产精品青草综合久久久久99| 亚洲日本电影| 亚洲综合色噜噜狠狠| aa国产精品| 一区二区三区四区五区视频| 亚洲电影免费| 欧美国产极速在线| 欧美激情a∨在线视频播放| 另类尿喷潮videofree | 亚洲大片在线| 亚洲日韩欧美一区二区在线| 亚洲一区二区三区四区视频| 午夜国产精品影院在线观看| 午夜久久一区| 亚洲国产小视频在线观看| 亚洲免费精品| 久久爱另类一区二区小说| 欧美成人免费va影院高清| 欧美劲爆第一页| 国产精品嫩草影院一区二区| 欧美日韩中文字幕在线视频| 欧美日韩中文字幕| 久久久综合免费视频| 麻豆精品一区二区av白丝在线| 欧美1区2区| 国产九区一区在线| 亚洲精品综合精品自拍| 亚洲午夜在线| 99这里有精品| 韩国一区二区三区在线观看 | 亚洲一区成人| 在线视频国产日韩| 亚洲国产国产亚洲一二三| 国产字幕视频一区二区| 狠狠狠色丁香婷婷综合激情| 国语自产在线不卡| 亚洲一区自拍| 亚洲一卡二卡三卡四卡五卡| 一区二区欧美国产| 美乳少妇欧美精品| 久久一区二区三区av| 国产精品每日更新| 这里只有精品视频| 欧美大色视频| 久久成人国产| 影音先锋成人资源站| 亚洲欧美清纯在线制服| 亚洲精品久久7777| 久久这里有精品15一区二区三区 | 国内精品久久久| 亚洲欧美日韩国产另类专区| 午夜精品一区二区三区在线播放| 国语精品中文字幕| 欧美大片免费观看| 精品福利av| 欧美岛国在线观看| 亚洲精品少妇网址| 国产性色一区二区| 日韩视频在线一区| 亚洲国产精品精华液2区45| 玖玖精品视频| 亚洲欧美综合精品久久成人| 久久精品免费| 激情综合在线| 99热免费精品| 亚洲视频在线观看网站| 国产精品人人爽人人做我的可爱| 一区二区精品| 午夜精品理论片| 一区二区亚洲精品| 国产精品久久一卡二卡| 先锋影音久久| 玖玖玖国产精品| 亚洲欧美一区二区视频| 欧美一区二区三区免费观看视频 | 欧美成人资源| 亚洲理论在线观看| 一本一本久久a久久精品综合麻豆| 欧美午夜视频在线观看| 久久av资源网站| 欧美色道久久88综合亚洲精品| 久久精品国亚洲| 国产精品成人免费精品自在线观看| 久久九九99| 国产一区二区精品| 亚洲一区欧美一区| 一区二区三区产品免费精品久久75| 久久精品人人爽| 美女性感视频久久久| 国产综合欧美| 久久久久久尹人网香蕉| 久久婷婷国产综合国色天香| 欧美激情精品久久久久久免费印度 | 亚洲美女免费视频| 欲色影视综合吧| 欧美与欧洲交xxxx免费观看| 性欧美1819sex性高清| 久久久久综合网| 亚洲欧美成人综合| 欧美午夜激情视频| 亚洲自拍啪啪| 国产精品成人免费精品自在线观看| 欧美高清你懂得| 一本一本久久a久久精品综合麻豆 一本一本久久a久久精品牛牛影视 | 亚洲综合国产激情另类一区| 欧美精品一区在线| 日韩西西人体444www| 亚洲免费在线播放| 国产亚洲毛片在线| 免费高清在线一区| 亚洲永久在线观看| 久久一区免费| 亚洲一区精品在线| 日韩一二在线观看| 午夜精品久久久久久久久| 韩国女主播一区| 美女网站久久| 一区二区不卡在线视频 午夜欧美不卡在 | 亚洲激情在线视频| 亚洲精品久久久久久久久久久久久| 欧美国产大片| 久久精品国产综合| 欧美视频在线观看一区| 亚洲免费在线看| 欧美日韩国产大片| 亚洲第一精品夜夜躁人人爽| 国产视频亚洲精品| 亚洲欧美激情在线视频| 亚洲欧美国产高清| 国产精品成人免费| 亚洲精品一区二区三区婷婷月 | 欧美日韩在线播放三区| 亚洲国产aⅴ天堂久久| 在线日韩欧美| 鲁大师成人一区二区三区| 欧美jizz19性欧美| 亚洲电影在线免费观看| 久久久久久久久岛国免费| 乱中年女人伦av一区二区| 精品白丝av| 蜜桃视频一区| 亚洲精品乱码久久久久| 中文一区二区在线观看| 欧美视频在线不卡| 亚洲欧美经典视频| 久久视频国产精品免费视频在线| 国产主播一区二区| 美女视频一区免费观看| 亚洲高清自拍| 亚洲欧美日韩区| 国内揄拍国内精品久久| 久热精品视频在线观看| 亚洲精品美女91| 欧美一区=区| 尤物99国产成人精品视频| 欧美成人免费一级人片100| 亚洲免费观看高清在线观看 | 国产日韩欧美在线看| 久久精品免费播放| 亚洲高清成人| 亚洲欧美日产图| 狠狠干狠狠久久| 欧美日韩大片一区二区三区| 亚洲一二三四区| 免费中文字幕日韩欧美| 一区二区三区视频在线| 国产网站欧美日韩免费精品在线观看 | 久久精品男女| 亚洲免费播放| 久久夜色精品亚洲噜噜国产mv| 亚洲欧洲午夜| 国产欧美日韩综合| 欧美激情按摩在线| 欧美亚洲一区| 亚洲理论电影网| 久久综合亚洲社区| 久久精品视频一| 亚洲一区二区综合| 国产精品美女久久久久久2018| 欧美一区久久| 99香蕉国产精品偷在线观看| 欧美 日韩 国产在线| 欧美亚洲网站| 亚洲一区二区三区777| 91久久黄色| 韩国在线视频一区| 国产精品看片资源| 欧美日韩不卡| 欧美大成色www永久网站婷| 欧美一区二视频在线免费观看| 亚洲毛片播放| 亚洲激情中文1区| 欧美韩日一区| 暖暖成人免费视频| 久久噜噜亚洲综合| 性欧美超级视频| 亚洲自拍另类| 亚洲综合色在线| 亚洲一区二区在线免费观看| 一本高清dvd不卡在线观看|