计算sin(20.86^100)

计算sin(20.86^100)

思路

首先20.86¹⁰⁰太大了，所以 考虑用高精度乘法来做。
但是浮点数的高精度处理起来比较麻烦，所以想到可以先把20.86¹⁰⁰变为2086¹⁰⁰，这样就是小数点向后移动了两百位。

接下来由于sinx是周期函数，周期为2π，小数点也同时后移两百位，所以需要模2π*10²⁰⁰。这里我选择直接上网搜了π的小数点后两百位，直接定义了一个字符串PI。

#include<bits/stdc++.h>
using namespace std;

inline int compare(string s1, string s2){
	if(s1.length()>s2.length()) return 1;
	else if(s1.length()<s2.length()) return -1;
	else return s1.compare(s2); 
}//s1>s2返回1,s1<s2返回-1,s1==s2返回0

//以下是用了一些高精度的模板
//⾼精度减法
string Sub(string s1, string s2){ 
	string str;
    int sign=1;
	int res=compare(s1, s2); 
    if(res==0)return "0";
	if(res<0){
		sign=-1;
		string tmp=s1;
        s1=s2;
		s2=tmp;
    }
	string::size_type len; 
    len=s1.size()-s2.size();
	for(int i=s2.size()-1;i>=0;i--) {
		if (s1[i+len]<s2[i]) {
			s1[i+len-1]=char(int(s1[i+len-1])-1); 
            str=char(s1[i+len]-s2[i]+58)+str;
		}
		else str=char(s1[i+len]-s2[i]+48)+str;
	}
    for (int i=len-1;i>=0;i--) str=s1[i]+str;

	str.erase(0,str.find_first_not_of('0')); 
    if (str.empty()) str="0";
	if ((sign==-1)&&(str[0]!='0')) str="-"+str;
	return str;

}

//⾼精度加法
string Add(string s1, string s2){
    int sign=1; 
	string str; 
    if(s1[0]=='-'){
		if(s2[0]=='-') {
			sign=-1;
			str=Add(s1.erase(0,1),s2.erase(0,1));
		}
		else str=Sub(s2,s1.erase(0,1));
	}
	else {
		if (s2[0]=='-')
            str=Sub(s1,s2.erase(0,1)); 
        else {
			string::size_type l1, l2; 
			l1=s1.size(); 
            l2=s2.size(); 
            if (l1<l2) 
                for(int i=1;i<=l2-l1;i++) s1="0"+s1;
			else 
                for(int i=1;i<=l1-l2;i++) s2="0"+s2;
			int int1=0,int2=0; //int2 记录进位
			for(int i=s1.size()-1;i>=0;i--) {
				int1=(s1[i]-'0'+s2[i]-'0'+int2)%10;
				int2=(s1[i]-'0'+s2[i]-'0'+int2)/10; 
                str=char(int1+48)+str;
			}
			if(int2!=0) str=char(int2+48)+str;
		}
	}
	if ((sign == -1)&&(str[0]!='0')) str="-"+str;
	return str;
}

//⾼精度乘法
string Mul(string s1, string s2)
{
	string str;
	string::size_type l1, l2;
	l1=s1.size(); 
    l2=s2.size();
	for(int i=l2-1;i>=0;i--) {  
		string tmp;
		int int1=0,int2=0,int3=s2[i]-'0';
		if (int3!=0) {
			for(int j=1;j<=(int)(l2-1-i);j++)
                tmp="0"+tmp;
			for(int j=l1-1;j>=0;j--) {
				int1=(int3*(s1[j]-'0')+int2)%10;
				int2=(int3*(s1[j]-'0')+int2)/10;
				tmp=char(int1+48)+tmp;
			}
			if(int2!=0)tmp=char(int2+48)+tmp;
		}
		str=Add(str,tmp);
	}
	str.erase(0,str.find_first_not_of('0'));
	if(str.empty()) str="0";
	return str;
}

//⾼精度除法求余数
string Div(string s1,string s2){
    string r;//余数
    int f=compare(s1,s2);
    if(f<0)r=s1;
    else if(f==0) r="0";
    else{
        string::size_type len1,len2;
        len1=s1.size();
        len2=s2.size();
        string tmp;
        tmp.append(s1,0,len2-1);
        for(int i=len2-1;i<len1;i++){
            tmp=tmp+s1[i];
            for(char ch='9';ch>='0';ch--){
                string str;
                str=str+ch;
                if(compare(Mul(s2,str),tmp)<=0){
                    tmp=Sub(tmp,Mul(s2,str));
                    break;
                }
            }
        }
        r=tmp;
    }
    return r;
}

void trans(string s){
    char c[400];
	strcpy(c,s.c_str());
	double ans=0;
	if (s.length()==200){ //是否有整数部分
		for(int i=0;i<16;i++) ans+=((c[i]-48)/pow(10,i+1));
	}else{
		for(int i=1;i<=16;i++) ans+=((c[i]-48)/pow(10,i));
		ans+=(c[0]-48);
	}
	printf("%.16lf",sin(ans));
}

int main() {
	string s0="2086";
	string st="2086";
	for (int i=0;i<99;i++) 
        st=Mul(st,s0);
    
    //π取小数点后两百位
	string PI="314159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196";
	string x="2";
	PI=Mul(PI,x);
	st=Div(st,PI);

    trans(st);//转化成浮点数并输出
    return 0;
}