Перевести код из СИ в Python

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <windows.h>
char bufRus[256];
char* ru(const char*text)
{
    CharToOem(text, bufRus);
    return bufRus;
}
 
int *input_matrix(int *n, int *m, char c);
void output_matrix(int *p, int n, int m, char c);
int sum(int *p, int n, int m);
int equation(int a, int b, int c);
 
int main()
{
    int *Y, *P, *Q;
    int ny, my, np, mp, nq, mq;
    int a, b, c, D;
    Y = input_matrix(&ny, &my, 'Y');
    P = input_matrix(&np, &mp, 'P');
    Q = input_matrix(&nq, &mq, 'Q');
    output_matrix(Y, ny, my, 'Y');
    output_matrix(P, np, mp, 'P');
    output_matrix(Q, nq, mq, 'Q');
    a = sum(Y, ny, my);
    b = sum(P, np, mp);
    c = sum(Q, nq, mq);
    printf(ru("Сумма элементов, лежащих ниже главной диагонали и на ней, матрицы Y, a = %d\n"), a);
    printf(ru("Сумма элементов, лежащих ниже главной диагонали и на ней, матрицы P, b = %d\n"), b);
    printf(ru("Сумма элементов, лежащих ниже главной диагонали и на ней, матрицы Q, c = %d\n"), c);
    D = equation(a, b, c);
    printf("D = %d", D);
    fflush(stdin);
    getchar();
    return (0);
}
 
//Вспомогательный алгоритм ввода массивов
int *input_matrix(int *n, int *m, char c)
{
     int i, j;
     int *p;
     printf(ru("Введите число строк и столбцов матрицы %c (n и m): "),c);  
     scanf("%d%d", n, m); 
     p=new int[*n**m];     
     for(i=0; i<*n; i++)
           for(j=0; j<*m; j++)
               {
                    printf("%c[%d][%d]=", c, i+1, j+1);
                    scanf("%d", &p[i**m+j]);
               }
     return (p);
}
 
//Вспомогательный алгоритм вывода массивов
void output_matrix(int *p, int n, int m, char c)
{
    int i, j;
    printf(ru("\tМатрица %c\n"), c);
    for(i=0; i<n; i++)
        {
            for(j=0; j<m; j++)
                   printf("%6d ",p[i*m+j]);
            printf("\n");
        }
 } 
 
 
//Вспомогательный алгоритм вычисления суммы элементов, лежащих ниже и на главной диагонали
int sum(int *p, int n, int m)
{ 
    int sum = 0;
    int i, j;
    for(i = 0; i < n; i++)
         for(j = 0; j < m; j++)
             if(i >= j) 
                sum += p[i*m+j];
    return sum;
}
 
int equation(int a, int b, int c)
{
    int D;
    float x, x1, x2;
    D = b*b-4*a*c;
    if(D>0)
       {
           x1 = (-b + sqrt(D))/(2*a);
           x2 = (-b - sqrt(D))/(2*a);
           printf("x1 = %f, x2 = %f ", x1, x2);
       }
    else if(D==0)
            {
                x=-b/(2*a);
                printf("x = %f ", x);
            }
         else
            printf(ru("D < 0 - нет корней "));
    return D;
}