UVa530 - Binomial Showdown

#include <stdio.h>

long gcd(long a, long b);
void divideGcd(long *a, long *b);
long c(long n, long k);

int main()
{
    long n, m;

    while (scanf("%ld%ld", &n, &m) && (n + m))
    {
        printf("%ld\n", c(n, m));
    }

    return 0;
}

long gcd(long a, long b)
{
    return b == 0?a:gcd(b, a % b);
}

void divideGcd(long *a, long *b)
{
    long g = gcd(*a, *b);
    *a = *a / g;
    *b = *b / g;
}

long c(long n, long k)
{
    long numerator = 1, denominator = 1, toMul, toDiv, i;

    if (k > n / 2)
        k = n - k;

    for (i = k; i > 0; i--)
    {
        toMul = n - k + i;
        toDiv = i;
        divideGcd(&toMul, &toDiv);
        divideGcd(&numerator, &toDiv);
        divideGcd(&toMul, &denominator);
        numerator *= toMul;
        denominator *= toDiv;
    }

    return numerator / denominator;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

kgduu

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值