2009年5月18日
今天組隊比賽時遇到了漢諾塔的題，開始時以為很難，但是仔細看題時，想到了曾今在杭電上做個一題，卻發現是水題一道，直接一個打表！
TJU
1731. Strange Towers of Hanoi
#include<stdio.h>
int main()
{
int a[13]={0,1,3,5,9,13,17,25,33,41,49,65,81};
for(int i=1;i<=12;i++)
printf("%d\",a[i]);
}真正的題目還是杭電上的
http://acm.hdu.edu.cn/showproblem.php?pid=1207
#include<stdio.h>
#include<math.h>
int main()
{
long a[65];
a[1]=1;
long i,j,n,m,r;
for(i=2,j=1,r=1,m=2;i<66;i++,j++)
if(j<=m)
a[i]=a[i-1]+pow(2,r);
else {
m++;
j=1;
r++;
a[i]=a[i-1]+pow(2,r);
}
while(scanf("%d",&n)==1){
printf("%d\n",a[n]);
}
}
規律：
a[1]=1;
a[2]=a[1]+2;a[3]=a[2]+2;(2個加2^1)
a[4]=a[3]+4;a[5]=a[4]+4;a[6]=a[5]+4;(3個加2^2);
…………………………………………(4個加2^3);
O(∩_∩)O~
posted @
2009-05-18 18:32 zhoubaozhong 閱讀(944) |
評論 (0) |
編輯 收藏【圖論】
1、Dijkstra算法
2、Floyd算法
3、Kruskal算法
4、Prim算法
5、歐拉回路
Dijkstra算法
Dijkstra算法是典型最短路算法，用于計算一個節點到其他所有節點的最短路徑。主要特點是以起始點為中心向外層層擴展，直到擴展到終點為止。Dijkstra算法能得出最短路徑的最優解，但由于它遍歷計算的節點很多，所以效率低。
Dijkstra算法是很有代表性的最短路算法，在很多專業課程中都作為基本內容有詳細的介紹，如數據結構，圖論，運籌學等等。
Dijkstra一般的表述通常有兩種方式，一種用永久和臨時標號方式，一種是用OPEN, CLOSE表方式，Drew為了和下面要介紹的 A* 算法和 D* 算法表述一致，這里均采用OPEN,CLOSE表的方式。
大概過程：
創建兩個表，OPEN, CLOSE。
OPEN表保存所有已生成而未考察的節點，CLOSED表中記錄已訪問過的節點。
1． 訪問路網中里起始點最近且沒有被檢查過的點，把這個點放入OPEN組中等待檢查。
2． 從OPEN表中找出距起始點最近的點，找出這個點的所有子節點，把這個點放到CLOSE表中。
3． 遍歷考察這個點的子節點。求出這些子節點距起始點的距離值，放子節點到OPEN表中。
4． 重復2，3，步。直到OPEN表為空，或找到目標點。
Floyd算法
Floyd-Warshall 算法用來找出每對點之間的最短距離。它需要用鄰接矩陣來儲存邊，這個算法通過考慮最佳子路徑來得到最佳路徑。
注意單獨一條邊的路徑也不一定是最佳路徑。
從任意一條單邊路徑開始。所有兩點之間的距離是邊的權，或者無窮大，如果兩點之間沒有邊相連。
對于每一對頂點 u 和 v，看看是否存在一個頂點 w 使得從 u 到 w 再到 v 比己知的路徑更短。如果是更新它。
不可思議的是，只要按排適當，就能得到結果。
// dist(i,j) 為從節點i到節點j的最短距離
For i←1 to n do
For j←1 to n do
dist(i,j) = weight(i,j)
For k←1 to n do // k為“媒介節點”
For i←1 to n do
For j←1 to n do
if (dist(i,k) + dist(k,j) < dist(i,j)) then // 是否是更短的路徑？
dist(i,j) = dist(i,k) + dist(k,j)這個算法的效率是O(V3)。它需要鄰接矩陣來儲存圖。
這個算法很容易實現，只要幾行。
即使問題是求單源最短路徑，還是推薦使用這個算法，如果時間和空間允許(只要有放的下鄰接矩陣的空間，時間上就沒問題)。
[編輯] 時間復雜度 O(N^3)
Kruskal算法
基本思想
假設WN=(V,{E})是一個含有n個頂點的連通網，則按照克魯斯卡爾算法構造最小生成樹的過程為：先構造一個只含n個頂點，而邊集為空的子圖，若將該子圖中各個頂點看成是各棵樹上的根結點，則它是一個含有n棵樹的一個森林。之后，從網的邊集E中選取一條權值最小的邊，若該條邊的兩個頂點分屬不同的樹，則將其加入子圖，也就是說，將這兩個頂點分別所在的兩棵樹合成一棵樹；反之，若該條邊的兩個頂點已落在同一棵樹上，則不可取，而應該取下一條權值最小的邊再試之。依次類推，直至森林中只有一棵樹，也即子圖中含有n-1條邊為止。</p>
Procedure kruskal(V,E);
begin
sort(E,1,m);//將邊按照權值排序
for t:=1 to n do begin
if getfather(edge[t].u)<>getfather(edge[t].v) then begin //利用并查集判斷兩個頂點是否在同一集合內
tot:=tot+edge[t].data;//計算權值和
union(edge[t].u,edge[t].v);//合并頂點
inc(k);//合并次數
end;
end;
if k=n-1 then 形成了一棵最小生成樹
else 不存在這樣的最小生成樹；
end;
優化：在判斷兩個頂點是否在同一集合內時可用并查集
Prim算法
基本思想
1. 在圖G=(V, E) （V表示頂點 ，E表示邊）中，從集合V中任取一個頂點（例如取頂點v0）放入集合 U中，這時 U={v0}，集合T(E)為空。
2. 從v0出發尋找與U中頂點相鄰（另一頂點在V中）權值最小的邊的另一頂點v1，并使v1加入U。即U={v0,v1 }，同時將該邊加入集合T(E)中。
3. 重復2，直到U=V為止。
這時T(E)中有n-1條邊，T = (U, T(E))就是一棵最小生成樹。
PASCAL代碼
procedure prim(v0:integer);
var
lowcost,closest:array[1..maxn] of integer;
i,j,k,min:integer;
begin
for i:=1 to n do begin
lowcost[i]:=cost[v0,i];
closest[i]:=v0;
end;
for i:=1 to n-1 do begin
{尋找離生成樹最近的未加入頂點k}
min:=maxlongint;
for j:=1 to n do
if (lowcost[j]<min) and (lowcost[j]<>0) then begin
min:=lowcost[j];
k:=j;
end;
lowcost[k]:=0; {將頂點k加入生成樹}
{生成樹中增加一條新的邊k到closest[k]}
{修正各點的lowcost和closest值}
for j:=1 to n do
if cost[k,j]<lowcost[j] then begin
lowcost[j]:=cost[k,j];
closest[j]:=k;
end;
end;
end;{prim}
〖返回頂部〗
歐拉回路
定義：
在一個圖中，尋找一條只通過每條邊一次的路徑，這叫做歐拉路徑，如果起點和終點是同一點，那么這條回路叫做歐拉回路．
判定一個圖中是否存在歐拉回路：并不是每個圖都存在歐拉回路．以下分三種情況：
無向圖：每個點的度數是偶數，則存在歐拉回路．
有向圖：每個結點的入度等于出度，則這個有向圖中存在歐拉回路．
總結：以上兩種情況很簡單，其原理歸根結底是每個結點進入和出去的路徑條數相等，就存在歐拉回路．還有一種更加復雜的情況．那就是混合圖．
混合圖：（有時邊是單向的，有時邊是無向的，就像城市交通網絡，有的街道是單向的，有的街道是雙向的）找一個給每個無向邊定向的策略，這樣就可以把圖轉化成為有向圖，使每個頂點的入度與出度是相等的，這樣就會存在歐拉回路．
具體過程如下：新建一個圖，對于原圖中的無向邊，在新圖中新加一個頂點e(i);對于原圖中的每一個頂點j，在新圖中建一個頂點v(i)，對于原圖中每一個頂點j和k之間有一條無向邊i，那么在新圖中從e(i)出發，添加兩條邊，分別連向v(j)和v(k)，容量都是1．
在新圖中，從源點向所有的e(i)都連一條容量為1的邊．. 對于原圖中每一個頂點j，它原本都有一個入度in、出度out和無向度un。顯然我們的目的是要把所有無向度都變成入度或出度，從而使它的入度等于總度數的一半，也就是(in + out + un) / 2（顯然與此同時出度也是總度數的一半，如果總度數是偶數的話）。當然，如果in已經大于總度數的一半，或者總度數是奇數，那么歐拉回路肯定不存大。如果in小于總度數的一半，并且總度數是偶數，那么我們在新圖中從v（j）到匯點連一條邊，容量就是(in + out + un) / 2 – in，也就是原圖中頂點j還需要多少入度。
按照這個網絡模型算出一個最大流，如果每條從v(j)到匯點的邊都達到滿流量的話，那么歐拉回路成立。
〖返回頂部〗
posted @
2009-05-18 11:06 zhoubaozhong 閱讀(263) |
評論 (0) |
編輯 收藏
2009年5月16日
http://acm.tju.edu.cn/toj/showp2794.html
#include<stdio.h>
long max(long a,long b)
{
if(a>b)return a;
return b;
}
int main()
{
long i,n,m,ma;
static long a[1000001],b[1000001];
scanf("%d",&m);
while(m--){
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%d",&a[i]);
for(i=0;i<n;i++)
scanf("%d",&b[i]);
a[1]+=a[0];
b[1]+=b[0];
for(i=2;i<n;i++){
a[i]+=max(a[i-1],b[i-2]);
b[i]+=max(b[i-1],a[i-2]);
}
ma=max(b[n-1],a[n-1]);
printf("%d\n",ma);
}
}
posted @
2009-05-16 10:25 zhoubaozhong 閱讀(142) |
評論 (0) |
編輯 收藏
2009年5月4日
#include<stdio.h>
long f(__int64 n)
{
while(n!=0){
if(n%10==7)return 1;
n/=10;
}
return 0;
}
int main()
{
static __int64 a[1000000],n,m,i,j=0,f1;
for(i=7;i<1000000000;i++){
if(f(i)||i%7==0)
a[j]++;
else {
f1=1;
for(int k=0;k<j;k++)
if(a[j]<=a[k])f1=0;
if(f1){printf("p=%I64d,x=%I64d\n",a[j],i-a[j]);
j++;}
else a[j]=0;
}
}
}
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3318比賽時沒有想到，就是靠人腦舉例，導致考慮不完整，出現很多錯誤！后來經高人指點，編學了個輔助程序，終于ac了！
O(∩_∩)O~，看來還是要好好利用電腦！！
主程序：
#include<stdio.h>
int main()
{
long n,t;
scanf("%d",&t);
while(t--){
scanf("%d",&n);
if(n==1)printf("7\n");
else if(n==2)printf("27\n");
else if(n<=10)printf("70\n");
else if(n<=11)printf("270\n");
else if(n<=100)printf("700\n");
else if(n<=101)printf("2700\n");
else if(n<=1000)printf("7000\n");
else if(n<=1002)printf("26999\n");
else if(n<=10000)printf("70000\n");
else if(n<=10001)printf("270000\n");
else if(n<=100000)printf("700000\n");
else if(n<=100001)printf("1699999\n");
else if(n<=1000000)printf("7000000\n");
else if(n<=1000001)printf("27000000\n");
else if(n<=10000000)printf("70000000\n");
else if(n<=10000001)printf("270000000\n");
else printf("700000000\n");
}
}
posted @
2009-05-04 22:14 zhoubaozhong 閱讀(456) |
評論 (0) |
編輯 收藏
2009年4月29日
#include<stdio.h>
int main()
{
int max,min,n,i,x,t,a[10001],sum1,sum2,s;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
max=0;
min=10001;
s=0;
for(i=0;i<n;i++)
{
scanf("%d",&a[i]);
if(a[i]>max)
max=a[i];
if(a[i]<min)
min=a[i];
if(a[i]>=n&&a[i]<0)s=1;
}
if(s==1||n==1&&max==0){ printf("Impossible!\n"); continue;}
sum1=sum2=0;
if(max==min+1){
for(i=0;i<n;i++){
if(max==a[i])sum1++;
else sum2++;
}
if(sum2==max)printf("%d\n",max);
else printf("Impossible!\n");
}
else if(max==min)
{ if(max==0)printf("0\n");
else {
if(max==n-1)printf("%d\n",n);
else printf("Impossible!\n"); }
}
else
printf("Impossible!\n");
}
}
這題是在組隊時候做的！需要考慮很多情況！特別 0 和 1 ！！
http://acm.hdu.edu.cn/showproblem.php?pid=2690
posted @
2009-04-29 21:01 zhoubaozhong 閱讀(250) |
評論 (0) |
編輯 收藏
2009年4月27日
http://acm.hdu.edu.cn/showproblem.php?pid=2687
#include<stdio.h>
int main()
{
int n,i,k,j,a[11][11],b[11][11];
while(scanf("%d",&n)!=EOF)
{
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
scanf("%d",&a[i][j]);
scanf("%d",&k);
for(i=1;i<=n;i++)
for(j=1;j<=n;j++) //首先求解每循環四次的總和*(k/4)
b[i][j]=(a[i][j]+a[j][n+1-i]+a[n+1-i][n+1-j]+a[n+1-j][i])*(k/4);
if(k%4==0)//再求解k%4的結果
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
b[i][j]+=a[i][j];
else if(k%4==1)
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
b[i][j]+=a[i][j]+a[n+1-j][i];
else if(k%4==2)
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
b[i][j]+=a[i][j]+a[n+1-j][i]+a[n+1-i][n+1-j];
else
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
b[i][j]+=a[i][j]+a[n+1-j][i]+a[n+1-i][n+1-j]+a[j][n+1-i];
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
printf("%d",b[i][j]);
if(j<n)
printf(" ");}
printf("\n");
}
}
}
這題只要分清對應的轉動位置就能求解！還要注意轉動順序！！！
posted @
2009-04-27 20:41 zhoubaozhong 閱讀(314) |
評論 (0) |
編輯 收藏
2009年4月21日
http://acm.hdu.edu.cn/showproblem.php?pid=1210開始尋找規律，找了好久想不出，最終想到了追蹤第一張牌來進行求解！O(∩_∩)O~
#include<stdio.h>
int main()
{
long i,n,sum;
while(scanf("%d",&n)==1){
sum=1;
i=2;//追蹤 1 的位置;
while(i!=1){
if(i>n)
i=(i-n)*2-1;
else i*=2;
sum++;//表示排序的次數
}
printf("%d\n",sum);
}
}
雖然做時不知道為什么！但是事后我在網上搜到了證明方法！
洗牌問題：
很多人都是跟蹤第一張牌來完成的，大家會不會證明呢？
定理1：當第一張牌(牌1)回到初始位置時，所有的牌都回到了初始位置。
證明：設有一次操作，某牌在操作前處于位置r(1<=r<=2*N)，那么，操作后，如果r原來是前一半的，顯然r'=r*2; 否則，r'=(r-n)*2-1，即r'='r*2-(N*2+1);
將兩個式子綜合，可以得到r'= (r*2)%(N*2+1);
根據同余定理之一 ((i%m)*j)%m=(i*j)%M，可以整理出公式：i次操作后，該牌的位置為：
(r*(2^i))%(N*2+1);其中^表示乘方。
現在，我們假設經過M次操作后，牌1回到了初始位置，即(1*(2^M))%(N*2+1)=1時，再次應用同余定理，可以證明對于任意的k(1<=k<=2*N)，有(k*(2^M))%(N*2+1)=k，翻譯成自然語言就是：當第一張牌回到初始位置時，所有的牌都回到了初始位置。命題得證。
定理2：一定存在某次操作M，這次操作后，所有的牌都回到了初始位置。
證明：我們已經證明了牌1回到初始位置時，所有的牌都回到初始位置，所以我們只需要證明牌1可以回到初始位置，即(2^M)%(N*2+1)=1一定有正整數解。證明這個定理前我們首先證明這樣一個定理：
定理2.1：(2*r)%(2*n+1)==t
當t確定為1到2*n中的某值時(1<t<=2*n)，r在1到2*n間有唯一解
因為2*n+1是奇數，我們只要看一下就會發現r到t是一個一一映射，當t為偶數時，顯然r=t/2，當t為奇數時，r=(t+1)/2+n，
現在我們來證明定理2。運用反正法，若牌1永遠不能回到初始位置，根據鴿籠定理，牌1必然陷入了某個不包含位置1的循環，因為下一狀態僅和當前狀態有關，當前狀態最多只有2*N種，所以顯然一定在不超過2*N+1次操作內出現重復狀態。而重復狀態則意味著循環。
因為我們假設這一循環不包括位置1，我們設f(i)為循環中某狀態，f(0)=1,f(n+1)=(f(n)*2)%(2*N+1)，f(j)為若干次循環后出現的重復狀態，即f(i)=f(j)。因為循環一直持續，不妨假設j>i+i。又因為循環不包括狀態1，即對于任意的k>=i，f(k)!=1
根據定理2.1，我們可以根據當前狀態推出上一狀態，換句話說，若f(i)=f(j)，則f(i-1)=f(j-1)。繼續套用定理2.1，可以得到f(i-i)=f(j-i)，即：f(j-i)=f(0)=1。又有假設j>i+i，即j-i>i，與另一假設對于任意的k>=i，f(k)!=1矛盾。
因此，可以證明，牌1所陷入的循環，必然是包含位置1的，即一定有某次操作，操作后牌1回到了初始狀態。而且顯然的，初始狀態就是這循環中的一個狀態。
因此命題得證。
posted @
2009-04-21 16:34 zhoubaozhong 閱讀(664) |
評論 (0) |
編輯 收藏
2009年4月18日
因為我是菜鳥！所以我一無是處!
因為我是菜鳥！所以我只會暢想！！！
因為我是菜鳥！所以我只是大家的累贅！！！！！
因為我是菜鳥！所以我不會在尋求改變！！！！！！！
因為我是菜鳥！所以我改變不了什么只有自己去承受！！！！！！！
不就是個人嗎！！不就是缺少討論嗎！！大家不需要,干嘛還要自作多情呢！！！！！！！！
我還年輕！我玩的起！！！！！！！！！！！！！！！！！！！！！
posted @
2009-04-18 21:26 zhoubaozhong 閱讀(224) |
評論 (10) |
編輯 收藏