首页 > 代码库 > HDU 4961 Boring Sum 暴力
HDU 4961 Boring Sum 暴力
题意:对于所有的A[I],同时找到左边和右边离它最近且是它的倍数的数相乘最后加起来求和。
解题思路:n*sqrt(n)的算法,开始以为过不了,wa了两发因为lld I64d对拍一个小时发现一个小时前交的代码没错只是没变I64d,。。具体思路是枚举每个a[i]的因子,找离它最近的那个更新,如果已经没更新就不用更新了。用两个辅助数组记录最近的就行。
解题代码:
1 // File Name: 1002.cpp 2 // Author: darkdream 3 // Created Time: 2014年08月19日 星期二 14时10分33秒 4 5 #include<vector> 6 #include<list> 7 #include<map> 8 #include<set> 9 #include<deque>10 #include<stack>11 #include<bitset>12 #include<algorithm>13 #include<functional>14 #include<numeric>15 #include<utility>16 #include<sstream>17 #include<iostream>18 #include<iomanip>19 #include<cstdio>20 #include<cmath>21 #include<cstdlib>22 #include<cstring>23 #include<ctime>24 #define LL long long25 26 using namespace std;27 int hs[100005];28 int hs1[100006];29 int ans[100005][2];30 int a[100005];31 int main(){32 //freopen("input.txt","r",stdin);33 //freopen("output1","w",stdout);34 int n ; 35 while(scanf("%d",&n) != EOF,n)36 {37 38 memset(hs,0,sizeof(hs));39 memset(ans,0,sizeof(ans));40 memset(hs1,0,sizeof(hs1));41 for(int i = 1;i <= n;i ++)42 {43 scanf("%d",&a[i]);44 int p = sqrt(a[i]+1e-8);45 for(int j = 1;j <= p ;j++)46 {47 if(a[i] % j == 0 )48 {49 if(!ans[hs[j]][1])50 ans[hs[j]][1] = a[i];51 if(!ans[hs[a[i]/j]][1])52 ans[hs[a[i]/j]][1] = a[i];53 }54 }55 hs[a[i]] = i;56 }57 for(int i = n;i >= 1;i --)58 {59 int p = sqrt(a[i]+1e-8);60 for(int j = 1;j <= p ;j ++)61 {62 if(a[i]%j == 0 )63 {64 if(!ans[hs1[j]][0])65 ans[hs1[j]][0] = a[i];66 if(!ans[hs1[a[i]/j]][0])67 ans[hs1[a[i]/j]][0] = a[i];68 }69 }70 hs1[a[i]] = i ; 71 }72 LL sum = 0 ; 73 for(int i = 1;i <= n;i ++)74 {75 for(int j = 0 ;j <= 1;j ++)76 if(!ans[i][j])77 ans[i][j] = a[i];78 //printf("%d %d\n",ans[i][0],ans[i][1]);79 sum += 1ll*ans[i][0] * ans[i][1];80 }81 printf("%I64d\n",sum);82 }83 return 0;84 }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。