hdu 4961 Boring Sum

最新推荐文章于 2022-02-07 21:25:06 发布

云起2013

最新推荐文章于 2022-02-07 21:25:06 发布

阅读量458

点赞数

CC 4.0 BY-SA版权

分类专栏：朴素的劳动者基础文章标签： hdu 简单题数组多校联合 SYSU

本文链接：https://blog.youkuaiyun.com/u012355629/article/details/38707683

朴素的劳动者同时被 2 个专栏收录

2 篇文章

订阅专栏

基础

2 篇文章

订阅专栏

本文详细解析了BoringSum问题的解题思路，包括遍历法的运用及优化，通过代码实例展示了如何解决该问题，并讨论了代码背后的逻辑与思考过程。

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

原题： hdu 4961 Boring Sum

来源： 2014 Multi-University Training Contest 9

参考：林燕同学的博客

做题感受：比赛时坚持了原始朴素的遍历法，不出意外的TLE，最后是小伙伴把它A掉了，但是他的代码我看不懂（忧桑）。赛后看了题解，

感觉智商不够用，最后还是百度了先A者的做题代码，敲完了，琢磨了许久，脑子才转过弯来。

参考代码如下：

#include<iostream>
#include<string>
#include<cstdio>
using namespace std;
#define SIZE 100005
#define m(a,b) memset(a,b,sizeof(a))
int a[SIZE],b[SIZE],c[SIZE],p[SIZE];
int main()
{
	int n,i,j;
	__int64 sum;
	while(scanf("%d",&n),n!=0)
	{
		sum=0;
		m(p,-1);
		for(i=0;i<n;i++)       //自左向右扫描
		{
			scanf("%d",&a[i]);
			if(p[a[i]]>-1)     //存在能被a[i]整除的数，该数必是离当前数最近且为a[i]的倍数
				b[i]=p[a[i]];
			else
				b[i]=i;        //a[i]左侧不存在能被a[i]整除的数
			for(j=1;j*j<=a[i];j++)
			{
				if(a[i]%j==0)  //更新a[i]各因子所对应的最大的数列元素下标，
					           //保证下次查询元素倍数时找到的元素是离查询元素最近的倍数
				{
					p[j]=i;
					p[a[i]/j]=i;
				}
			}
		}
		m(p,-1);
		for(i=n-1;i>=0;i--)    //自右向左扫描
		{
			if(p[a[i]]>-1)
				c[i]=p[a[i]];
			else
				c[i]=i;        //a[i]右侧不存在能被a[i]整除的数
			for(j=1;j*j<=a[i];j++)
			{
				if(a[i]%j==0)
				{
					p[j]=i;
					p[a[i]/j]=i;
				}
			}
		}
		for(i=0;i<n;i++)
			sum=sum+(__int64)a[b[i]]*a[c[i]];
		printf("%I64d\n",sum);
	}
	return 0;
}