cogs 2521. 首遇lancer

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_34004420/article/details/53046374

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=1N [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;
}



展开阅读全文

没有更多推荐了,返回首页