洛谷 p1625

高精度

我以为这题必有高论，怎么想也想不出来，没想到竟是如此粗鄙做法。

我们写一个高精度模拟一下，然后枚举约数看是否能约分，由于我不会高精度除法，就抄了一发

其实这种两项之比和项数有关的数列是不能推通项的，只能暴力模拟

#include<bits/stdc++.h>
using namespace std;
const int N = ;
int n, m;
struct Bigint {
    int a[N], len;
    Bigint friend operator + (Bigint &A, Bigint &B)
    {
        Bigint C;
        memset(C.a, , sizeof(C.a));
        int lim = max(A.len, B.len);
        for(int i = ; i <= lim; ++i)
        {
            C.a[i] += A.a[i] + B.a[i];
            C.a[i + ] += C.a[i] / ;
            C.a[i] %= ;
        }
        C.len = lim;
        if(C.a[lim + ]) ++C.len;
        return C;
    }
    Bigint friend operator * (Bigint A, int x)
    {
        for(int i = ; i <= A.len; ++i) A.a[i] *= x;
        for(int i = ; i <= A.len; ++i)
        {
            A.a[i + ] += A.a[i] / ;
            A.a[i] %= ;
        }
        while(A.a[A.len + ] >= )
        {
            ++A.len;
            A.a[A.len + ] += A.a[A.len] / ;
            A.a[A.len] %= ;
        }
        if(A.a[A.len + ] > ) ++A.len;
        return A;
    }
    Bigint friend operator / (Bigint &A, int x)
    {
        Bigint B;
        int tot = ;
        memset(B.a, , sizeof(B.a));
        for(int i = A.len; i; --i)
        {
            tot = tot *  + A.a[i];
            B.a[i] = tot / x;
            tot %= x;
        }
        B.len = A.len;
        while(B.a[B.len] == ) --B.len;
        return B;
    }
    void print()
    {
        for(int i = len; i; --i) printf("%d", a[i]);
        puts("");
    }
} down, up, base;
bool can_div(Bigint &A, int x)
{
    int tot = ;
    for(int i = A.len; i; --i) tot = (tot *  + A.a[i]) % x;
    return (tot == );
}
int main()
{
    scanf("%d%d", &n, &m);
    down.len = ;
    down.a[] = ;
    base.len = ;
    base.a[] = ;
    for(int i = m + ; i < n + m; ++i) base = base * i;
    for(int i = ; i <= n; ++i)
    {
        up = up + base;
        base = base * i;
        base = base / (m + i);
    }
    for(int i = ; i < n + m; ++i)
    {
        if(can_div(up, i)) up = up / i;
        else down = down * i;
    }
    for(int i = ; i < n + m; ++i) if(can_div(up, i) && can_div(down, i))
    {
        up = up / i;
        down = down / i;
    }
    up.print();
    down.print();
    return ;
}