fork download
copy 
  1. #include <stdio.h>
  2. #include <math.h>
  3.  
  4. typedef unsigned int Uint;
  5. #define N 300
  6. Uint DP=1;
  7. //RADIX == 10^LOG
  8. #define LOG 1
  9. #define RADIX 10
  10. typedef Uint Liczba[N+1];
  11.  
  12. void pisz( Liczba a, int ctrl){ Uint i;
  13. 	for( i=1; i<DP && !a[i]; i++ ) printf("%*s", LOG,"");
  14. 	printf("%*u", LOG, a[i]);
  15. 	for( i++; i<=DP; i++ ) printf("%0*u", LOG, a[i]);
  16. 	putchar('.');
  17. 	for(    ; i<=N; i++ ) printf("%0*u", LOG, a[i]); 
  18. 	if( ctrl < 0 ) puts("");}
  19.  
  20. void plus( Liczba a, Liczba b){                  // a[]:=a[]+b[]
  21. 	for( Uint c=0, i=N; i>0; i-- ){
  22. 		c += a[i]+b[i];
  23. 		a[i] = c % RADIX;
  24. 		c /= RADIX;}}
  25.  
  26. void minus( Liczba a, Liczba b ) {               // a[]:=a[]-b[]
  27. 	for( Uint c=0, i=N; i>0; i-- ){
  28. 		c += b[i];
  29. 		if( a[i]<c ){
  30. 			a[i] = RADIX + a[i] - c; c=1; }
  31. 		else {
  32. 			a[i] -= c; c=0; }}}
  33.  
  34. void razyU( Liczba a, Uint b ) {                 // a[]:=a[] * b
  35. 	for( Uint c=0,i=N; i>0; i-- ){
  36. 		c += a[i]*b;
  37. 		a[i]= c % RADIX;
  38. 		c /= RADIX; }}
  39.  
  40. int podzielU( Liczba w, Liczba a, Uint b ) {     // w[]:=a[] / b
  41. 	Uint z=0;                                // return w[] != 0
  42. 	for( Uint c=0, i=1; i<=N; i++ ) {
  43. 		c = c*RADIX+a[i];
  44. 		w[i]= c/b % RADIX; if( w[i] ) z=1;
  45. 		c = c - w[i]*b;}
  46. 	return z;}
  47.  
  48. void wstawU( Liczba a, Uint b ) {                // a[]:= b
  49. 	Uint i;
  50. 	for( i=0; i<=N; i++ ) a[i]=0;
  51. 	i=DP;
  52. 	while( b ) {
  53. 		a[i--] = b%RADIX; b /= RADIX; }}
  54.  
  55. void atg( Liczba a, Uint x ) {                  // a[]:= atan( 1/x )
  56. 	Liczba px;
  57. 	wstawU( px, 1 ); podzielU( px, px, x );
  58. 	wstawU( a,  0 ); plus    ( a, px );
  59. 	for( Uint n=1;   ; n++) {
  60. 		Liczba b;
  61. 		podzielU( px, px, x*x );
  62. 		if( ! podzielU( b, px, 2*n + 1 ) ) 
  63. 			break;
  64. 		if( n & 1 ) minus( a, b );
  65. 		else plus( a, b );}}
  66.  
  67. void silnia( Liczba a, Uint n ){                 // a[]:= n!
  68. 	wstawU( a, 1 );
  69. 	while( n>1 )
  70. 		razyU( a, n-- );}
  71.  
  72. int bbp( Liczba pi ){
  73. 	wstawU( pi, 0 );
  74. 	Liczba a, b;
  75. 	for( Uint k=0;  ; k++ ) {
  76. 		wstawU( a, 4 ); podzielU( a, a, 8*k + 1 );
  77. 		wstawU( b, 2 ); podzielU( b, b, 8*k + 4 ); minus( a, b );
  78. 		wstawU( b, 1 ); podzielU( b, b, 8*k + 5 ); minus( a, b );
  79. 		wstawU( b, 1 ); podzielU( b, b, 8*k + 6 ); minus( a, b );
  80. 		for( int i=1; i<=k; i++ )
  81. 			if( ! podzielU( a, a, 16 ) )
  82. 				return k;
  83. 		plus( pi, a ); }}
  84.  
  85. int byat( Liczba pi ) {
  86. 	wstawU( pi, 1 );
  87. 	Uint n=1, d=3;
  88. 	Liczba a; wstawU( a, 1 );
  89. 	while( 1 ) {
  90. 		razyU( a, n ); n++;
  91. 		if( ! podzielU( a, a, d ) ) {
  92. 			razyU( pi, 2 );
  93. 			return n;
  94. 		d += 2;
  95. 		plus( pi, a );}} 
  96.  
  97. int main(void) {
  98. 	Liczba a;
  99. 	DP=1;
  100. 	if( 0 ) {
  101. 		Liczba b;
  102. 		atg  ( a, 18 ); 
  103. 		razyU( a, 12 );
  104. 		atg  ( b, 57 );
  105. 		razyU( b,  8 );
  106. 		plus ( a,  b );
  107. 		atg  ( b,239 );
  108. 		razyU( b,  5 );
  109. 		minus( a,  b );
  110. 		razyU( a,  4 );
  111. 		pisz ( a, -1 );
  112. 	} else {
  113. 		printf("%d\n", byat  ( a ) );
  114. 		pisz ( a, -1 );
  115. 	}
  116. 	DP=N;
  117. 	silnia(a, 57 );
  118. 	pisz ( a, -1 );
  119. 	return 0;
  120. }
  121.  
Compilation error #stdin compilation error #stdout 0s 0KB
comments (?)
stdin
copy 
Standard input is empty
compilation info 
cc1: warnings being treated as errors
prog.c: In function ‘byat’:
prog.c:97: error: ISO C forbids nested functions
prog.c:97: error: ‘main’ is normally a non-static function
prog.c:120: error: expected declaration or statement at end of input
stdout
copy 
Standard output is empty
https://ideone.com/Fr4A5
language:
C99 (gcc 8.3)
created:
14 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!