[已解决]救救孩子吧

llf123456

根据公式
arctanx(x)=x - x^3/3 + x^5/5 - x^7/7 +......(-1)^(n+1) x^(2n-1)/(2n-1)
和π=6arctanx(1/√3)，求当最后一项小于 10^(-t)时π的值。 t 由用户从键盘输入，t
等于 1 时，表示 arctanx(x)的数列的最后一项要小于 0.1；t 等于 2 时，表示最后一项要小
于 0.01。 输出的π需要保留 8 位小数。


#include <stdio.h>
#include <math.h>
int main()
{
        int i,j,e,t;
        double T=0.0,sum=0.0,fenzi,fenmu,x=1.0/sqrt(3);
        scanf("%d",&t);
        for(i=e=1;;i++,e=-e){
                for(j=1,fenzi=1;j<=(2*i-1);j++)fenzi*=x;
            //for(j=1,fenmu=-1;j<=(2*i-1);j++)fenmu+=2；//fenmu的运算为什么不能用这个循环体，而要用楼下的才行**************
            fenmu = 2*i-1;
            T = e*fenzi/fenmu;
            sum+=T;
            if(fenzi/fenmu<pow(0.1,t))break;
        }
        printf("%.8lf",6*sum);
}
输入样例 5
输出样例 3.14159977

最佳答案

月排行榜 / 总排行榜

jackz007

2020-12-27 16:34:51

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. double fun(int t)
   
5. {
   
6.         int e , i , j                                    ;
   
7.         double d , r , s , x                             ;
   
8.         for(r = 1 , i = 0 ; i < t ; i ++) r /= 10        ;
   
9.         for(x = 1.0 / sqrt(3) , s = 0 , e = i = 1 ; ; i ++ , e = - e) {
   
10.                 for(d = j = 1 ; j < 2 * i ; j ++) d *= x ;         // d = pow(x , 2 * i - 1)
    
11.                 d = d / (2 * i - 1)                      ;
    
12.                 s += e * d                               ;
    
13.                 if(d < r) break                          ;
    
14.         }
    
15.         return s * 6                                     ;
    
16. }
    
17. 
    
18. int main(void)
    
19. {
    
20.         float x                                          ;
    
21.         int t                                            ;
    
22.         scanf("%d" , & t)                                ;
    
23.         printf("%.8lf\n" , fun(t))                       ;
    
24. }

复制代码

        编译、运行实况

1. D:\00.Excise\C>g++ -o x x.c
   
2. 
   
3. D:\00.Excise\C>x
   
4. 5
   
5. 3.14159977
   
6. 
   
7. D:\00.Excise\C>

复制代码

跳转到最佳答案楼层

jackz007

        楼主，有测试样例吗？

llf123456

jackz007 发表于 2020-12-27 15:33
楼主，有测试样例吗？

输入样例 5
输出样例 3.14159977

jackz007

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. double fun(int t)
   
5. {
   
6.         int e , i , j                                    ;
   
7.         double d , r , s , x                             ;
   
8.         for(r = 1 , i = 0 ; i < t ; i ++) r /= 10        ;
   
9.         for(x = 1.0 / sqrt(3) , s = 0 , e = i = 1 ; ; i ++ , e = - e) {
   
10.                 for(d = j = 1 ; j < 2 * i ; j ++) d *= x ;         // d = pow(x , 2 * i - 1)
    
11.                 d = d / (2 * i - 1)                      ;
    
12.                 s += e * d                               ;
    
13.                 if(d < r) break                          ;
    
14.         }
    
15.         return s * 6                                     ;
    
16. }
    
17. 
    
18. int main(void)
    
19. {
    
20.         float x                                          ;
    
21.         int t                                            ;
    
22.         scanf("%d" , & t)                                ;
    
23.         printf("%.8lf\n" , fun(t))                       ;
    
24. }

复制代码

        编译、运行实况

1. D:\00.Excise\C>g++ -o x x.c
   
2. 
   
3. D:\00.Excise\C>x
   
4. 5
   
5. 3.14159977
   
6. 
   
7. D:\00.Excise\C>

复制代码

pa1st

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. int main()
   
5. {
   
6.     double x;
   
7.     int n;
   
8.     scanf("%lf %d",&x,&n);
   
9.     int i;
   
10.     double sin,sum=0;
    
11.     for(i=1;i<=n;i++)
    
12.     {
    
13.             sin=pow(-1,i-1)*pow(x,2*i-1)/fact(2*i-1);
    
14.             sum=sum+sin;
    
15.         }
    
16.         printf("%.8lf",sum);
    
17.         return 0;
    
18. }
    
19. 
    
20. int fact(int i)
    
21. {
    
22. 
    
23.          int pro=1;
    
24. 
    
25.          int m=1;
    
26. 
    
27.          for(;m<=i;++m)
    
28. 
    
29.                    pro=pro*m;
    
30. 
    
31.          return(pro);
    
32. }

复制代码