关于约数问题

几天前做了一个求约数方面的题，因为数据不大采用了打表法。个人觉得这种方案，有点无耻。于是上网查找了一下求约数的方法，现在总结一下，分成如下几类：

1.求出所有的质约数；

2.求所有约数的和；

3.求出所有约数的个数；

4.求出所有质约数指数的和；

求约数个人总结了一个还算比较高效的方案模板如下：

/*ps：这段程序默认不求1和负数，思路：连续分解成质因数，
使用贪心的法则，将质因数和对应的指数存入map中*/
map <int,int> m;
map <int,int>::iterator it;
void submultiple(int n)
{
	int i;
	m.clear();
	for(i=2 ; i*i<=n ; i++)
	{
		while(n%i==0) m[i]++,n/=i;
	}
	if(n!=1) m[n]++;
}

这个模板里面包含了，所有的质约数，也就是it->first；而质约数指数则是it->second；解决1，4两份！

对于所有的约数和和所有的约数的个数其实是有公式的：

首先我们定义一下约数——约数指的是a%b==0即b是a的约数！这里我们先假设质因数的指数用bi表示，质因数本身用ai表示；

约数的个数=(a1+1)*(a2+1)*(a3+1)*......*(an+1)；

约数的和=(1+a1^1+a1^2+......a1^b1)*(1+a2^1+a2^2+.....a2^b2)*........*(1+an^1+an^2+.....+an^bn)；

如果知道这两个公式，然后套用前面的程序所求出来的值，估计都是小菜一碟了吧。