Dwumiany - SPOJ

#include <iostream>
#include <iomanip>

using namespace std;

int t, liczbaA, liczbaB;

int main()
{
    cin >> t;
    for(int i=1; i<=t;i++)
    {
        cin >> liczbaA >> liczbaB;
        long double silniaA=1;
        long double silniaB=1;
        long double silniaC=1;
        if(liczbaA==0) silniaA=1;
        else
        {
            for (int j=1; j<=liczbaA;j++)
            {
                silniaA = silniaA * j;
            }
        }
        if(liczbaB==0) silniaB=1;
        else
        {
            for (int k=1; k<=liczbaB;k++)
            {
                silniaB = silniaB * k;
            }
        }
        if ((liczbaA==0) && (liczbaB == 0))silniaC=1;
        else
        {
           for (int l=1; l<=liczbaA-liczbaB;l++)
            {
                silniaC = silniaC * l;
            }
        }

        long double wynik = silniaA/(silniaB*silniaC);
        cout << setprecision(10000);
        cout << wynik << endl;
    }

    return 0;
}

import std.algorithm;
import std.bigint;
import std.conv;
import std.format;
import std.functional;
import std.math;
import std.stdio;
import std.string;


auto binomialCoefficientImpl(T,U)(U k, U n)
{
	T result = 1;
	auto loopMax = min(k-n, n);
	for(U i = 0; i < loopMax; ++i){
		result *= k-i;
		result /= i+1;
	}
	return result;
}

alias binomialCoefficient = memoize!(binomialCoefficientImpl!(ulong, uint));

auto extract(T, R)(R r)
{
	T t;
	if(r.readf(" %s", &t) != 1){
		throw new Error("readf failed");
	}
	//"extracted: '%s'".format(t).writeln;
	return t;
}

void main()
{
	int numberOfTests = stdin.extract!uint;
	for(int i = 0; i < numberOfTests; ++i) {
		int k, n;
		k = stdin.extract!uint;
		n = stdin.extract!uint;
		binomialCoefficient(k, n).writeln;
	}
}