OpenJudge9278:旅行

总时间限制: 

10000ms

单个测试点时间限制: 

1000ms

内存限制: 

131072kB

描述

转眼毕业了，曾经朝夕相处的同学们不得不都各奔东西，大家都去了不同的城市开始新的生活。在各自城市居住了一段时间后，他们都感到了一些厌倦，想去看看其他人的生活究竟如何，于是他们都选择到另一个同学所在城市去旅游，并且希望旅游的城市各不相同，他们想知道有多少种不同的方案，可是数量实在太多了，他们无法计算出来，你能帮助他们吗。

输入

一个正整数n(n<200)，表示人数。

输出

一个数，表示有多少不同的方案。

样例输入

3

样例输出

2

提示

有如下两种方案：

同学1去同学2的城市，同学2去同学3的城市，同学3去同学1的城市

同学1去同学3的城市，同学3去同学2的城市，同学2去同学1的城市

题解：错排+高精

错排递推式：f[i]=(n-1)*(f[i-1]+f[i-2])

最后一个元素n，如果到第k个元素的位置，有n-1种方法

k这个元素到第n个元素原来的位置，那么剩下的是f[i-2]，否则f[i-1]。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#define B 10000
using namespace std;

int n;

struct Big{
    int a[], n;
    int& operator [](int x) {return a[x];}
    Big():n() {memset(a, , sizeof(a));}
    void ini(int x) {a[]=x; n=;}
}ans[];

Big operator *(Big a, int b) {
    int g=;
    for(int i=; i<=a.n; i++)
        g += a[i]*b, a[i] = g%B, g/=B;
    if(g) a[++a.n] = g;
    return a;
}

Big operator *(Big a, Big b) {
    Big c;
    for(int i=; i<=a.n; i++) {
        int g=;
        for(int j=; j<=b.n; j++)
            g += c[i+j-]+a[i]*b[j], c[i+j-] = g%B, g/=B;
        c[i+b.n] = g;
    }
    c.n = a.n + b.n;
    while(c.n> && c[c.n]==) c.n--;
    return c;
}

Big operator +(Big a, Big b) {
    int g=, n=max(a.n, b.n);
    for(int i=; i<=n; i++) {
        g += i<=a.n ? a[i] : ;
        g += i<=b.n ? b[i] : ;
        a[i] = g%B, g/=B;
    }
    a.n = n;
    if(g) a[++a.n] = g;
    return a;
}

Big operator -(Big a, Big b) {
    for(int i=; i<=b.n; i++) {
        if(a[i]<b[i]) a[i]+=B, a[i+]--;
        a[i] -= b[i];
    }
    int p=b.n+;
    while(a[p]<) a[p]+=B, a[++p]--;
    while(a.n> && a[a.n]==) a.n--;
    return a;
}

void Print(Big &a) {
    printf("%d", a[a.n]);
    for(int i=a.n-; i>=; i--) printf("%04d", a[i]);
}
int main(){
    scanf("%d",&n);
    ans[].ini();ans[].ini();
    for(int i=;i<=n;i++){
        Big a;
        a=ans[i-]+ans[i-];
        ans[i]=a*(i-);
    }
    Print(ans[n]);
    return ;
}