USACO 2.1 分析

第一題 The Castle

算法：求連通分量的個數

題目給出的形式很好，沒有以強的形式給出，而是以1，2，4，8來對應二進制給出，還要打牆~~~

 

難點1：怎樣判斷有沒有牆？

位運算可以發揮作用了！

 

難點2：怎樣敲牆？

由於題目中有西是最優的，其次是南，所以枚舉的順序有關係；其次，要先敲N，再去敲E，可以寫兩個FOR，也可以把它過程化；

小細節：敲的時候兩邊要屬於不同的連通分量啊！

 

這道題目比以前的代碼大大縮減，只有60++，很高興！
#include <stdio.h> #include <string.h> const int max = 50 + 3, dx[4] = {0, -1, 0, 1}, dy[4] = {-1, 0, 1, 0}; int num[max][max], map[max][max], s[max*max], qx[max*max], qy[max*max]; int n, m; int newMaxArea = 0, newx, newy, newk; void calc(int j, int k) { for (int i = n; i > 0; i--) // wall second farthest to the south if ( (map[i][j] & (1 << k)) ) { int nowx = i+dx[k], nowy = j+dy[k]; if (num[nowx][nowy] != num[i][j] && s[num[i][j]] + s[ num[i+dx[k]][j+dy[k]] ] > newMaxArea) { newx = i; newy = j; newMaxArea = s[num[i][j]] + s[ num[nowx][nowy] ]; newk = k; } } } int main() { freopen("castle.in", "r", stdin); freopen("castle.out", "w", stdout); scanf("%d%d", &m, &n); for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) scanf("%d", &map[i][j]); /*number the adjusted-parts and get the separate area*/ int number = 0, maxarea = 0; s[0] = - 1 << 20; // let the area outside the border be the minimal area memset(num, 0, sizeof(num)); for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) if (!num[i][j]) { number++; num[i][j] = number; qx[0] = i; qy[0] = j; int first = 0, last = 0; while (first <= last) { for (int k = 0; k < 4; k++) if ( !(map[qx[first]][qy[first]] & (1 << k)) ) { int nowx = qx[first] + dx[k], nowy = qy[first] + dy[k]; if (!num[nowx][nowy]) { qx[++last] = nowx; qy[last] = nowy; num[nowx][nowy] = number; } } first++; } s[number] = last + 1; if (s[number] > maxarea) maxarea = s[number]; } /*get the removed wall*/ for (int j = 1; j <= m; j++) {// wall first farthest to the west calc(j,1); calc(j,2); } /*print the answer*/ printf("%d/n%d/n%d/n%d %d %c/n", number, maxarea, newMaxArea, newx, newy, newk == 1 ? 'N':'E'); return 0; }

小問題：maxarea >?= s[number] 在USACO上編譯不能通過？

 

第二題：Ordered Fractions

題目意思很明白，第一反映排序，而且NLOGN排序可以過，於是沒想其他方法，就上了。

看了Analysis後，Russ有一種遞歸直接生成解的方法，很厲害啊。

We notice that we can start with 0/1 and 1/1 as our ``endpoints'' and recursively generate the middle points by adding numerators and denominators. 0/1 1/1 1/2 1/3 2/3 1/4 2/5 3/5 3/4 1/5 2/7 3/8 3/7 4/7 5/8 5/7 4/5 Each fraction is created from the one up to its right and the one up to its left. This idea lends itself easily to a recursion that we cut off when we go too deep. #include <stdio.h> #include <stdlib.h> #include <string.h> #include <assert.h> int n; FILE *fout; /* print the fractions of denominator <= n between n1/d1 and n2/d2 */ void genfrac(int n1, int d1, int n2, int d2) { if(d1+d2 > n) /* cut off recursion */ return; genfrac(n1,d1, n1+n2,d1+d2); fprintf(fout, "%d/%d/n", n1+n2, d1+d2); genfrac(n1+n2,d1+d2, n2,d2); } void main(void) { FILE *fin; fin = fopen("frac1.in", "r"); fout = fopen("frac1.out", "w"); assert(fin != NULL && fout != NULL); fscanf(fin, "%d", &n); fprintf(fout, "0/1/n"); genfrac(0,1, 1,1); fprintf(fout, "1/1/n"); }

 

 第三題：Sorting A Three-Valued Sequence

這道題很水，先兩個儘量交換，剩下的只能三個交換了，由於最終的排列是給定的，所以只需O(N)的時間

#include <stdio.h> #include <string.h> const int maxn = 1000; int num[4] = {0,0,0,0}, s[4][4]; int a[maxn]; void calc(int start, int end, int sym) { for (int i = start; i <= end; i++) s[sym][a[i]]++; } int min(int a, int b) { if (a > b) return b; else return a; } int main() { freopen("sort3.in", "r", stdin); freopen("sort3.out", "w", stdout); int n; scanf("%d", &n); for (int i = 1; i <= n; i++) { scanf("%d", &a[i]); num[a[i]]++; } memset(s,0,sizeof(0)); for (int i = 1; i <= 3; i++) { num[i] = num[i]+num[i-1]; calc(num[i-1]+1,num[i],i); } int tot = 0; for (int i = 1; i <= 3; i++) for (int j = 1; j <= 3; j++) if (i != j) { int Min = min(s[i][j], s[j][i]); tot += Min; s[i][j] -= Min; s[j][i] -= Min; } tot += (s[1][2]+s[1][3]) << 1; printf("%d/n", tot); return 0; }

 

第四題：Healthy Holstein

最多也只有15種飼料，0/1枚舉也就32768種，能承受，如果用二進制來構造方案統計的效率比較低，所以還是用DFS了來構造選擇方案比較划算

小技巧：記錄選擇的方案不需要一個數組，只需一個二進制數來記錄，方便，快速，給力！

#include <stdio.h> #include <string.h> const int maxneed = 25+3, maxtype = 15+3; int minchoose = 1 << 20, path, choose = 0, minpath, n, type; int map[maxtype][maxneed], get[maxneed], need[maxneed]; bool check() { for (int i = 0; i < n; i++) if (get[i] < need[i]) return false; return true; } void dfs(int now) { if (minchoose < choose) return; //小小的剪枝 if (now == type) { if (choose < minchoose && check()) { minpath = path; minchoose = choose; } return; } //choose path = (path << 1)+1; choose++; for (int i = 0; i < n; i++) get[i] += map[now][i]; dfs(now+1); for (int i = 0; i < n; i++) get[i] -= map[now][i]; choose--; path = path >> 1; //not choose path = path << 1; dfs(now+1); path = path >> 1; } int main() { freopen("holstein.in", "r", stdin); freopen("holstein.out", "w", stdout); scanf("%d", &n); for (int i = 0; i < n; i++) scanf("%d", &need[i]); scanf("%d", &type); for (int i = 0; i < type; i++) for (int j = 0; j < n; j++) scanf("%d", &map[i][j]); path = 0; memset(get,0,sizeof(get)); dfs(0); printf("%d ", minchoose); for (int i = type-1; i >= 0; i--) if (minpath & (1 << i)) printf("%d%c", type-i, --minchoose?' ':'/n'); return 0; }

 

第五題：Hamming Codes

要用到位運算來獲取第N位二進制數的操作；另外可以先預處理出所有數對之間的距離，我想數據本身就不太大，不預處理了，還浪費空間，時間還是挺可以的

#include <stdio.h> #include <string.h> const int max = 256, maxN = 64+3; short int dis[max][max]; int N, B, D; int path[maxN]; bool check(int a, int b) { int dif = 0; for (int i = 0; i < B; i++) if ((a & (1 << i)) != (b & (1 << i))) dif++; if (dif >= D) return true; else return false; } int dfs(int now) { if (now == N) { for (int i = 0; i < N; i++) printf("%d%c", path[i], i % 10 == 9 || i == N-1?'/n':' '); return 1; } int temp = path[now-1]; while (true) { temp++; bool flag = true; for (int i = 0; i < now; i++) if (!check(path[i], temp)) { flag = false; break; } if (flag) { path[now] = temp; if (dfs(now+1)) return 1; } if (temp > (1 << B)-1) { printf("Impossible"); return 1; } } } int main() { freopen("hamming.in", "r", stdin); freopen("hamming.out", "w", stdout); scanf("%d%d%d", &N, &B, &D); memset(path,0,sizeof(path)); path[0] = 0; path[1] = (1 << D) - 1; dfs(2); return 0; }

 

 

好久閒着了，今天下午+晚上間歇得把2.1寫完了，還是有收穫的，明天晚上就要動身去學校了，短暫的寒假結束了，準備迎接新一輪的挑戰吧！