C++编程-高精度算法（STL表标程）

#include <cstdlib>
#include <cstdio>
#include <algorithm>
int array[103];
int out[103];
using namespace std;
int main()
{
    int n;
    scanf("%d",&n);
    int n_n=n;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&array[i]);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(array[i]==array[j]&&i!=j)
            {
                n_n--;
                array[i]=-1;
            }
        }
    }
    printf("%d\n",n_n);
    int c=1;
    for(int i=1;i<=n;i++)
    {
        if(array[i]!=-1)
        {
            out[c]=array[i];
            c++;
        }
    }
    sort(out+1,out+n_n+1);
    for(int i=1;i<=n_n;i++)
    {
        printf("%d ",out[i]);
    }
    //system("pause");
    return 0;
}

二：剔除相关的数

本题暂无标准程序

高精度算法介绍

先言

利用计算机进行数值计算，有时会遇到这样的问题：有些计算要求精度高，希望计算的数的位数可达几十位甚至几百位，虽然计算机的计算精度也算较高了，但因受到硬件的限制，往往达不到实际问题所要求的精度，且long long或者无字符long long也无法满足需求。我们可以利用程序设计的方法去实现这样的高精度计算。介绍常用的几种高精度计算的方法。

高精度计算中需要处理好以下几个问题：

(1)数据的接收方法和存贮方法

数据的接收和存贮：当输入的数很长时，可采用字符串方式输入，这样可输入数字很长的数，利用字符串函数和操作运算，将每一位数取出，存入数组中。另一种方法是直接用循环加数组方法输入数据。

	void init(int a[])                            //传入一个数组
	{
            string s; 
	       cin>>s;                                  //读入字符串s 
	       len=s.length();                     //用len计算字符串s的位数 
	       for(i=1;i<=len;i++)
	          a[i]=s[len-i] -‘0’;               //将数串s转换为数组a，并倒序存储 
	}

另一种方法是直接用循环加数组方法输入数据。

(2)高精度数位数的确定

位数的确定：接收时往往是用字符串的，所以它的位数就等于字符串的长度。

(3)进位,借位处理

加法进位：

c[i]=a[i]+b[i];
if (c[i]>=10) { c[i]%=10; ++c[i+1]; }

减法借位：

if (a[i]<b[i]) { --a[i+1]; a[i]+=10; }
c[i]=a[i]-b[i];

乘法进位：

c[i+j-1]= a[i]*b[j] + x + c[i+j-1];
x = c[i+j-1]/10;
c[i+j-1] %= 10;

(4)商和余数的求法

商和余数处理：视被除数和除数的位数情况进行处理.

例题一：高精度加法

本题选自【信息学奥赛一本通高精度算法】

题目描述

高精度加法。输入两个正整数，求它们的和。

分析题目

输入两个数到两个变量中，然后用赋值语句求它们的和，输出。但是，我们知道，在C++语言中任何数据类型都有一定的表示范围。而当两个被加数很大时，上述算法显然不能求出精确解，因此我们需要寻求另外一种方法。在读小学时，我们做加法都采用竖式方法，如图1。这样，我们方便写出两个整数相加的算法。

如果我们用数组A、B分别存储加数和被加数，用数组C存储结果。则上例有A[1]=6，A[2]=5， A[3]=8，B[1]=5，B[2]=5，B[3]=2，C[4]=1，C[3]=1，C[2]=1，C[1]=1，两数相加如图2所示。

根据上方分析，我们可以写出代码的框架来了，框架如下

int c[100];
void add(int a[],int b[])                 //a,b,c都为数组，分别存储被加数、加数、结果
{
    int  i=1,x=0;                             //x是进位
    while ((i<=a数组长度)||(i<=b数组的长度))
　{
　　　　c[i]=a[i]+b[i]+x;    	//第i位相加并加上次的进位
　　　　x=c[i]/10;	             //向高位进位
　　　　c[i]%=10;                      //存储第i位的值
　　　　i++;                               //位置下标变量
　}
}

在通常的情况下，我们会将读入的两个整数放到两个char数组中，因此，标准程序如下

标准程序

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
    char a1[101],b1[101];
    int a[101],b[101],c[10001],lena,lenb,lenc,i,j,x;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    gets(a1);
    gets(b1);               //输入加数与被加数
    lena=strlen(a1);
    lenb=strlen(b1);
    for (i=0;i<=lena-1;i++) a[lena-i]=a1[i]-48;//加数放入a数组
　  for (i=0;i<=lenb-1;i++) b[lenb-i]=b1[i]-48;//加数放入b数组
    lenc =1;
    x=0;
	while (lenc <=lena||lenc <=lenb)
	{
	　　	c[lenc]=a[lenc]+b[lenc]+x;     //两数相加
	　　	x=c[lenc]/10;
	　　	c[lenc]%=10;
	       	lenc++;
	}
	c[lenc]=x; 
	if (c[lenc]==0)
	    lenc--;                                        //处理最高进位
	for (i=lenc;i>=1;i--) 
	   cout<<c[i];                                   //输出结果
	cout<<endl;
	return 0;
}

例题二：高精度减法

题目描述

高精度减法。输入两个正整数，求它们的差。

分析题目

高精减法类似加法，可以用竖式求减法。在做减法运算时，需要注意的是：被减数必须比减数大，同时需要处理借位。高精度减法的标准程序如下

标准程序

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
	int a[256],b[256],c[256],lena,lenb,lenc,i;
	char n[256],n1[256],n2[256];
	memset(a,0,sizeof(a));
	memset(b,0,sizeof(b));
	memset(c,0,sizeof(c));
	printf("Input minuend:");    gets(n1);   //输入被减数
	printf("Input subtrahend:"); gets(n2);   //输入减数
 	if (strlen(n1)<strlen(n2)||(strlen(n1)==strlen(n2)&&strcmp(n1,n2)<0))
                                         //strcmp()为字符串比较函数，当n1==n2, 返回0;
　　　　　　　　　　　 //n1>n2时，返回正整数；n1<n2时，返回负整数
	{                                        //处理被减数和减数，交换被减数和减数
　　   		strcpy(n,n1);        //将n1数组的值完全赋值给n数组
　　   		strcpy(n1,n2);
　　   		strcpy(n2,n);
　　   		cout<<"-";            //交换了减数和被减数，结果为负数
　　	}   
　　
　　  	lena=strlen(n1); lenb=strlen(n2);
　　  	for (i=0;i<=lena-1;i++) a[lena-i]=int(n1[i]-'0');  //被减数放入a数组
　　  	for (i=0;i<=lenb-1;i++) b[lenb-i]=int(n2[i]-'0');  //减数放入b数组
　　 	i=1;
	while (i<=lena||i<=lenb)
	{
		if (a[i]<b[i])
		{
			a[i]+=10;               //不够减，那么向高位借1当10
			a[i+1]--;
		}
		c[i]=a[i]-b[i];                        //对应位相减
		i++;
	}
	lenc=i;
	while ((c[lenc]==0)&&(lenc>1)) lenc--;   //最高位的0不输出　　  
	for (i=lenc;i>=1;i--) cout<<c[i];                //输出结果
	cout<<endl;
	return 0;
}

例题三：高精度乘法

题目描述

高精度乘法。输入两个正整数，求它们的积。

分析题目

高精度乘法类似加法，但有点难度，可以用竖式求乘法。在做乘法运算时，同样也有进位，同时对每一位进行乘法运算时，必须进行错位相加，如图1-3、图1-4。

分析c数组下标的变化规律，可以写出如下关系式：ci = c’i +c”i +…由此可见，c i跟a[i]*b[j]乘积有关，跟上次的进位有关，还跟原c i的值有关，分析下标规律，有c[i+j-1]= a[i]*b[j]+ x + c[i+j-1]; x=c[i+j-1]/10 ; c[i+j-1]%=10;

本题的标准程序如下

标准程序

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int main()
{
    char a1[101],b1[101];
    int a[101],b[101],c[10001],lena,lenb,lenc,i,j,x;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    scanf("%s",a1);                       //gets改成scanf
    scanf("%s",b1); 
    lena=strlen(a1);lenb=strlen(b1);
    for (i=0;i<=lena-1;i++) a[lena-i]=a1[i]-48;
    for (i=0;i<=lenb-1;i++) b[lenb-i]=b1[i]-48;
	for (i=1;i<=lena;i++)
	{
	x=0;                                              //用于存放进位
	     for (j=1;j<=lenb;j++)                    //对乘数的每一位进行处理
	     {
		     c[i+j-1]=a[i]*b[j]+x+c[i+j-1];  //当前乘积+上次乘积进位+原数	    	   
             x=c[i+j-1]/10;
		     c[i+j-1] %= 10;
	     }
	     c[i+lenb]=x;                                 //进位
	}
	lenc=lena+lenb;
	while (c[lenc]==0&&lenc>1)            //删除前导0
		lenc--;
	for (i=lenc;i>=1;i--)
		cout<<c[i];
	cout<<endl;
	return 0;
}

例题四：高精度除法（整除）

题目描述

高精除以低精。输入两个正整数，求它们的商（做整除）。

分析题目

做除法时，每一次上商的值都在０～９，每次求得的余数连接以后的若干位得到新的被除数，继续做除法。因此，在做高精度除法时，要涉及到乘法运算和减法运算，还有移位处理。当然，为了程序简洁，可以避免高精度除法，用0～9次循环减法取代得到商的值。这里，我们讨论一下高精度数除以单精度数的结果，采取的方法是按位相除法。

标准程序

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int main()
{
	char a1[100];
  	int a[100],c[100],lena,i,x=0,lenc,b;
  	memset(a,0,sizeof(a));
  	memset(c,0,sizeof(c));
  	gets(a1);
  	cin>>b;
  	lena=strlen(a1);
  	for (i=0;i<=lena-1;i++)
　　	    a[i+1]=a1[i]-48;
 	for (i=1;i<=lena;i++)                        //按位相除
	{
		c[i]=(x*10+a[i])/b;
	    x=(x*10+a[i])%b;
	}
　 	lenc=1;
    while (c[lenc]==0&&lenc<lena) lenc++;        //删除前导0
    for (i=lenc;i<=lena;i++) 
        cout<<c[i];
    cout<<endl;
    return 0;
}

例题五：高精度除法（商和余数）

题目描述

高精除以高精，求它们的商和余数。

分析题目

高精除以低精是对被除数的每一位（这里的“一位”包含前面的余数，以下都是如此）都除以除数，而高精除以高精则是用减法模拟除法，对被除数的每一位都减去除数，一直减到当前位置的数字（包含前面的余数）小于除数（由于每一位的数字小于10，所以对于每一位最多进行10次计算）具体实现程序如下：

标准程序

#include<iostream>
#include<cstring>
using namespace std;
int a[101],b[101],c[101],d,i;   
void init(int a[]) 
{    string s; 
	cin>>s;                        //读入字符串s 
	a[0]=s.length();           //用a[0]计算字符串 s的位数 
	for(i=1;i<=a[0];i++)
	a[i]=s[a[0]-i]-'0';           //将数串s转换为数组a，并倒序存储． 
}
void print(int a[])              //打印输出
{
	if (a[0]==0){cout<<0<<endl;return;}
	for(int i=a[0];i>0;i--) cout<<a[i];
	cout<<endl;
	return ;
}
int compare (int a[],int b[])  
			//比较a和b的大小关系，若a>b则为1，a<b则为-1,a=b则为0 
{    int i;
	if(a[0]>b[0]) return 1;                    //a的位数大于b则a比b大 
	if(a[0]<b[0]) return -1;                   //a的位数小于b则a比b小 
	for(i=a[0];i>0;i--)                           //从高位到低位比较 
	{
		if (a[i]>b[i]) return 1; 
		if (a[i]<b[i]) return -1;
	} 
	return 0;                                        //各位都相等则两数相等。 
} 
void numcpy(int p[],int q[],int det)      //复制p数组到q数组从det开始的地方
{
	for (int i=1;i<=p[0];i++) q[i+det-1]=p[i];
	q[0]=p[0]+det-1;
}
void jian(int a[],int b[])               //计算a=a-b
{ 
	int flag,i; 
	flag=compare(a,b);              //调用比较函数判断大小 
	if (flag==0) {a[0]=0;return;}   //相等 
	if(flag==1)                             //大于   
	{
		for(i=1;i<=a[0];i++) 
		{
			if(a[i]<b[i]){ a[i+1]--;a[i]+=10;}    //若不够减则向上借一位 
			a[i]-=b[i];
		} 
		while(a[0]>0&&a[a[0]]==0) a[0]--;           //修正a的位数 
		return;
	} 
}
void chugao(int a[],int b[],int c[])
{
	int tmp[101]; 
	c[0]=a[0]-b[0]+1;
	for (int i=c[0];i>0;i--)
	{
		memset(tmp,0,sizeof(tmp));                              //数组清零 
		numcpy(b,tmp,i);
		while(compare(a,tmp)>=0){c[i]++;jian(a,tmp);}  //用减法来模拟
	}
	while(c[0]>0&&c[c[0]]==0)c[0]--;
	return ;
}
int main()
{
  	memset(a,0,sizeof(a));
  	memset(b,0,sizeof(b));
  	memset(c,0,sizeof(c));
  	init(a);init(b);
  	chugao(a,b,c);
  	print(c);
  	print(a);
  	return 0;
}

例题六：回文数

题目描述

【问题描述】

若一个数（首位不为零）从左向右读与从右向左读都是一样，我们就将其称之为回文数。例如：给定一个 10进制数 56，将 56加 65（即把56从右向左读），得到 121是一个回文数。又如，对于10进制数87.

STEPl： 87＋78= 165

STEP2： 165＋561= 726

STEP3：726＋627＝1353

STEP4：1353+3531=4884

在这里的一步是指进行了一次N进制的加法，上例最少用了4步得到回文数4884。　　

写一个程序，给定一个N（2＜N＜＝10或N=16）进制数 M．求最少经过几步可以得到回文数。如果在30步以内（包含30步）不可能得到回文数，则输出“Impossible”

【输入样例】

9 87

【输出样例】

分析题目

N进制运算

1、当前位规范由%10改为% n

2、进位处理由/10改为/n

3、其他运算规则不变

标准程序

#include<iostream>
#include<cstring>
using namespace std;
int n,a[101],b[101],ans,i;
void init(int a[])                                  //将数串s转化为整数数组a
{ 
	string s; 
  	cin>>n>>s;                                   //读入字符串s 
  	memset(a,0,sizeof(a));                 //数组a清0 
  	a[0]=s.length();                             //用a[0]计算字符串s的位数 
  	for(i=1;i<=a[0];i++)
    	    if(s[a[0]-i]>='0'&&s[a[0]-i]<='9')   a[i]=s[a[0]-i]-'0';
         else a[i]=s[a[0]-i]-'A'+10;
}
bool check(int a[])                           //判别整数数组a是否为回文数
{ 
	for(i=1;i<=a[0];i++)
    		if(a[i]!=a[a[0]-i+1])return false;
  	return true;
}
void  jia(int a[])                    //整数数组a与其反序数b进行n进制加法运算
{
	for(int i=1;i<=a[0];i++)b[i]=a[a[0]-i+1];    //反序数b
	for(int i=1;i<=a[0];i++) a[i]+=b[i];             //逐位相加
	for(int i=1;i<=a[0];i++)                             //处理进位
	{a[i+1]+=a[i]/n;
		a[i]%=n;
	} 
	if(a[a[0]+1]>0) a[0]++;    //修正新的a的位数（a+b最多只能的一个进位） 
}
int main()
{	init(a);
	if(check(a)){cout<<0<<endl;return 0;}
	ans=0;                                               //步数初始化为0
	while(ans++<=30)
	{       jia(a);
		if(check(a)){cout<<ans<<endl;return 0;}
	}
	cout<<"Impossible";                           //输出无解信息
	return 0;
}