计算sin(20.86^100)
思路
首先20.86100太大了,所以 考虑用高精度乘法来做。
但是浮点数的高精度处理起来比较麻烦,所以想到可以先把20.86100变为2086100,这样就是小数点向后移动了两百位。
接下来由于sinx是周期函数,周期为2π,小数点也同时后移两百位,所以需要模2π*10200。这里我选择直接上网搜了π的小数点后两百位,直接定义了一个字符串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;
}