高精度浮点数算法 计算1-1/2+1/3-1/4+··· ···1/999-1/1000

s=1-1/2+1/3-1/4+……+1/999-1/1000，要求精确到小数点后50位。

float的有效精度是8位，double是16位，都不满足我们的要求，所以只能用数组来处理这个问题。加减每次从最低为开始操作，每次计算处理进位，除法从最高位计算，每次计算处理进位

PS:这个题我写的所有代码只适用于同类型的题，比如你可以让精度到100位或者加到1/10000。因为我的除法，加法，减法都只考虑了正数的情况，显示中可能有负数的情况，同时减法考虑的是被减数必定大于减数（这是这个题里必定具有的性质）。

基本没用啥c++的特性，如果要用纯c把输入输出和头文件改下就好了

//高精度算法
#include<iostream>
using namespace std;

//精度，修改精度直接修改size就好 
const int SIZE=51; 
int num[SIZE]={1};

//除法
int* div(int *arr,int n)
{
	int t=0,s;
	for(int i=0;i<SIZE;i++)
	{
		//s是上面一个数取模得来的值（t）+本来位置的值
		s=arr[i]+t*10;
		arr[i]=(s/n)%10;
		//除不尽的要留给下一位
		t=s%n;
	}
	return arr;
}

//加减法一起写了，key是1加法，0减法
int* add(int* a,int* b,int key)
{
	int t=0;
	if(key)
	{

		for(int i=SIZE-1;i>=0;--i)
		{

			a[i]+=(b[i]+t);
			t=a[i]/10;
			a[i]%=10;
		}

	}else
	{
		for(int i=SIZE-1;i>=0;--i)
		{
			//去掉后一个数借的位，同时把t重置回0，因为不一定会借1
			a[i]-=(b[i]+t);
			t=0;
			//如果当前位置不够减，向前借1
			if(a[i]<0)
			{
				a[i]+=10;
				t=1;
			}
		}
	}
	return a;
	
}

//结果
int* get_result(int *arr,int n)
{
	for(int i=2;i<=n;i++)
	{
		//从第二项开始处理，因为第一个是1无需处理
		arr=add(arr,div(num,i),i%2);
		num[0]=1;
		for(int i=1;i<SIZE;i++)
		{
			num[i]=0;
		}
	}
	return arr;
}

int main()
{
	int arr[SIZE]={1};
	//计算分母到1000的值 
	int*ans=get_result(arr,1000);
	for(int i=0;i<SIZE-1;++i)
	{
		if(i==0)
		{
			cout<<arr[i]<<'.';
		}else		
		cout<<ans[i];
	}
	return 0;
 }