2、穷举（枚举、暴力）

一、算法思想

    列举出所有可能的情况，逐个判断哪些符合问题的要求，从而得到问题的全部解答。

    穷举法的关键问题：确定穷举范围、确定判断条件。

二、入门

    1、百鸡问题

    鸡翁一值钱五，鸡母一值钱三，鸡雏三值钱一。百钱买百鸡，问鸡翁、鸡母、鸡雏各几何？

#include<iostream> using namespace std; int main( ) { int i,j,k; for(i=1;i<=20;i++) //鸡翁 for(j=1;j<=33;j++) //鸡母、 穷举范围 { k=100-i-j; //鸡雏 if(k%3==0 && i*5+j*3+k/3==100) //判断条件 cout<<"鸡翁："<<i<<" 鸡母："<<j <<" 鸡雏："<<k<<endl; } return 0; }

三、例题

    1、生理周期

    http://poj.org/problem?id=1006

    思路：

    (1)输入四个整数：d1,d2,d3和d。p,e,i分别表示第一次体力、情感和智力高峰出现的时间（即距当年第一天的天数，如：0表示当年第一天，1表示当年第2天……）。d是给定的时间，可能小于d1,d2,d3。

#include <stdio.h> int main() { int x, tcase, d, d1, d2, d3; tcase = 0; //FILE *p=fopen("1.txt","r"); while (1) { scanf("%d%d%d%d", &d1, &d2, &d3, &d); if (d1 == -1 && d2 == -1 && d3 == -1 && d == -1) break; for (x = d+1; x <= 21252; ++x) if ((x-d1) % 23 == 0) //找到23的第一个周期 break; for (; x <= 21252; x += 23)//找23和28相同的周期 if ((x-d2) % 28 == 0) break; for (; x <= 21252; x += 23*28)//找33和23、28相同的周期 if ((x-d3) % 33 == 0) break; printf("Case %d: the next triple peak occurs in %d days./n", ++tcase, x-d); } return 0; }

    2、假币问题

    http://poj.org/problem?id=1013

    思路：

    (1)12枚硬币real[12],标志着如果是even时，它们都标志为true

    (2)flag[12]当up时，左边的++，右边的--

    (3)当down时，左边的--，右边的++

    (4)当real=true的时候必定为真币，其他的偏离值越大，它就是假的

#include<iostream> #include<cmath> #include<string> using namespace std; int main() { int i,n,j; bool real[12]; int flag[12]; string a,b,c; int la,lb; int max,p,absmax; cin>>n; while(n--) { memset(real,false,sizeof(real)); memset(flag,0,sizeof(flag)); for(j=0;j<3;j++) { cin>>a>>b>>c; la=a.length(); lb=b.length(); switch(c[0]) { case 'e': for(i=0;i<la;i++) real[a[i]-'A']=true; for(i=0;i<lb;i++) real[b[i]-'A']=true; break; case 'u': for(i=0;i<la;i++) flag[a[i]-'A']++; for(i=0;i<lb;i++) flag[b[i]-'A']--; break; case 'd': for(i=0;i<la;i++) flag[a[i]-'A']--; for(i=0;i<lb;i++) flag[b[i]-'A']++; break; } } max=0; p=0; for(i=0;i<12;i++) { if(!real[i]) { absmax=abs(flag[i]); if(absmax>max) { max=absmax; p=i; } } } printf("%c is the counterfeit coin and it is ",'A'+p); if(flag[p]>0) printf("heavy./n"); else printf("light./n"); } return 0; }

    3、完美立方

    http://poj.org/problem?id=1543

    思路：直接暴力吧

    #include<iostream> using namespace std; int main() { int a,b,c,d,n; int temp,temp2; cin>>n; for(a=2;a<=n;a++) { for(b=2;b<=n;b++) { for(c=b+1;c<=n;c++) { for(d=c+1;d<=n;d++) { if(a*a*a==b*b*b+c*c*c+d*d*d) { cout<<"Cube = "<<a<<", Triple = ("<<b<<","<<c<<","<<d<<")"<<endl; } } } } } return 0; }

    4、画家问题

    http://poj.org/problem?id=1681

    思路：

    (1)我们按照从上到下、从左到右的顺序确定需要粉刷的砖块。

    (2)考虑相邻两行间的相互影响，下一行的每一格只能影响上一行正对的一格。因此，当上一行的粉刷方案确定后，下一行的粉刷方案也随之确定。

例如：

   

    第一行粉刷后墙壁状态如上左图显示，则粉刷第二行时，为了使第一行都变成黄色，必须粉刷第1、3、4、7格，粉刷后如上右图所示。

    总的算法是枚举第一行的粉刷方案，然后用上述方法依次确定下面各行的粉刷方案。在由于第k行的粉刷方案保证第k-1行为黄色。n行全部完成后，只需检查第n行是否全为黄色，就可以判断这种方案是否合法。

    用一个bool数组map[n][n]表示整个墙壁的状态，map[i][j]=true表示第i行，j列为白色

     代码：

#include<iostream> using namespace std; char board[16][16]; bool map[16][16]; int cnt;//涂块的次数 int t,n,i,temp,j; void click(int x,int y) { ++cnt; map[x][y]=!map[x][y]; if(x-1>=0) map[x-1][y]=!map[x-1][y]; if(y-1>=0) map[x][y-1]=!map[x][y-1]; if(x+1<n) map[x+1][y]=!map[x+1][y]; if(y+1<n) map[x][y+1]=!map[x][y+1]; } int check(int s) { cnt=0; int i; for(i=0;i!=n;++i) //初始状态 for(j=0;j!=n;++j) map[i][j]=board[i][j]=='w'; for(i=0;i!=n;++i) //确定第一行的具体情况 if((s&(1<<i))!=0) click(0,i); for(i=1;i!=n;++i) //涂第二行至最后一行 for(j=0;j!=n;++j) if(map[i-1][j]) click(i,j); for(j=0;j!=n;++j)//检查最后一行是否全是黄色 if(map[n-1][j]) return 999999999; return cnt; } int main() { cin>>t; while(t--) { cin>>n; for(i=0;i<n;i++) cin>>board[i]; int ans=999999999; for(i=0;i!=(1<<n);++i)//第一行有2^n种变化，对每一种变化尝试 { temp=check(i); ans=temp<ans?temp:ans; } if(ans<999999999) cout<<ans<<endl; else cout<<"inf"<<endl; } return 0; return 0; }