【C++】20191002 奇怪数二进制简便算法

数据结构：找出1~1000内所有奇怪数

【盈数】
盈数（丰数，过剩数，abundant number)是一种特殊的自然数，除去它本身以外的一切正约数的和大于它本身。

【奇怪数】
奇怪数是盈数，同时满足：
在其约数构成的所有子集中，找不到任何一个子集的和等于它本身

【思路要点】
用当前数的所有真因子之和减去其本身，缩小遍历范围

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int* Narrow(int num)//定义函数Narrow用于缩小查找范围
{
	int i, k, p = 2, t = 0, sum = 0, rest;//定义i,j,k作为循环变量，t,p分别用于记录数组a,b中的元素个数
	int a[1000], b[1000];

	for (i = 1; i < num; i++)//利用一个for循环求num的所有真因子，存入数组a中，求和，并用t记录真因子个数
	{
		if (num % i == 0)
		{
			a[t] = i;
			sum = sum + a[t];
			t = t + 1;
		}
	}
	rest = sum - num;//用所有真因子之和减去num本身，记作rest，缩小遍历范围
	b[0] = rest;//将rest值存入数组b中
	
	for (k = 0; k < t; k++)//利用一个for循环选出小于rest的真因子，存入数组b中，并用p记录b中元素个数
	{
		if (a[k] <= rest)
		{
			b[p] = a[k];
			p = p + 1;
		}
	}
	b[1] = p-2;//将真因子个数（p-1-1）存入数组b中

	return b;//返回数组b
}

int Select(int b[], int s, int p)//定义函数Select，利用二进制法找到由rest的真因子组成的集合的全部子集，并对每个子集求和
{
	int i, sum = 0;

	for (i = 0; i < p; i++)//利用for循环遍历二进制的每一位，其中p为数组b的元素个数
	{
		if (s & (1 << i))//判断s的二进制的第i位是否存在（是否为1）（二进制移位运算、与运算）（注意：二进制从第0位开始）
			sum = sum + b[i+2];//求和
	}

	return sum;//返回每个子集元素的和
}

int main()
{
	int i, j, p, total, rest, sum, flag;
	int *b;

	for (i = 71; i < 1000; i++)
	{
		b = Narrow(i);//缩小范围后的数组b
		rest = b[0];//目标值rest
		p = b[1];//数组b中真因子个数
		total = pow(2,p)-1;//数组b的真子集个数

		if (rest > 0)//当该数为盈数时
		{
			flag = 1;
			for (j = 0; j <= total; j++)//遍历b的每一个真子集
			{
				sum = Select(b, j, p);//得到每一个真子集的和
				if (sum == rest)
				{
					printf("%d is not an odd number\n", i);//若某个真子集的和等于rest，则此时的i不是奇怪数
					flag = 0;//若某个真子集的和不等于rest，改变flag的值
					break;
				}
			}
        
			if (flag == 1)
				printf("%d is an odd number\n", i);//若flag的值未被改变，说明b的所有真子集的和都不等于rest，此时的i即为所求的奇怪数
		}
		    else
				printf("%d is not an odd number\n", i);	//若该数不是盈数，则该数必然不是奇怪数
	}

	return 0;
}

【编程总结】

1. 如何找到数num的真因子？
   （1）思路一：

for (i=1; i<num; i++)//遍历1~num所有数
{
   If (num % i== 0)//若num整除i
      {……}//将i存入数组
}

（2）思路二：（更简便）

 for (i = sqrt(num); i > 0; i--)//遍历1~num的开平方数
 {
     if (num % i == 0)//若num整除i
	{
		b = k / i;//i与b都是num的真因子
		if (b ! = i)//若i与b不相等
		{……}//将i与b都存入数组
		else {……}//若i与b相等，则只将i存入数组
	 }
 }

2. 如何将找到的真因子存起来？
   （1）利用数组：多加一个变量用于记录真因子个数，即：

int a[1000];
for(i....)
{
   a[t]=i;  
   t++; 
}

（2）利用向量：

 vector<int> s;//定义一个向量s
 s.insert(s.begin(), i)//将找到的真因子插入向量前端

3. （1）如何在自定义函数中返回一个数组？
   将函数定义为指针形式：

int*ABC(...)
{ 
  int a[1000]; 
  ...... ; 
  return a;
}

（2）自定义函数中返回的数组如何在主函数中进行调用？
在主函数中将数组定义为指针形式：

 int*a；
 a=ABC(...);

（3）在自定义函数中，除了数组外还有其他值需要返回怎么办？
将这些值存入数组中一起返回

4. 二进制枚举子集法

   for(int i=0; i<(1<<n); i++) //从0～2^n-1个子集
    {
        for(int j=0; j<n; j++) //遍历二进制的每一位
        {
            if(i&(1<<j))//判断二进制第j位是否存在
            {
                printf("%d ",j);//如果存在输出第j个元素
            }
        }
        printf("\n");
    }

5. 遍历：

for(auto&i:j) 

等价于

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