ZOJ 2095 数论,因子之和

本文介绍了一种高效的因子分解算法,用于解决特定数学问题。通过观察数的因子规律,并结合牛顿二项式原理,实现了快速计算因子的方法。文章还提供了一个使用素数筛选法分解因式的C语言实现示例。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目连接

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1095

刚刚看到这个题的数据规模的时候本来想用暴力,因为题目的数据量只有50w,但是很悲剧的TLE了

于是开始认认真真的想算法

可以很容易观察到,一个如下的规律

比如:20可以拆分为:2^2*5^1

那么就是1,2,4和1,5的组合相乘就是每一个因子

联想到牛顿二项式

得到:

(1+2+4)*(1+5)可以得到该答案为42,由于这个算式的和包括了一个特殊的因子,即该数n本身,于是我们必须将它减去就得到42-20=22

那么我们就可以很容易的写出程序来了(利用素数筛选法进行分解因式加速可以在500MS左右解决问题)

我的代码:

#include<stdio.h> #include<string.h> int prime[1000]; bool flag[1000]; void init() { int i,j; int num=0; memset(flag,0,sizeof(flag)); for(i=2;i<1000;i++) { if(!flag[i]) { prime[num++]=i; for(j=i*i;j<1000;j=j+i) flag[j]=true; } } } void slove(int n) { int i,time,ans=1,N=n,k1,k2; for(i=0;prime[i]*prime[i]<=n;i++) { time=0; k1=1; k2=0; if(n%prime[i]==0) { k1=k1*prime[i]; k2=k2+k1; n=n/prime[i]; while(n%prime[i]==0) { n=n/prime[i]; k1=k1*prime[i]; k2=k2+k1; } ans=ans*(k2+1); } if(n==1) break; } k1=1; k2=0; if(n>1) { k1=n; k2=k2+k1; ans=ans*(k2+1); } ans=ans-N; printf("%d/n",ans); } int main() { int n,t,ans; init(); scanf("%d",&t); while(t--) { ans=0; scanf("%d",&n); slove(n); } return 0; }

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值