【FJOI2007】轮状病毒，以及一类普通的生成树计数

原来只接触过盾盾油头之类的、最小生成树计数。并且同边权的边数很小。

然后看到FJOI2007轮状病毒和JSOI2010巨额奖金，都可以用一个神奇的东西——基尔霍夫矩阵。

介绍可以看周冬的《生成树的计数及其应用》

于是我们可以得到一个灰常神奇的东西，这个东西可以n^3的计算无向图的生成树个数。其中n为图的点数。

此外这道题还有一个Trick：

因为算n-1阶主子式的时候我们显然是把中间那个点对应的行和列删掉。

那么得到行列式：

-1 3 -1 0 0 0 0 0 .. 0 0

0 -1 3 -1 0 0 0 0 .. 0 0

…………………………………

0 0 0 0 0 0 0 .. -1 3 -1

-1 0 0 0 0 0 0 0 .. -1 3

3 -1 0 0 0 0 0 0 0 .. -1

算出来之后发现了一个神奇的事情：f[n] = 3f[n - 1] - f[n - 2] + 2、

虽不知怎么证明但觉厉

代码如下：（此题的找规律解法：f[n]=C(n,n)*1+C(n,n-1)*2+...+C(n,2)*n-1+C(n,1)*n，很神奇是不。）

（另外此题是有DP解法的，也准备写一下，但是还是按照惯例开坑了...对不起）

#include<iostream>
#include<cstdio> 
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<climits>
#include<map>
#include<set>
#define ot "%d"
#define kg " "
#define kh "\n"
#ifdef WIN32
#define otl "%I64d"
#else
#define otl "%lld"
#endif
#define max(a, b) ({int _ = (a), __ = (b); _ > __ ? _ : __;})
#define min(a, b) ({int _ = (a), __ = (b); _ < __ ? _ : __;})
#define swap(a, b) ({int _ = (a); (a) = (b); (b) = _;})
#define maxn 105
#define fort(a, b, c) for (int a = (b); a <= (c); ++a)
#define ford(a, b, c) for (int a = (b); a >= (c); --a)
#define ms(a, b) memset(a, b, sizeof a)
#define termin() {puts("No solution");exit(0);}
#define tp typedef
#define mo 100000000
#define he 10000000
#define cps CLOCKS_PER_SEC
#define ct clock_t

using namespace std;

tp int gjd[505];

int n;
gjd f[maxn];

void init()
{
	freopen("input.txt", "r", stdin);
	freopen("output.txt", "w", stdout);
}

int main()
{
	init();
	scanf(ot, &n);
	if (n == 1) puts("1"); else if (n == 2) puts("5");
	else
	{
		int t;
		memset(f, 0, sizeof f);
		f[1][0] = f[1][1] = 1;
		f[2][0] = 1, f[2][1] = 5;
		fort(i, 3, n)
		{
			f[i][1] = 2;		
			fort(j, 1, f[i - 1][0])
			{
				f[i][j] = f[i - 1][j] * 3 - f[i - 2][j] + f[i][j];
				while (f[i][j] >= mo) f[i][j] -= mo, ++f[i][j + 1];
				if (f[i][j] < 0) --f[i][j + 1], f[i][j] += mo;
			}
			if (f[i][(f[i][0] = f[i - 1][0]) + 1]) ++f[i][0];
		}
		printf(ot, f[n][f[n][0]]);
		ford(i, f[n][0] - 1, 1)
		{
			t = f[n][i];
			while (f[n][i] < he) printf("0"), f[n][i] *= 10;
			printf(ot, t);
		}
	}
	return 0;
}

ps:模板越来越长...求破..说到底还是我太懒了..什么ot olt tp cps ms..