两点式直线方程找对称点，一般式不支持特殊位置直线（平行于x或者y轴的直线），下面是全部代码：
在这里插入代码片
```#-*- coding: utf-8 -*-
import matplotlib.pyplot as plt
import numpy as np
from pylab import * #让绘图显示中文的库
import sys

mpl.rcParams['font.sans-serif'] = ['SimHei']  #让绘图显示中文的命令

def caculate1(x1,x3,y3):
    """计算特殊情况下的直线对称点，输入的两点坐标X相同，即关于平行于Y轴的直线的对称点"""
    x4=2*x1-x3
    y4=y3
    plt.scatter([x4],[y4],color='blue',marker='*',label='所求对称点')
    plt.legend()
    plt.show()

def caculate2(x1,x3,y3):
    """计算特殊情况下的直线对称点，输入的两点坐标Y相同，即关于平行于X轴的直线的对称点"""
    x4=x3
    y4=2*y1-y3
    plt.scatter([x4],[y4],color='blue',marker='*',label='所求对称点')
    plt.legend()
    plt.show()

def caculate3(A,B,C,x3,y3):
    """计算一般情况的直线对称点，根据斜率关系推导的数学关系式"""
    x4=x3-2*A*((A*x3+B*y3+C)/(A*A+B*B))
    y4=y3-2*B*((A*x3+B*y3+C)/(A*A+B*B))
    plt.scatter([x4],[y4],color='blue',marker='*',label='所求对称点')
    plt.legend()
    plt.show()
    
while True:   
    start=eval(input("\n请选择直线方程形式(1 or 2)\n1 一般式\n2 两点式\n\n"))
    
#一般式还不能画特殊情况直线

    if start==3:
            break
            sys.exit()
    
    if start==1:

        plt.xticks([0,10,20,30,40,50,60,70,80,90,100])
        plt.yticks([0,10,20,30,40,50,60,70,80,90,100])
        plt.axis('scaled') #注意

        A,B,C=eval(input("请输入一般式方程三系数A,B,C(以逗号分隔):"))
        h=np.linspace(0,8,10)
        z=(-A*h-C)/B
        plt.plot(h,z,color='red',linewidth=2)
        
        x3,y3=eval(input("\n请输入已知任意对称点X,Y坐标:"))
        
        #plt.text(x3,y3,str((x3,y3)))
        
        plt.scatter([x3],[y3],color='red',marker='+')
        caculate3(A,B,C,x3,y3)
        
#两点式支持所有直线
        
    elif start==2:
        plt.xticks([0,10,20,30,40,50,60,70,80,90,100])
        plt.yticks([0,10,20,30,40,50,60,70,80,90,100])
        plt.xlabel("X")
        plt.ylabel("Y")
        plt.title(" 找对称点")
        plt.axis('scaled') #注意

        x1,y1=eval(input("\n请输入X1,Y1坐标(以逗号分隔下同):"))
        x2,y2=eval(input("请输入X2,Y2坐标:"))
        
        x3,y3=eval(input("\n请输入已知任意对称点X,Y坐标:"))

        if x1==x2:
            plt.plot([x1,x2],[y1,y2],color='red',linewidth=2)
            plt.scatter([x3],[y3],color='red',marker='+',label='已知点')
            plt.legend()
            caculate1(x1,x3,y3)
        elif y1==y2:
            plt.plot([x1,x2],[y1,y2],color='red',linewidth=2)
            plt.scatter([x3],[y3],color='red',marker='+',label='已知点')
            plt.legend()
            caculate2(x1,x3,y3)
        else:
            plt.plot([x1,x2],[y1,y2],color='red',linewidth=2)
            plt.scatter([x3],[y3],color='red',marker='+',label='已知点')
            plt.legend()
            A=y1-y2
            B=x2-x1
            C=x1*y2-y1*x2
            caculate3(A,B,C,x3,y3)
    else:
        print("请选择输入1或者2")

#include <stdio.h>
 int main()
{
    int n;
    float a,b,A,B,C,a1,b1;
    scanf("%d\n",&n);
    while(n--)
    {  scanf("%f %f %f %f %f",&a,&b,&A,&B,&C);
     int a1=int (a-2*A*(A*a+B*b+C)/(A*A+B*B));
     int b1=int (b-2*B*(A*a+B*b+C)/(A*A+B*B));
     printf("%d %d\n",a1,b1);
    }
}
Ax+By+C=0
double dist(POINT p1,POINT p2)                // 返回两点之间欧氏距离 
{ 
	return( sqrt( (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y) ) ); 
}

Problem Description
The light travels in a straight line and always goes in the minimal path between two points, are the basic laws of optics.


Now, our problem is that, if a branch of light goes into a large and infinite mirror, of course，it will reflect, and leave away the mirror in another direction. Giving you the position of mirror and the two points the light goes in before and after the reflection,
 calculate the reflection point of the light on the mirror.

  

You can assume the mirror is a straight line and the given two points can’t be on the different sizes of the mirror.
Input
The first line is the number of test case t(t<=100).

  

The following every four lines are as follow:

  X1 Y1

  X2 Y2

  Xs Ys

  Xe Ye


  (X1,Y1),(X2,Y2) mean the different points on the mirror, and (Xs,Ys) means the point the light travel in before the reflection, and (Xe,Ye) is the point the light go after the reflection.


  The eight real number all are rounded to three digits after the decimal point, and the absolute values are no larger than 10000.0.
Output
  Each lines have two real number, rounded to three digits after the decimal point, representing the position of the reflection point.
Sample Input

1
0.000 0.000
4.000 0.000
1.000 1.000
3.000 1.000

Sample Output

2.000 0.000 

Source
2009 Multi-University Training Contest 5 - Host by NUDT 
 
模板：
const double eps=1e-6;
struct point
{
    double x,y;
};
struct line//直线ax+by+c=0
{
    double a,b,c;
};
line LineFromSegment(point p1,point p2)//两个点求直线
{
    line temp;
    temp.a=p2.y-p1.y;
    temp.b=p1.x-p2.x;
    temp.c=p2.x*p1.y-p1.x*p2.y;
    return temp;
}
point LineInter(line l1,line l2)//两条直线求交点
{
    point temp;
    double a1=l1.a;
    double b1=l1.b;
    double c1=l1.c;
    double a2=l2.a;
    double b2=l2.b;
    double c2=l2.c;
    if(fabs(b1)<eps)
    {
        temp.x=-c1/a1;
        temp.y=(-c2-a2*temp.x)/b2;
    }
    else
    {
        temp.x=(c1*b2-b1*c2)/(b1*a2-b2*a1);
        temp.y=(-c1-a1*temp.x)/b1;
    }

    return temp;
}
point symmetrical(point p, line L)//求一个点关于一条直线的对称点
{
    point p2;
    double d;
    d = L.a * L.a + L.b * L.b;
    p2.x = (L.b * L.b * p.x - L.a * L.a * p.x -
            2 * L.a * L.b * p.y - 2 * L.a * L.c) / d;
    p2.y = (L.a * L.a * p.y - L.b * L.b * p.y -
            2 * L.a * L.b * p.x - 2 * L.b * L.c) / d;
    return p2;
}


本题代码：
#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
const double eps=1e-6;
struct point
{
    double x,y;
};
struct line//直线ax+by+c=0
{
    double a,b,c;
};
line LineFromSegment(point p1,point p2)//两个点求直线
{
    line temp;
    temp.a=p2.y-p1.y;
    temp.b=p1.x-p2.x;
    temp.c=p2.x*p1.y-p1.x*p2.y;
    return temp;
}
point LineInter(line l1,line l2)//两条直线求交点
{
    point temp;
    double a1=l1.a;
    double b1=l1.b;
    double c1=l1.c;
    double a2=l2.a;
    double b2=l2.b;
    double c2=l2.c;
    if(fabs(b1)<eps)
    {
        temp.x=-c1/a1;
        temp.y=(-c2-a2*temp.x)/b2;
    }
    else
    {
        temp.x=(c1*b2-b1*c2)/(b1*a2-b2*a1);
        temp.y=(-c1-a1*temp.x)/b1;
    }

    return temp;
}
point symmetrical(point p, line L)//求一个点关于一条直线的对称点
{
    point p2;
    double d;
    d = L.a * L.a + L.b * L.b;
    p2.x = (L.b * L.b * p.x - L.a * L.a * p.x -
            2 * L.a * L.b * p.y - 2 * L.a * L.c) / d;
    p2.y = (L.a * L.a * p.y - L.b * L.b * p.y -
            2 * L.a * L.b * p.x - 2 * L.b * L.c) / d;
    return p2;
}
int main()
{
    int t;
    cin>>t;
    for(int i=1;i<=t;i++)
    {
        point p[4];
        for(int j=0;j<4;j++)
            cin>>p[j].x>>p[j].y;
        point tmp;
        line l1,l2;
        l1=LineFromSegment(p[0],p[1]);//两点求直线
        tmp=symmetrical(p[3],l1);//直线外一点与一条直线，求该点关于直线的对称点,这里是求的反射点的对称点
        l2=LineFromSegment(tmp,p[2]);//对称点与入射点求直线
        tmp=LineInter(l1,l2);//求反射点（两条直线相交的点）
        cout<<setiosflags(ios::fixed)<<setprecision(3)<<tmp.x<<" "<<tmp.y<<endl;
    }
    return 0;
}



辅助图（根据所写代码画的）

链接：https://www.nowcoder.com/acm/contest/86/B
来源：牛客网

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
Special Judge, 64bit IO Format: %lld
题目描述
   开始涉猎几何领域了。他现在正在研究小马喝水问题。 
 
  众所周知，这个问题中有一匹口渴的小马，一条笔直的河，以及小马的家。小马需要去河边喝水，然后再去家里。它需要走最短的路径。 
 
  
 
  解决这个问题也很简单，其中有一个步骤是要做小马家关于河水的对称点。 
 
   正对此感到一些烦恼。他不会做这个。他想请你帮他作一条过小马家且垂直于河水的线，然后告诉  垂足的位置。
输入描述:
第一行一个整数  ，表示  的询问个数。
接下去  行，每行  个实数 ，表示小马家在点  ，河水为直线 
输出描述:
输出共  行，每行两个实数 ， 表示答案垂足点的坐标 。 
当你的答案与标准输出的误差小于  时，视为答案正确。