二进制是计算技术中广泛采用的一种数制。二进制数据是用0和1两个数码来表示的数。利用二进制的特点,可以用于枚举一个集合中各元素的所有组合情况。

例如,设某集合M中有3个元素A、B和C,即M={A,B,C}。可以用3位二进制数来枚举3个元素的各种组合情况(也可以称为子集),每一位二进制数字对应一个元素,1表示选中,0表示不选。

二进制数 001,表示选中A元素,B和C元素不选择,即子集为 { A };

二进制数 010,表示选中B元素,A和C元素不选择,即子集为 { B };

二进制数 011,表示选中A和B元素,C元素不选择,即子集为 { A,B };

二进制数 100,表示选中C元素,A和B元素不选择,即子集为 { C };

二进制数 101,表示选中A和C元素,B元素不选择,即子集为 { A,C };

二进制数 110,表示选中B和C元素,A元素不选择,即子集为 { B,C };

二进制数 111,表示选中A、B和C这3个元素均选中,即子集为 { A,B,C };

当然,二进制数000表示3个元素均不选择,就是子集是空集。

这样,用循环 for (int i = 0; i < (1 << n); i++) 就可以遍历每一个子集,即遍历各种组合情况。其中 (1 << n)表示 2n ,2n-1就对应n为二进制数,各位上的数字均为1。

然后,可以利用二进制的位运算,来判断二进制数的某一位上是0还是1。

二进制位运算的左移运算 a << b就表示把a转为二进制后左移b位(即在后面添b个0)。

这样, 1<<0=1(0) ; 1<<1=2 (10); 1<<2=4 (100); 1<<3=8 (1000);

……

因此,要判断一个二进制数x的第k位上是0还是1,只需计算 x & (1<<K),若计算出来的值为0,则x的第k位为0;若计算出来的值非0,则x的第k位为0。

这样,再利用循环

For (int j = 0; j < n; j++)

if (i & (1 << j)) { … }

就可以对每一个子集中所选取的具体元素进行处理。

编写如下的示例程序。

#include int main(){    char s[6]="ABCDE";    int n;    scanf("%d",&n);    int i,j;    for (i=0; i<(1<<n); i++)    // 枚举从0~2^n-1中组合情况    {        printf("[  ");        for (j = 0; j < n; j++) // 遍历二进制的每一位        {            if (i & (1 << j))   // 判断二进制第j位是否为1,若为1,表示元素选中            {                printf("%c  ",s[j]); // 如果第j个元素选中,则输出            }        }        printf("]\n");    }    return 0;}

该程序运行后,输入3,可以输出A、B、C这3个元素的8中组合情况。运行结果如下所示。

下面通过一些实例来看二进制枚举的应用方法。

【例1】最接近的和

问题描述

给出两个整数n(n≤20)和b,然后给出n个整数h1、h2、…、hn。问如何在这n个整数里面选择若干个数使得它们的和在满足大于等于b的情况下最小,并输出这个和与b的差。

输入

第1行:2个用空格隔开的整数:N和B。

第2…N+1行:第i+1行是1个整数:Hi。

输出

输出1个非负整数,即满足大于等于b的情况下的最小和值与b的差。

输入样例

5 16

3

1

3

5

6

输出样例

1

(1)编程思路。

直接套用上面的二进制枚举的模板对h1、h2、…、hn这n个元素的各种组合情况进行处理。只是在具体处理时,将每次挑选中的元素累加起来。当和值大于等于b时,记录最小的和值即可。

(2)源程序。

#include int main(){    int n,b;    scanf("%d%d",&n,&b);    int i,j,h[21];    for (i=0;i<n;i++)        scanf("%d",&h[i]);    int ans=0x3f3f3f3f;    for (i=0; i<(1<<n); i++)    // 枚举从0~2^n-1中组合情况    {        int sum=0;        for (j = 0; j < n; j++) // 遍历二进制的每一位        {            if (i & (1 << j))   // 判断二进制第j位是否为1,若为1,表示元素选中            {                sum+=h[j];      // 如果第j个元素选中,则加上它            }        }        if (sum>=b)        {            if (sum<ans) ans=sum;        }    }    printf("%d\n",ans-b);    return 0;}

将上面的源程序提交给 洛谷题库P2677 [USACO07DEC]Bookshelf 2 B(https://www.luogu.com.cn/problem/P2677),可以Accepted。

【例2】权利指数

问题描述

在选举问题中,总共有n个小团体,每个小团体拥有一定数量的选票数。如果其中m个小团体的票数和超过总票数的一半,则此组合为“获胜联盟”。n个团体可形成若干个获胜联盟。一个小团体要成为一个“关键加入者”的条件是:在其所在的获胜联盟中,如果缺少了这个小团体的加入,则此联盟不能成为获胜联盟。一个小团体的权利指数是指:一个小团体在所有获胜联盟中成为“关键加入者”的次数。请你计算每个小团体的权利指数。

输入

输入数据的第一行为一个正整数T,表示有T组测试数据。每一组测试数据的第一行为一个正整数n(0<n<=20)。第二行有n个正整数,分别表示1到n号小团体的票数。

输出

对每组测试数据,在同一个行按顺序输出1到n号小团体的权利指数。

输入样例

2

1

10

7

5 7 4 8 6 7 5

输出样例

1

16 22 16 24 20 22 16

(1)编程思路。

同样直接套用上面的二进制枚举的模板对n个小团体的各种组合情况进行处理。若某个团体组合的得票数超过了总票数的一半,则依次对该组合中的每个小团队进行判断,若去掉了某个小团体后,得票数小于总票数的一半,则去掉的小团体就是关键加入者,进行计数。

(2)源程序。

#include #include <string.h>int main(){    int t;    scanf("%d", &t);    while (t--)    {        int a[21],key[21], ans[21];        memset(ans, 0, sizeof(ans));        int n,total = 0;        scanf("%d", &n);        int i;        for (i = 0; i < n; i++)        {            scanf("%d", &a[i]);            total += a[i];        }        total /= 2;        for (i = 1; i < 1 << n; i++) // 二进制枚举各团体组合        {            int j, k = 0, sum = 0;            for (j = 0; j < n; j++)  // 遍历二进制的每一位            {                if (i & (1 << j))    // 判断二进制第j位是否为1                {                   key[k++] = j; sum += a[j];                }            }            if (sum > total)         // 如果团体的票数超过总票数的一半            {                for (j = 0; j < k; j++)                  if (sum - a[key[j]] <= total) // 编号为key[j]的元素为 关键加入者                      ans[key[j]]++;            }        }        printf("%d", ans[0]);        for (i = 1; i < n; i++)            printf(" %d", ans[i]);        printf("\n");    }    return 0;} 

将上面的源程序提交给 HDU题库 HDU 1557 权利指数(http://acm.hdu.edu.cn/showproblem.php?pid=1557),可以Accepted。

【例3】目标和习惯

问题描述

小明有很多习惯,比如看电影、听音乐、玩电脑游戏和足球等等。此外,他还有很多目标,因为不是每个习惯都有助于实现他的目标,他陷入了在目标和习惯之间做出选择的困境。

现在,我们将一个习惯对每个目标的影响定义为向量,向量的元素是整数,例如。向量(100,90,-10,80)表示它对目标1有100分的影响,对目标2有90分的影响、对目标3有-10分的影响和对目标4有80分的影响(正表示良好效果,负表示不良效果),每个目标的给定要求表示为整数。请帮助小明实现他的目标,并尽可能多地保持他的习惯。

输入

有多种情况,读取数据直到EOF。(不超过10例)

第1行:小明的目标数N(0<N<=20)

第2行:每个目标的要求(0<w<=1000)

第3行:小明的习惯数M(0<M<=16)

第4-M+4行:每行包含N个整数,第i个整数表示对第i个目标的影响(-1000<=data<=1000)。

输出

对于每种情况:输出是一行,其中包含:

小明可以保持的最大习惯数,其次是:

小明可以保持的习惯的排序列表(从最小到最大)。如果有一组以上的习惯可以满足要求,请选择习惯数最少的一组。

如果无法实现小明的目标,请输出0。

输入样例

4

100 200 300 400

3

100 100 400 500

100 -10 50 300

100 100 -50 -50

输出样例

2 1 3

(1)编程思路。

采用二进制枚举对m个习惯的全部组合情况进行处理。对每种习惯组合,将组合中的各习惯对目标的影响累加到数组g中。若数组g中每个目标均可实现,则是一种满足要求的习惯组合。

(2)源程序。

#include #include <string.h>int main(){    int goal[21];    int f[17][21],g[21];    int n;    while (scanf("%d",&n)!=EOF)    {        int ans1=-1,ans2;        int i,j,k;        for (i=0;i<n;i++)           scanf("%d",&goal[i]);        int m;        scanf("%d",&m);        for (i=0;i<m;i++)            for(j=0;j<n;j++)               scanf("%d",&f[i][j]);        for (i=0;i<(1<<m);i++)         // 对所有m个习惯的组合进行枚举        {            for (k=0;k<n;k++) g[k]=0;  // 各目标的达到值初始化为0            int num=0;                 // 当前组合中的习惯个数            for (j=0;j<m;j++)                if (i&(1<<j))          // 若组合中有第j个习惯                {                    num++;                    for (k=0;k<n;k++)  // 第j个习惯对每个目标k的影响                        g[k]+=f[j][k];                }            for (k=0;k<n;k++)                if (g[k]<goal[k]) break;   // 第k个目标不能实现            if (k==n)                  // 当前习惯组合中,所有目标可达成            {                if (num>ans1) { ans1=num; ans2=i; }            }        }        if (ans1==-1)            printf("0\n");        else        {            printf("%d",ans1);            for (i=0;i<n;i++)               if (ans2&(1<<i))                  printf(" %d",i+1);            printf("\n");        }    }    return 0;}

将上面的源程序提交给 HDU题库 HDU 4152 ZZY’s Dilemma(http://acm.hdu.edu.cn/showproblem.php?pid=4152),可以Accepted。

【例4】调查问卷

问题描述

度度熊为了完成毕业论文,需要收集一些数据来支撑他的论据,于是设计了一份包含 m 个问题的调查问卷,每个问题只有 ‘A’ 和 ‘B’ 两种选项。

将问卷散发出去之后,度度熊收到了 n 份互不相同的问卷,在整理结果的时候,他发现可以只保留其中的一部分问题,使得这 n 份问卷仍然是互不相同的。这里认为两张问卷是不同的,当且仅当存在至少一个被保留的问题在这两份问卷中的回答不同。

现在度度熊想知道,存在多少个问题集合,使得这 n 份问卷在只保留这个集合的问题之后至少有 k 对问卷是不同的。

输入

第一行包含一个整数 T,表示有 T 组测试数据。

接下来依次描述 T 组测试数据。对于每组测试数据:

第一行包含三个整数 n,m 和 k,含义同题目描述。

接下来 n 行,每行包含一个长度为 m 的只包含 ‘A’ 和 ‘B’ 的字符串,表示这份问卷对每个问题的回答。

保证 1≤T≤100,1≤n≤103,1≤m≤10,1≤k≤106,给定的 n 份问卷互不相同。

输出

对于每组测试数据,输出一行信息 “Case #x: y”(不含引号),其中 x 表示这是第 x 组测试数据,y 表示满足条件的问题集合的个数,行末不要有多余空格。

输入样例

2

2 2 1

AA

BB

2 2 2

AA

BB

输出样例

Case #1: 3

Case #2: 0

(1)编程思路。

先把每张问卷对10个问题回答的A、B字符串转换成二进制数保存在数组a中,转换时,A转换为1,B转换为0。

采用二进制枚举的方法依次枚举每个可能的子集(选择部分问题的情况,下面将枚举到的某个二进制数i称为当前状态),与每张问卷的这个子集情况(当前状态i)下的答案匹配。具体做法是:将当前状态i与每张试卷的答案(也已转换为了二进制数并保存在数据元素a[j]中)进行与运算,即now=i & a[j],并将cnt[now]加1。显然,若问卷j和问卷k的答案相同,则a[j]与a[k]的值一定相同,这样now的值也相同,cnt[now]的值会加上2。

当匹配到第j张问卷时,第j张问卷与前面的j-1张问卷在当前选择状态i下一定有cnt[i&a[j]]张(包括第j张)的答案相同,也就是有j- cnt[i&a[j]]张的答案存在不同,将它累加到和值sum上。当某个当前状态i与每张问卷匹配完成后,如果和值sum大于或等于k,这种当前状态 i 就是符合题意的一种子集选择。遍历所有可能的子集,就得到最终的答案。

(2)源程序。

#include #include <string.h>int main(){    int t,iCase=0;    scanf("%d",&t);    while (t--)    {        int n, m, k;        scanf("%d%d%d",&n,&m,&k);        int a[1001];        int i,j;        for (i=1; i<=n; i++)        {            char s[15];            scanf("%s",s);            int num=0;            for (j=0; s[j]!='\0'; j++)  // 第i张问卷对10个问题回答的A、B字符串转换成二进制数保存在数组a[i]            {                if (s[j]=='A')                    num=num*2+1;                else  num=num*2;            }            a[i] = num;        }        int ans = 0;        for (i=1; i<(1<<m); i++)       // 对m个问题的全部组合情况进行枚举        {            int cnt[1050];            memset(cnt, 0, sizeof(cnt));            int sum = 0;            for (j=1; j<=n; j++)            {                int now = i & a[j];                cnt[now]++;                sum += j-cnt[now];            }            if(sum>=k) ans++;        }        printf("Case #%d: %d\n",++iCase,ans);    }    return 0;}

将上面的源程序提交给 HDU题库 HDU 6344 调查问卷(http://acm.hdu.edu.cn/showproblem.php?pid=6344),可以Accepted。

【例5】连续求余

问题描述

有一个整数a和n个整数b1,…,bn。在按任意顺序从b1、…、bn中选择一些数字后,例如c1、…、cr,我们希望确保a mod c1 mod c2 mod…mod cr=0(即,a将成为每次除以ci的余数,最后,我们希望a变为0)。请确定r的最小值。如果无法实现目标,请输出−1代替。

输入

第一行包含一个整数T≤5,表示测试用例的数量。

对于每个测试用例,有两行:

第一行包含两个整数n和a(1≤n≤20,1≤a≤106)。

第二行包含n个整数b1,…,bn(1≤bi≤106

输出

每个测试用例输出1行,该行中输出所求得的最小的r值。

输入样例

2

2 9

2 7

2 9

6 7

输出样例

2

-1

(1)编程思路。

采用二进制枚举的方法对b1、…、bn中的数的各种选择组合进行处理。只是在本题处理时,需要先将数组b按从大到小的顺序进行排列。因为,若先选择一个小的数进行求余,则后面较大的数的选择毫无意义。​​

(2)源程序。

#include int main(){    int t;    scanf("%d",&t);    while (t--)    {        int n,a,b[25]={0};        scanf("%d%d",&n,&a);        int i,j;        for (i=0;i<n;i++)            scanf("%d",&b[i]);        for (i=0;i<n-1;i++)         // 将数组b按从大到小的顺序排列            for (j=0;j<n-1-i;j++)               if (b[j]<b[j+1])               {                  int temp;                  temp=b[j];  b[j]=b[j+1];  b[j+1]=temp;               }        int ans=25;    // 所以数全部选,不会超过20个数        for (i=0;i<(1<<n);++i)        {            int res=a,cnt=0;            for (j=0;j<n;++j)            {                if(i>>j&1)                {                    res=res%b[j];                    cnt++;                    if(res==0)break;                }            }            if (res==0)                ans=(cnt<ans)?cnt:ans;        }        if(ans==25) printf("-1\n");        else  printf("%d\n",ans);    }    return 0;}

将上面的源程序提交给 HDU题库 HDU 5339 Untitled(http://acm.hdu.edu.cn/showproblem.php?pid=5339),可以Accepted。

在通过上面的5道例题的学习,掌握了二进制枚举方法的应用的基础上,可以刷一下HDU题库中的下面几道题目,也均可采用二进制枚举的方法来完成。

【练习1】HDU 1882 Strange Billboard (http://acm.hdu.edu.cn/showproblem.php?pid=1882)。

#include #include <string.h>int main(){    int r,c;    while (scanf("%d%d",&r,&c)!=EOF)    {        if (r==0 && c==0) break;        int i,j,k;        char map[20][20],temp[20][20];        for (i=0;i<r;++i)            scanf("%s",map[i]);        if (r<c)   // 若行数小于列数,则使矩阵转置。        {            memcpy(temp,map,sizeof(map));            for (i=0;i<r;++i)            {                for (j=0;j<c;++j)                {                    map[j][i]=temp[i][j];                }            }            k=r; r=c; c=k;        }        int hash[20],key[20];        for(i=0;i<r;++i)    // 将第i行数据以二进制形式保存在数组元素hash[i]中,‘X’符号位置1        {            hash[i]=0;            for (j=0;j<c;++j)            {                if (map[i][j]=='X')                    hash[i]=(hash[i]<<1)^1;                else                    hash[i]=(hash[i]<<1);            }        }        int ans=1000;   // 按题意,最多256个方格,因此翻转次数肯定不会超过1000        for (i=0; i<(1<<c); i++)    // 对第0行的所有翻转组合进行二进制枚举        {            int cnt=0;            memcpy(key,hash,sizeof(hash));  // 将初始画面情况数组hash暂存到key数组中            for (j=0;j<c;j++)     // 处理第0行(首行)            {                if (i & (1<<(c-1-j)))                {                    key[0]^=(1<<(c-1-j));  // 当前方块翻转                    if (j-1>=0)   // 当前方块左边(前一列)的方块翻转                        key[0]^=(1<<(c-1-j+1));                    if (j+1<c)   // 当前方块右边(后一列)的方块翻转                        key[0]^=(1<<(c-1-j-1));                    if (1<r)     // 当前方块下方(下一行)的方块翻转                        key[1]^=(1<<(c-1-j));                    cnt++;       // 翻转次数加1                }            }            for (j=1;j<r;++j)         // 上1行确定了,下一行只能根据上一行进行处理。依次处理后面的r-1行            {                for (k=0;k<c;++k)     // 对当前行的每列的方块进行判断,根据正上方的方块状态看是否需要翻转                {                    if (key[j-1]&(1<<(c-1-k))) // 当前方块正上方方块为'X‘,需要翻转当前方块                    {                        cnt++;                        key[j]^=(1<<(c-1-k));  // 当前方块翻转                        if (k-1>=0)            // 当前方块左边(前一列)的方块翻转                            key[j]^=(1<<(c-1-k+1));                        if (k+1<c)             // 当前方块右边(后一列)的方块翻转                            key[j]^=(1<<(c-1-k-1));                        if (j-1>=0)            // 当前方块上方(上一行)的方块翻转                            key[j-1]^=(1<<(c-1-k));                        if (j+1<r)             // 当前方块下方(下一行)的方块翻转                            key[j+1]^=(1<<(c-1-k));                    }                }            }            if (!key[r-1])  // 由于下一行必受上一行约束,则只要判断最后行是否达到要求即可            {                if (ans>cnt) ans=cnt;            }        }        if (ans<1000)            printf("You have to tap %ld tiles.\n",ans);        else            printf("Damaged billboard.\n");    }    return 0;}

参考源程序

【练习2】HDU 3118 Arbiter (http://acm.hdu.edu.cn/showproblem.php?pid=3118)。

#include int main(){   int t;   scanf("%d",&t);   while (t--)   {       int n,m;       scanf("%d%d",&n,&m);       int map[305][2];       int i,j;       for (i=0;i<m;i++)       {           scanf("%d%d",&map[i][0],&map[i][1]);       }       int ans=0x3f3f3f3f;       for (i=0;i<(1<<n);i++)       {           int cnt=0;           for (j=0;j<m;j++)           {               if (((i>>map[j][0]) & 1)==((i>>map[j][1])&1))                   cnt++;           }           if (cnt<ans) ans=cnt;       }       printf("%d\n",ans);   }   return 0;}

参考源程序

【练习3】 HDU 4462 Scaring the Birds (http://acm.hdu.edu.cn/showproblem.php?pid=4462) 。

#include #include <string.h>int abs(int x){    return x>0?x:-x;}int  main(){    int n;    while (scanf("%d", &n) && n!=0)    {        int m;        scanf("%d", &m);        int i,j;        int r[51],c[51],f[51];        for (i = 0; i < m; i++)            scanf("%d%d", &r[i], &c[i]);        for (i = 0; i < m; i++)            scanf("%d", &f[i]);        int ans = 100;        for (i = 0; i < (1 << m); i++)   // 二进制枚举各位置放置稻草人的组合情况        {            int num = 0;                 // 放稻草人的个数            int a[10],vis[51][51];            memset(vis, 0, sizeof(vis));            for (j = 0; j < m; j++)                if (i&(1 << j))                    a[num++] = j;        // 保存所放稻草人的编号            int  flag = 1;            int x,y,z;            for (x = 1; x <= n; x++)     // 放置了num个稻草人后,看稻田所有位置是否全部覆盖            {                for (y = 1; y <= n; y++)                {                    for (z = 0; z < m; z++)                        if (r[z] == x && c[z] == y)                            vis[x][y] = 1;                    if (vis[x][y])                        continue;                    for (z = 0; z < num; z++)                    {                        if (abs(x - r[a[z]]) + abs(y - c[a[z]]) <= f[a[z]])                        {                            vis[x][y] = 1;                            break;                        }                    }                    if (vis[x][y] == 0)    // 当前位置(x,y)不能覆盖                        flag = 0;                    if (!flag)                        break;                }                if (!flag)                    break;            }            if (flag)                if (ans>num) ans=num;        }        if (ans == 100)            printf("-1\n");        else            printf("%d\n", ans);    }    return 0;}

参考源程序

【练习4】HDU 5616 Jam’s balance(http://acm.hdu.edu.cn/showproblem.php?pid=5616)。

#include #include <string.h>int main(){    int t;    scanf("%d",&t);    while (t--)    {        int vis[2005];        memset(vis,0,sizeof(vis));        int n;        scanf("%d",&n);        int i,j,a[21];        for (i=0; i<n; i++)            scanf("%d",&a[i]);        for (i=0; i<(1<<n); i++)   // 二进制有n位,枚举二进制数        {            int sum=0;             // 当前砝码组合可称出的总重量            for (j=0; j<n; j++)            {                if (i&(1<<j))                    sum+=a[j];            }            vis[sum]=1;            for (j=0; j<n; j++)    // 将每个砝码放在天平另一边            {                if (sum-a[j]>=0)                    vis[sum-a[j]]=1;            }        }        int m,x;        scanf("%d",&m);        while (m--)        {            scanf("%d",&x);            if (vis[x]) printf("YES\n");            else  printf("NO\n");        }    }    return 0;}

参考源程序

【练习5】HDU 6435 Problem J. CSGO (http://acm.hdu.edu.cn/showproblem.php?pid=6435)。

#include int main(){    int t;    scanf("%d",&t);    while (t--)    {        int n,m,K;        scanf("%d%d%d",&n,&m,&K);        long long mw[35],sw[35];    // K最大为5,2^5=32。性能组合情况不超过32种        int i,j,k;        for (i=0;i<32;i++)            mw[i] = sw[i] = -1e16;  // 为求最大值,主、副武器的每种性能组合情况均预置一个很小的数        int s,x[6];        long long ret;        for (i=0;i<n;i++)           // 依次处理每个主武器        {            scanf("%d",&s);            for (j=0;j<K;j++)                scanf("%d",&x[j]);            for (j=0;j<(1<<K);j++)   // 对当前主武器的性能组合进行二进制枚举            {                ret = s;                for (k=0;k<K;k++)                    ret += x[k]*((1&(j>>k))?1:-1);                if (mw[j]<ret) mw[j]=ret;            }        }        for (i=0;i<m;i++)            // 依次处理每个副武器        {            scanf("%d",&s);            for (j=0;j<K;j++)               scanf("%d",&x[j]);            for (j=0;j<(1<<K);j++)   // 对当前副武器的性能组合进行二进制枚举            {                ret = s;                for (k=0;k<K;k++)                    ret += x[k]*((1&(j>>k))?1:-1);                if (sw[j]<ret) sw[j]=ret;            }        }        long long ans = 0;        for (i=0;i<(1<<K);i++)        {            ret=mw[i]+sw[(1<<K)-1-i];  // 主、副武器二进制位要求互补            if (ans<ret) ans=ret;        }        printf("%lld\n",ans);    }    return 0;}

参考源程序