HDU-2139 Calculate the formula

本文详细介绍了如何解决HDU-2139问题,通过数学归纳法推导出计算奇数平方和的通式,并提供了一段简洁高效的C++代码实现。

HDU-2139 Calculate the formula

题目描述

You just need to calculate the sum of the formula: 1^2+3^2+5^2+……+ n ^2.

Input

In each case, there is an odd positive integer n.

Output

Print the sum. Make sure the sum will not exceed 2^31-1

Sample Input
3

Sample Output
10

思路

我绝对不会说最后就算找到公式,还是TLE;
记住,这一道题目的输入很大,不要用cin & cout!!!;

一开始就想到模拟做肯定白搭,所以用到了数学归纳法求解通式;

先给出求前n项平方和的证明(经典解法):

    设S=1^2+2^2+.+n^2
    已知(n+1)^3-n^3 = 3n^2+3n+1
    n^3-(n-1)^3 = 3(n-1)^2+3(n-1)+1
    ...
    ...
    ...
    2^3-1^3 = 3*1^2+3*1+1
    把上面n个式子相加得:(n+1)^3-1 = 3* [1^2+2^2+...+n^2] +3*[1+2+.+n] +n
    所以S= (1/3)*[(n+1)^3-1-n-(1/2)*n*(n+1)] = (1/6)*n*(n+1)*(2n+1)
    即 1^2+2^2+3^2+...+n^2=n*(n+1)*(2n+1)/6;

    我们把每一项都乘上2,那么上式就变为:2^2+4^2+6^2+...+(2n)^2=4*(1^2+2^2+3^2+...+n^2)=4*n*(n+1)*(2n+1)/6;

    即前2n项中偶数的平方和,

    由第一个式子,我们又可以知道前(2n+1)项数的平方和为:1^2+2^2+3^2+...+n^2+(n+1)^2+...+(2n+1)^2=(2n+1)*(2n+1+1)*(4n+2+1)/6;

    那么用第三个式子减去第二个式子,便得到了前2n+1个数中奇数的平方和,即前n个奇数的平方和:

    1^2+3^2+5^2+...+(2n+1)^2=(n+1)*(2n+1)*(2n+3)/3,

    所以就得到了前n个数中奇数的平方和公式:n*(n+1)*(n+2)/6.

代码

#include <iostream>
#include <algorithm>

using namespace std;
typedef unsigned long long ULL;

int main()
{
    long long n;
    while(~scanf("%lld",&n))
    {
        printf("%lld\n",n*(n+1)*(n+2) / 6);
    }
    return 0;
}
### 关于HDU - 6609 的题目解析 由于当前未提供具体关于 HDU - 6609 题目的详细描述,以下是基于一般算法竞赛题型可能涉及的内容进行推测和解答。 #### 可能的题目背景 假设该题目属于动态规划类问题(类似于多重背包问题),其核心在于优化资源分配或路径选择。此类问题通常会给出一组物品及其属性(如重量、价值等)以及约束条件(如容量限制)。目标是最优地选取某些物品使得满足特定的目标函数[^2]。 #### 动态转移方程设计 如果此题确实是一个变种的背包问题,则可以采用如下状态定义方法: 设 `dp[i][j]` 表示前 i 种物品,在某种条件下达到 j 值时的最大收益或者最小代价。对于每一种新加入考虑范围内的物体 k ,更新规则可能是这样的形式: ```python for i in range(n): for s in range(V, w[k]-1, -1): dp[s] = max(dp[s], dp[s-w[k]] + v[k]) ``` 这里需要注意边界情况处理以及初始化设置合理值来保证计算准确性。 另外还有一种可能性就是它涉及到组合数学方面知识或者是图论最短路等相关知识点。如果是后者的话那么就需要构建相应的邻接表表示图形结构并通过Dijkstra/Bellman-Ford/Floyd-Warshall等经典算法求解两点间距离等问题了[^4]。 最后按照输出格式要求打印结果字符串"Case #X: Y"[^3]。 #### 示例代码片段 下面展示了一个简单的伪代码框架用于解决上述提到类型的DP问题: ```python def solve(): t=int(input()) res=[] cas=1 while(t>0): n,k=list(map(int,input().split())) # Initialize your data structures here ans=find_min_unhappiness() # Implement function find_min_unhappiness() res.append(f'Case #{cas}: {round(ans)}') cas+=1 t-=1 print("\n".join(res)) solve() ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值