上两道省选的高精吧！

T1HAOI2016 放棋子

题目描述

给你一个N*N的矩阵，每行有一个障碍，数据保证任意两个障碍不在同一行，任意两个障碍不在同一列，要求你在这个矩阵上放N枚棋子（障碍的位置不能放棋子），要求你放N个棋子也满足每行只有一枚棋子，每列只有一枚棋子的限制，求有多少种方案。

输入输出格式

输入格式：

第一行一个N，接下来一个N*N的矩阵。N<=200,0表示没有障碍，1表示有障碍，输入格式参考样例

输出格式：

一个整数，即合法的方案数。

输入输出样例

输入样例#1：

2
0 1
1 0

输出样例#1：

1


额，就是递推的错排，证明网上一大坨。。。难打的是高精。。。

#include<cstdio>
#include<iostream>
#include<iomanip>
#define maxn 5050
#define mod 10000000000
#define _ 10
#define ll long long
using namespace std;
struct bignum{
    int sz;
    ll v[maxn];
}s[3100];
typedef bignum bg;
int n;

inline bg operator +(bg a , bg b)
{
    if(a.sz < b.sz) swap(a,b);
    ll jin = 0;
    for(int i=1;i<=a.sz;i++)
    {
        a.v[i] += b.v[i] + jin;
        jin = a.v[i] / mod;
        a.v[i] %= mod;
    }
    while(jin)
    {
        a.v[++a.sz] = jin % mod;
        jin /= mod;
    }
    return a;
}

inline bg operator *(bg a , ll b)
{
    ll jin = 0;
    for(int i=1;i<=a.sz;i++)
    {
        a.v[i] = a.v[i] * b + jin;
        jin = a.v[i] / mod;
        a.v[i] %= mod;
    }
    while(jin)
    {
        a.v[++a.sz] = jin % mod;
        jin /= mod;
    }
    return a;
}

inline ostream& operator << (ostream &os,bg x)
{
    ll i;
    os<<x.v[x.sz];
    for(i=x.sz-1;i;i--)
        os<<setfill('0')<<setw(_)<<x.v[i];
    return os;
}

int main(){
    cin >> n;
    s[1].sz = 1;
    s[1].v[1] = 0;
    s[2].sz = 1;
    s[2].v[1] = 1;    
    for(int i=3;i<=n;i++)
        s[i] = (s[i - 1] + s[i - 2]) * (i - 1);
    cout << s[n];
}

T2HNOI2004 高精度开根

晓华所在的工作组正在编写一套高精度科学计算的软件，一些简单的部分如高精度加减法、乘除法早已写完了，现在就剩下晓华所负责的部分：实数的高精度开m次根。

因为一个有理数开根之后可能得到一个无理数，所以这项工作是有较大难度的。现在要做的只是这项工作的第一步：只对自然数进行开整数次根，求出它的一个非负根，并且不考虑结果的小数部分，只要求把结果截断取整即可。

程序需要根据给定的输入，包括需要开根的次数，以及被开根的整数；计算出它的非负根取整后的结果。

输入输出格式

输入格式：

输入文件共有两行，每行都有一个整数，并且输入中没有多余的空格：

第一行有一个正整数m (1 <= m <= 50)，表示要开的根次；

第二行有一个整数n (0<=n <= 10^10000)，表示被开根的数。

输出格式：

输出文件只有一行，包括一个数，即为开根取整后的结果。

输入输出样例

输入样例#1：

 3
1000000000

输出样例#1：

1000


靠，这个题目网上的奆佬们都用PY，，，，QAQ，，，，不会PY的小蒟蒻难受啊。只有打一波高精了...TAT...

#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<iomanip>
#define maxn 10050
#define mod 100000000
#define _ 8
#define ll long long
using namespace std;
struct bg{
    int sz;
    ll v[maxn];
    bg()
    {
        memset(v,0,sizeof(v));
    }
};

bg ans;
int n;
char shu[100100];

inline ostream& operator << (ostream &os,bg x)
{
    ll i;
    os<<x.v[x.sz];
    for(i=x.sz-1;i;i--)
        os<<setfill('0')<<setw(_)<<x.v[i];
    return os;
}

inline bool operator <(bg a , bg b)
{
    if(a.sz > b.sz) return false;
    if(a.sz < b.sz) return true;
    for(int i=a.sz;i>=0;--i)
    {
        if(a.v[i] < b.v[i]) return true;
        if(a.v[i] > b.v[i]) return false;
    }
    return false;
}

inline bg operator +(bg a , bg b)
{
    if(a.sz < b.sz) swap(a , b);
    ll jin = 0;
    for(int i=1;i<=a.sz;++i)
    {
        ll val = a.v[i] + b.v[i] + jin;
        a.v[i] = val % mod;
        jin = val / mod;
    }
    while(jin)
    {
        a.v[++a.sz] = jin % mod;
        jin /= mod;
    }
    return a;
}

inline bg operator -(bg a , int b)
{
    for(int i=a.sz;i>1;--i)
    {
        --a.v[i];
        a.v[i-1] += mod;
    }
    a.v[1]--;
    ll jin = 0;
    for(int i=1;i<=a.sz;++i)
    {
        a.v[i] += jin;
        jin = a.v[i] / mod;
        a.v[i] %= mod;
    }
    return a;
}

inline bg operator *(bg a , bg b)
{
    bg c;
    for(int i=1;i<=a.sz;++i)
        for(int j=1;j<=b.sz;++j)
                c.v[i + j - 1] += a.v[i] * b.v[j] ;
    ll ce = a.sz + b.sz - 1,jin = 0;
    for(int i=1;i<=ce;++i)
    {
        c.v[i] += jin;
        jin = c.v[i] / mod;
        c.v[i] %= mod;
    }
    while(jin)
    {
        c.v[++ce] = jin % mod;
        jin /= mod;
    }
    c.sz = ce;
    return c;
}

inline bg operator /(bg a , int b)
{
    for(int i=a.sz;i>=1;--i)
    {
        if(a.v[i] & 1) a.v[i - 1] += mod;
        a.v[i] >>= 1;
    }
    a.v[0] = 0;
    if(!a.v[a.sz]) --a.sz;
    return a;
}

inline bg read()
{
    bg ans;
    scanf("%s",shu + 1);
    int len = strlen(shu + 1);
    int sq = len;
    int qwer = 0;
    for(int i=1;;++i)
    {
        if(qwer)
        {
            ans.sz = i - 1;
            break;
        }
        for(int j=0;j<_;++j)
        {
            ans.v[i] += (shu[sq--] - '0') * pow(10 , j);
            if(sq < 1)
            {
                qwer = 1;
                break;
            }
        }
    }
    return ans;  
}

inline bg find(bg a)
{
    bg k;
    k.sz = 1;
    k.v[1] = 1;
    int ce = n;
    while(ce)
    {
        if(ce & 1) k = k * a;
        a = a * a;
        ce >>= 1;
    }
    return k;
}

int main(){
    cin >> n;
    getchar();
    ans = read();
    
    bg l;
    bg r;
    r.sz = 1;
    r.v[1] = 1;
    
    bg ce;
    ce.sz = 1;
    ce.v[1] = 1;
    
    bg c;
    c.sz = 1;
    c.v[1] = 2;
    
    while(find(r) < ans)
    {
        l = r;
        r = r * c;
    }
    while(l < r)
    {
        bg mid = (l + r) / 2;
        if(find(mid) < ans) l = mid + ce;
        else r = mid;
    }
    
    bg k = l;
    for(int i=1;i<n;++i)
     k = k * l;
    if(ans < k)
        l = l - 1;
    cout << l;
}

在BZOJ上是可以轻松跑过的，，，但是洛谷里面被卡了（洛谷里面时限是1.5s，BZOJ5s，程序2.7s）...无语了。。。