计算一个 32 位无符号整数有多少个位为1

Counting out the bits

    可以很容易的判断一个数是不是2的幂次:清除最低的1位(见上面)并且检查结果是不是0.尽管如此,有的时候需要直到有多少个被设置了,这就相对有点难度 了。

GCC有一个叫做__builtin_popcount的内建函数,它可以精确的计算1的个数。尽管如此,不同于__builtin_ctz,它并没有被 翻译成一个硬件指令(至少在x86上不是)。相反的,它使用一张类似上面提到的基于表的方法来进行位搜索。这无疑很高效并且非常方便。

其他语言的使用者没有这个选项(尽管他们可以重新实现这个算法)。如果一个数只有很少的1的位,另外一个方法是重复的获取最低的1位,并且清除它。

查看引用来源

为了在 VC 上实现 __builtin_popcount (unsigned u) 的功能,自己写了两个函数,分别是 popcnt (unsigned u), popcount (unsigned u) 。

前者是通过清除 u 最低的 bit 1 ,直至 u 为 0 ,每次都为计数器加 1 。时间复杂度为 O (m) , m 为 bit 1 的个数。

后者是使用二分法,比较巧妙,跟踪调试一下就知道原理了。时间复杂度为 O (lg N) , N 为位数。

不过最高效的还是使用查表的方式来计算。

但是需要弄一个很大的表,不然随着位数的增长,查表的速度还是比不上二分法的速度。例如 64 位整数,保存 所有 8 位整数的结果 (256 个)。查表需要操作 8 次,而二分法需要 6 次而已。

测试代码如下:

// 计算一个 32 位无符号整数有多少个位为 1

// By Rappizit@yahoo.com.cn

#include <iostream>

#include <ctime>

using namespace std;

unsigned popcnt (unsigned u)

{

    unsigned ret = 0;

    while (u)

        {

        u = (u & (u - 1));    // 将 u 最右边的 1 清除

        ret ++;

        }

    return ret;

}

unsigned popcount (unsigned u)

{

    u = (u & 0x55555555) + ((u >> 1) & 0x55555555);

    u = (u & 0x33333333) + ((u >> 2) & 0x33333333);

    u = (u & 0x0F0F0F0F) + ((u >> 4) & 0x0F0F0F0F);

    u = (u & 0x00FF00FF) + ((u >> 8) & 0x00FF00FF);

    u = (u & 0x0000FFFF) + ((u >> 16) & 0x0000FFFF);

    return u;

}

int main ()

{

     // 先测试一下函数是否正常工作

    for (int i = 0; i <= 1000; i ++)

        {

        int k = popcnt (i);

        cout << k << " ";

        }

    cout << endl;

    

    for (int i = 0; i <= 1000; i ++)

        {

        int k = popcount (i);

        cout << k << " ";

        }

    cout << endl;

    

    // 比较速度     

    int test;

    test = (1U << 31) - 1;

    //test = 1000000;

    clock_t start, finish;

start = clock ();

    for (int i = 0; i < test; i ++)

        {

        popcnt (i);

        }

    finish = clock ();

    cout << (double)(finish - start) / (double)(CLOCKS_PER_SEC) << " s" << endl;

start = clock ();

    for (int i = 0; i < test; i ++)

        {

        popcount (i);

        }

    finish = clock ();

    cout << (double)(finish - start) / (double)(CLOCKS_PER_SEC) << " s" << endl;

start = clock ();

    for (int i = 0; i < test; i ++)

        {

        //__builtin_popcount (i);    // G++内建函数, Dev-C++ 可支持

        }

    finish = clock ();

    cout << (double)(finish - start) / (double)(CLOCKS_PER_SEC) << " s" << endl;

    

    system ("pause");

return 0;

}

/*

使用 VC 2005 ,在 test = (1U << 31) - 1 时, release 状态是 0 s 出结果的,但是 debug 状态就很慢, 换 test = 1000000 试试。

使用其他 IDE, 即使是 release 状态下,速度也不理想。但是可以测得 __builtin_popcount (通过查表来计算)比较快。

*/

说到底,这个函数到底有什么实际用处呢?当然有了,使用一个二进制数字表示一个集合的时候,枚举一个组合(子集),需要判断这个数字里面的 1 的个数是不是和子集的大小相等。这种方法通常是属于暴力法,如果不是使用二进制数字表示集合,很可能就计算超时了。

下面有个例子(sicily 1158),期待更好的解法,但是暴力法的效果还不错^_^而且实现简单。。

====================================帅气的分割 线====================================

Pick numbers

Total Submit : 371    Accepted Submit : 135

Problem

Given a matrix of size M*N, the elements of which are integer numbers from -10 to 10. You task is to go from 

the top-left corner (1, 1) to the bottom-right corner (M, N). You can only move right or down, and you can 

not go out of the matrix. The number in the grid you passed must be picked. The sum of numbers you picked must 

be positive and as minimal as possible.

Input

The input contains several test cases.

The first line of each test case are two integer numbers M, N (2<=M<=10, 2<=N<=10), indicating the number 

      of rows and columns of the matrix. The following M lines, each contains N integer numbers.

Output

For each test case, output the sum of numbers you picked on a single line. If you can't get a positive sum, 

output -1.

Sample input

2 2

0 2

1 0

3 3

0 0 0

0 0 0

0 0 0

Sample output

1

-1

Problem Source

2006 Algorithm Course Examination

====================================无聊的分割 线====================================

分析:

用一个二进制数字 path 表示所走的路径,从低位开始,如果某位为 1 ,表示向下走,否则向右走。

path 中的位 1 必须有 m - 1 个,使用 G++ 内建函数 __builtin_popcount (path) 来计算。

path 的起始值为 first = (1 << (m - 1)) - 1 ,表示先向下走 m - 1 步。

path 的终止值为 last = first << (n - 1) ,表示先向右走 n - 1 步。

从 first 到 last ,枚举所有可行的 path ,然后计算对应的 sum ,目标是找一个最小的 sum > 0 ,

如果找不到,sum = -1 。

注意此题不能用动态规划来求最小值,因为该最小值可能是负数。

枚举的次数是 2m + n - 2 - 2n - 1, 当 m, n <= 10 时为 512 。

#include <cstdio>

using namespace std;

const int M = 10;

const int N = 10;

int d [M] [N];

int m, n;

void f ()

...{

    for (int i = 0; i < m; i ++)

        ...{

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

            ...{

            scanf ("%d", &d [i] [j]);

            }

        }

    int first = (1 << (m - 1)) - 1;

    int last = first << (n - 1);

    int cnt = m - 1;

    int MASK = 1 << (m + n - 2);

    int max = (1U << 31) - 1;

for (int path = first; path <= last; path ++)

        ...{

        if (__builtin_popcount (path) == cnt)

            ...{

            int sum = d [0] [0];

            for (int i = 0, j = 0, mask = 1; mask < MASK; mask <<= 1)

                ...{

                if (mask & path)

                    ...{

                    i ++;                    

                    } else    ...{

                        j ++;

                        }

                sum += d [i] [j];

                }

            if (sum > 0 && sum < max)

                ...{

                max = sum;

                }

            }

        }

    printf ("%d ", (max == ((1U << 31) - 1) ? -1 : max));

}

int main ()

...{

    while (scanf ("%d%d", &m, &n) != EOF)

        ...{

        f ();

        }

    return 0;

}

提交结果:

Run ID      User Name      Problem      Language      Status      Run Time      Run Memory      Submit Time

84848       rappizit     1158     C++     Accepted     0.01 sec     260 KB     2007-09-17 13:38:57

===================================遗忘的分割 线=====================================

以下内容为新添加的:

经测试,该函数的速度是 G++ 的一半,所以有理由认为 G++ 保存的表的大小是 65536 的(或者使用汇编)。但是不可能在自己的代码里面加上那么一大段数据吧?权衡一下,表的大小设为 256 是适当的,可以在 VC 里面使用以下代码。

char poptable [256] = 

{

0, 1, 1, 2, 1, 2, 2, 3, 1, 2, 2, 3, 2, 3, 3, 4,

1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,

1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,

1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,

2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,

3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,

3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,

4, 5, 5, 6, 5, 6, 6, 7, 5, 6, 6, 7, 6, 7, 7, 8

};

unsigned __builtin_popcount (unsigned u)

{

    return poptable [u & 0xFF] + poptable [(u >> 8) & 0xFF] + poptable [(u >> 16) & 0xFF] + poptable [(u >> 24) & 0xFF];

}=====================20080725添加========================

unsigned popcount (unsigned u)

{

    u = (u & 0x55555555) + ((u >> 1) & 0x55555555);

    u = (u & 0x33333333) + ((u >> 2) & 0x33333333);

    u = (u & 0x0F0F0F0F) + ((u >> 4) & 0x0F0F0F0F);

    u = (u & 0x00FF00FF) + ((u >> 8) & 0x00FF00FF);

    u = (u & 0x0000FFFF) + ((u >> 16) & 0x0000FFFF);

    return u;

}

这个二分法的原理:

用8位二进制数来做示范好了,例如 u = 10110011。

10110011

00010001    //每两位取1位,即取偶数位, u & 01010101

01010001    //取奇数位并右移一位, (u >> 1) & 01010101

---------------

01100010    //上面两数相加,赋值给u,注意每两列相加的结果不会进位到第三列

00100010    //每四位取低两位, u & 00110011

00010000    //每四位取高两位并右移两位, (u >> 2) & 00110011

---------------

00110010    //上面两数相加,赋值给u

00000010    //每八位取低四位, u & 00001111

00000011    //每八位取高四位并右移四位,(u >> 4) & 00001111

---------------

00000101    //上面两数相加,赋值给u

最终结果 u = 5。

__builtin_popcount()的更多相关文章

  1. __builtin_popcount(n)

    Gcc提供的内建函数__builtin_popcount(n),可以精确计算n表示成二进制时有多少个1.借助这个函数可以快速判断一个数是否是2的幂. bool isPowerOfTwo(int n) ...

  2. Codeforces 839E Mother of Dragons【__builtin_popcount()的使用】

    E. Mother of Dragons time limit per test:2 seconds memory limit per test:256 megabytes input:standar ...

  3. C/C++中__builtin_popcount()的使用及原理

    __builtin_popcount()用于计算一个 32 位无符号整数有多少个位为1 Counting out the bits     可以很容易的判断一个数是不是2的幂次:清除最低的1位(见上面 ...

  4. __builtin_popcount() 函数

    详解 该函数的主要作用是计算一个数字的二进制中有多少个1,返回值就是其中1的个数. 它使用一张基于表的方法来进行位搜索,因此这个操作的执行效率很高 此处举一题 P1582 倒水 #include &l ...

  5. 洛谷P1582 倒水 二进制 lowbit __builtin_popcount

    P1582 倒水:https://www.luogu.org/problemnew/show/P1582 题意: 给定n瓶装有1升的水瓶,每次可以把两瓶装水量相同的水和成一瓶,问最少还要增加几瓶装有1 ...

  6. USACO 2.1 海明码 Hamming Codes (模拟+位运算+黑科技__builtin_popcount(n))

    题目描述 给出 N,B 和 D,要求找出 N 个由0或1组成的编码(1 <= N <= 64),每个编码有 B 位(1 <= B <= 8),使得两两编码之间至少有 D 个单位 ...

  7. GNU的__builtin_popcount函数

    用来计算32位的unsigned int中的1的个数, 其内部实现是根据查表法来计算的.

  8. POJ2104 K-th Number [分块做法]

    传送:主席树做法http://www.cnblogs.com/candy99/p/6160704.html 做那倒带修改的主席树时就发现分块可以做,然后就试了试 思想和教主的魔法差不多,只不过那个是求 ...

  9. [LeetCode] Counting Bits 计数位

    Given a non negative integer number num. For every numbers i in the range 0 ≤ i ≤ num calculate the ...

随机推荐

  1. Android 第三方授权(微信篇)

    0.申请开发者: https://open.weixin.qq.com/cgi-bin/frame?t=home/app_tmpl&lang=zh_CN 1.下载sdk包: https://o ...

  2. sqlserver关于事务

    事务的四个特性: 原子性,一致性,持久性,隔离性 原子性: 原子性:表示事务执行是作为原子,不可分割,整个语句要么执行,要么不执行sqlserver中每一个单独的语句可以看做是包含在事务中每一句本身具 ...

  3. 第三篇、微信小程序-网络请求API

    wx.request(OBJECT)发起的是https请求.一个微信小程序,同时只能有5个网络请求连接. OBJECT参数说明: 效果图: net.js Page({ data:{ result:{} ...

  4. Aisen仿新浪微博客户端项目源码

    新浪目前已经限制了第三方微博的很多API接口,加上平常时间不够,所以后续可能不会面向产品的去维护Aisen,不过也有了一些新的方向,例如引入最新Android-support-library,在一个完 ...

  5. C# WinForm TextBox添加水印效果

    1.新建项目添加WatermarkTextBox类: using System; using System.Collections.Generic; using System.Text; using ...

  6. 【知识分享】UIButton setTitle 设置为空 失效

    今天开发练习超级猜图,但是碰到了一个奇怪的问题 困扰我一个晚上,低效的夜晚 可恨~ 示例说明1 [button setTitle:@"" forState:UIControlSta ...

  7. css 笔记——设置禁用中文输入法

    ime-mode ime-mode的语法解释如下: ime-mode : auto | active | inactive | disabled 取值: auto  : 默认值.不影响IME的状态.与 ...

  8. C++对象的JSON序列化与反序列化探索

    一:背景 作为一名C++开发人员,我一直很期待能够像C#与JAVA那样,可以轻松的进行对象的序列化与反序列化,但到目前为止,尚未找到相对完美的解决方案. 本文旨在抛砖引玉,期待有更好的解决方案:同时向 ...

  9. Linux ---> 简单socket

    #include <stdio.h> #include <stdlib.h> #include <errno.h> #include <string.h> ...

  10. [DevExpress]GridControl 列头绘制Checkbox

    关键代码: /// <summary> /// 为列头绘制CheckBox /// </summary> /// <param name="view" ...