用python写程序：有算式ABCD-CDC=ABC，其中A、B、C、D均为非负整数，编写程序，求A、B、C、D的值。

这个问题实际上是一个线性代数方程组的问题，可以用Python中的数学库（如numpy或scipy）来解决。我们首先可以将算式转化为矩阵形式，然后利用解线性方程组的方法来找到A、B、C、D的值。以下是使用numpy的一个简单示例：

import numpy as np

# 设定已知等式：[A, B, C] * [D, -1, -1] = [A, B, C]
# 这里假设A、B、C已经给定，例如A=10, B=20, C=30
known_values = [10, 20, 30]

# 创建系数矩阵和常数向量
coeff_matrix = np.array([[1, 0, -1], [0, 1, -1]])
constant_vector = known_values

# 使用numpy的linalg.solve函数求解线性方程
solution = np.linalg.solve(coeff_matrix, constant_vector)

# A, B, D的值就是解
A, B, D = solution

# 根据题目描述，D是非负整数，我们需要取解中的最大非负整数部分
D = max(0, int(D))

# 现在我们需要计算C，因为C = (ABCD - CDC) / A
# 我们先尝试直接用已知的A和D，如果结果不是整数，则需要调整B，使得C也是非负整数
try:
    C = (known_values[0]*D - known_values[2]) // known_values[0]
except TypeError:
    # 如果除法出错，说明B需要调整
    B -= 1  # 减一试一下
    C = (known_values[0]*D - known_values[2]) // known_values[0]
    
assert C >= 0 and isinstance(C, int), "无法得到非负整数解"

(A, B, C, D)