Để tối giản phân số M/N, bạn cần tìm ước số chung lớn nhất (GCD) của M và N, và sau đó chia cả hai số cho GCD đó. Dưới đây là mã code Python để thực hiện điều này:

```python
def gcd(a, b):
while b != 0:
a, b = b, a % b
return a

def simplify_fraction(M, N):
g = gcd(M, N)
P = M // g
Q = N // g
return P, Q

def main():
M = int(input("Nhập vào số M (nguyên dương, nhỏ hơn 1000): "))
N = int(input("Nhập vào số N (nguyên dương, nhỏ hơn 1000): "))

if 0 < M < 1000 and 0 < N < 1000:
P, Q = simplify_fraction(M, N)
print(f"Tử số P là: {P}, Mẫu số Q là: {Q}")
else:
print("Vui lòng nhập M và N là số nguyên dương và nhỏ hơn 1000.")

if __name__ == "__main__":
main()
```

### Hướng dẫn sử dụng:
1. Chạy đoạn mã trên.
2. Nhập vào hai số nguyên dương M và N (cả hai đều phải nhỏ hơn 1000).
3. Kết quả sẽ hiển thị tử số P và mẫu số Q của phân số M/N sau khi đã tối giản.

### Giải thích mã:
- Hàm `gcd(a, b)` sử dụng thuật toán Euclid để tính ước số chung lớn nhất của `a` và `b`.
- Hàm `simplify_fraction(M, N)` sẽ gọi hàm `gcd` để tìm GCD và sau đó chia M và N cho GCD để có phân số đã tối giản.
- Chương trình chính (`main()`) nhận đầu vào và gọi hàm tối giản, đồng thời kiểm tra điều kiện đầu vào.