Giải tích tổ hợp

Lập trình tính các công thức của giải tích tổ hợp:
       - Số hoán vị của các phần tử P(n)=n!
       - Số chỉnh hợp chập k từ n phần tử A(k,n)= n!/(n-k)!
       - Số tổ hợp chập k từ n phần tử C(k,n)= n!/(k!(n-k)!)
       - Tổng số các tập con của tập có n phần tử =C(0,n)+…+C(b,n)
   Với k,n là các số tự nhiên nhập từ bàn phím.
                                                                 Giải
Var  n,k: Integer; P,A,C,T:Real;
Function GT(n:Integer):Real;
   Var  i: Integer; tmp: Real;
   Begin
       tmp:=1;
       For i:=1 to n Do Tmp:=tmp*i;
       GT:=tmp;
   End;
BEGIN
       Write(' Cho n, k '); Readln(n,k);
       P:=Gt(n); A:=Gt(n)/gt(n-k); C:=gt(n)/(gt(n-k)*gt(k)); T:=0;
       For k:=0 to n Do T:=T+gt(n)/(gt(n-k)*gt(k));
       Writeln('=> Hoan vi cua ',n,' phan tu la P = ',p:12:0);
       Writeln('=> Chinh hop chap ',k,' cua ',n,' phan tu la A = ',A:12:0);
       Writeln('=> To hop chap ',k,' cua ',n,' phan tu la C = ',c:12:0);
       Writeln('=> Tong cua cac to hop chap k cua n phan tu(k chay tu 0 den ',n,') la T = ',t:9:0);
       Readln;
END.