cogs 2521. 首遇lancer

2521. 首遇lancer

★☆   输入文件： lancer.in   输出文件： lancer.out    简单对比
时间限制：1 s   内存限制：256 MB

【题目背景】

第五次圣杯之战。kito乱入。

【题目描述】

    一个宁静的夜晚，kito在学校的楼顶乘凉。

   “多么美好的星空啊，可惜只有在学校才能看到。”kito发表感慨。

    天上闪过一道星光，划过天际，kito正以为看到了少有的流星，谁知一声巨响降落在楼顶。

    “哟，年轻人。”烟雾散去，是lancer。

    “lancer？天哪，你找我干什么，虽然我乱入了，但是我还没有servant。”kito立马认怂。

    “没事，今天不是来打架的，”lancer帅气的一笑，kito不禁看呆。

    “saber还没出现，还没到我的戏份，所以闲的无聊，”lancer说，“于是就去破解教主的密码箱，但是密码太恶心了，像我这种帅气的英灵自然是不会的啦。听说你最近刚复习完数论，我来请教一下一个正整数的因数个数”。

    “枚举一下不就好了。”

    “但是它是求和啊，要求n的所有因数的的因数个数之和。”

    “……”

    “由于教主讨厌大于等于109+7的数，所以你只需要把答案对109+7取模后的值告诉我。如果你答对了，我就和你一起分享密码箱里的宝物，同时陪你一天哦～。”

   一听lancer会陪他一天，kito非常的激动，但是kito实在是太弱了，所以他只能拜托你，为了kito的幸福生活，希望你能把答案告诉他。

【输入格式】

   输入一共一行，有一个整数，为n。

【输出格式】

   共一行，为 ∑i=1 N [Nmod i==0]F[i]对109+7取模的值，F[i]表示i的因数个数。

【输入样例1】

   1

【输出样例1】

   1

【输入样例2】

10

【输出样例2】

9

【数据范围】

   30%的数据满足n<=5000

   70%的数据满足n<=100000

   90%的数据满足n<=10000000
   100%的数据满足n<=1000000000

【来源】

   Hzoi 2016

这种数学题， 一般一眼看不出答案， 那就手动模拟几个数， 再找一找规律， 一般答案就出来啦；

我们先举个例子： 把18分解：

eg: 分解 18;
primes    2, 3, 5, 7 ,,,,,
e[]           :1,  2, 0, 0  ,,,,,  


  因子          因子数 
   1  :0, 0;     1*1       （0， 0）指的是2^0 *3^0==1 下同
   3  :0, 1;     1*2
   9  :0, 2;     1*3
   2  :1, 0;     2*1
   6  :1, 1;     2*2
   18 :1, 2;    2*3

解释一下： 18的因子总共有（1， 3， 9， 2， 6， 18）这么几个；

1的因子数为1 （1*1）；

3的因子数为2 （1*2）；

9的因子数为3 （1*3）；

慢慢地就会发现答案  ans=(0+1, 1+1...e[2]+1)*(0+1, 1+1...e[3]+1);

再由个体推广到一般   ans=ans*((1+e[i]+1)/2 * (e[i]+1) (i=0, 1, 2。。。。);  （等差数列求和公式首项1，加上末项e[i]+1, 除以2， 再乘项数e[i]+1）

完毕；

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

#define LL long long 

const int maxn=31623;
const int INF=0x3f3f3f3f;

LL e[maxn], ans=1;
int n, cnt;

void fj(int n) {
	for(int i=2; i<=n; i++) {
		while(n%i==0) {
			n/=i;
			e[cnt]++;
		}
		cnt++;
	}
}

int main() {
	freopen("lancer.in","r",stdin);
    freopen("lancer.out","w",stdout);
	scanf("%d", &n);
	fj(n);
	for(int i=0; i<cnt; i++) ans=ans*(e[i]+2)*(e[i]+1)/2;       //因数据太水， 此处不取模也没事；
	cout<<ans;
	return 0;
}