关于自然数的高次幂的问题

欠你一滴泪

请教一个关于高次幂的问题:
我们要求计算1 + 2 + 3 + … + n, 而n很大（例如1000000）时，大概我们会用求和公式 S1(n) = 1 + 2 + 3 +… + n = n (n+1)/2 计算，而不会直接去累加。
对于 1^2 + 2^2 + 3^2 + … + n^2 , 我们有求和公式
S2 (n) = 1^2 + 2^2 + 3^2 + … + n^2 = n (n+1) (2n+1)/6
现在考虑较一般点的问题
Sp (n) = 1^p + 2^p+ 3^p + … + n^p （p是正整数）
输入 n，p；
输出结果；
???
编程思想???

爱琴海

那就用高精度



by 雨中飞燕 QQ:78803110 C/C++讨论群:46520219
Blog: http://yzfy.programfan.com

请大家不要用TC来学习C语言，点击此处查看原因 请不要写出非int声明的main函数
C++编写的Windows界面游戏

天之骄灵

输出的不要科学计数，要精确的值啊；
比如：
8 9999
111061117777777731111111333333333000

爱琴海

那你要输出什么？



by 雨中飞燕 QQ:78803110 C/C++讨论群:46520219
Blog: http://yzfy.programfan.com

请大家不要用TC来学习C语言，点击此处查看原因 请不要写出非int声明的main函数
C++编写的Windows界面游戏

▄愛變鎖ゞ

请教一下...怎么解决输出是科学计数的问题呢？

‖蠻忝謃☆

#include<iostream>
using namespace std;
double fun(int n,int p)
{
double i;
double m=1;
for(i=1;i<=p;i++)
m*=n;
return m;
}
int main()
{
double n,p,i;
long ans=0;
cin>>p>>n;
for(i=1;i<=n;i++)
ans+=fun(i,p);
cout<<ans<<endl;
return 0;
}
怎么解决输出是科学计数的问题呢？？？？

无民

1. 
   
2. void main()
   
3. {
   
4. double sum = 1, total = 0;
   
5. int n, p;
   
6. scanf("%d%d", &n, &p);
   
7. 
   
8. for (int i=1; i<=n; i++)
   
9. {
   
10. 
    
11. for (int j=1; j<=p; j++)
    
12. {
    
13. sum = sum * i;
    
14. 
    
15. }
    
16. total = total + sum;
    
17. sum = 1;
    
18. }
    
19. printf("sum = %.1f \n", total);
    
20. }
    
21. 
    

复制代码

肖伊紫轩

但是好像是不行的哦...

℡可可℡

谢谢，我想想..

泪在爱上雕刻

如果单纯从数学公式来说，
进一步推导是用二项公式做的
比如说吧：
1. （1+n)^p-(n)^p=a0+a1s(1)+a2s(2)+...+ap-1s(p-1)
2. (n)^p-(n-1)^p=a0+a1s(1)+a2s(2)+...+ap-1s(p-1)
3. (n-1)^p-(n-2)^p=a0+a1s(1)+a2s(2)+...+ap-1s(p-1)
..
n. (2)^p-1^p=a0+a1s(1)+a2s(2)+...+ap-1s(p-1)

再把所有的项相加就可以了，实际上就是利用前s(1),s(2),...,s(n-1)的公式来求s(n）
对于（1+n)^p-(n)^p=（c(p,0)+c(p,1)*n+c(p,2)*n^2,...+c(p,p)*n^p-n^p)
就可以求得：
（1+n)^p-1^p=c(p,0)s(0)+c(p,1)*s(1)+c(p,2)*s(2)+c(p,p-1)*n^p;
所以很容易就算得了三阶，四阶，五阶。。。的公式了
如果是编的话，可以用这个递推公式来做了，
但是，我觉得还是直接算来得方便，而且简单
不用递推公式：
#include <iostream.h>
long fnsub(n,p)
{ long fresult=n;
long i=1;
for(i=1;i<=p;i++)
fresult*=i;
return fresult;
}

void main()
{
int p;
long num;
cout<<"num:"<<endl;
cin>>num>>p;
for（int i=1;i<=num;i++)
fs=fs+fnsub(i,p) ;
cout<<fs<<endl;
}