两条直线的夹角

编写一个c++程序，计算四个点所构成的两条直线的夹角
setprecision(n)描述：可以控制输出流显示浮点数的数字个数。C++默认的输出流数值有效位是6。
平面上给定4个点的坐标A B C D，分别表示直线AB和CD。坐标均为绝对值不超过100的整数。计算出这两条直线所成的锐角是多少度。（保留两位小数）
注意：当两直线平行或者重合，答案为0。
测试举例：
测试输入：0 0 1 0 0 0 0 1
预期输出：90.00
测试输入：16 34 21 87 98 23 12 9
预期输出：75.36

#include<iostream>
#include<cmath>
#include<iomanip>
#define pi 3.1415926;
using namespace std;
int main()
{    
    double x1,y1,x2,y2,x3,y3,x4,y4;
    long double x,y;
    cin>>x1>>y1>>x2>>y2>>x3>>y3>>x4>>y4;
    x=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
    y=(x3-x4)*(x3-x4)+(y3-y4)*(y3-y4);
    long double s1,s2;
    s1=sqrt(x);
    s2=sqrt(y);
    long double s3;
    s3=(x2-x1)*(x4-x3)+(y2-y1)*(y4-y3);
    long double s4;
    s4=abs(s3);
    long double ans,ans2;
    ans2=s4/(s1*s2);
    ans=acos(ans2);
    long  double ans3;
    ans3=(180*ans)/pi;
    cout<<fixed<<setprecision(2)<<(long double)ans3;
    return 0;
}