POJ出栈序列统计

2:出栈序列统计

总时间限制: 

100ms  

内存限制: 

64kB

描述

栈是常用的一种数据结构，有n个元素在栈顶端一侧等待进栈，栈顶端另一侧是出栈序列。你已经知道栈的操作有两种：push和pop，前者是将一个元素进栈，后者是将栈顶元素弹出。现在要使用这两种操作，由一个操作序列可以得到一系列的输出序列。请你编程求出对于给定的n，计算并输出由操作数序列1，2，…，n，经过一系列操作可能得到的输出序列总数。

输入

就一个数n(1≤n≤15)。

输出

一个数，即可能输出序列的总数目。

样例输入

3

样例输出

5

提示

先了解栈的两种基本操作，进栈push就是将元素放入栈顶，栈顶指针上移一位，等待进栈队列也上移一位，出栈pop是将栈顶元素弹出，同时栈顶指针下移一位。
　　 用一个过程采模拟进出栈的过程，可以通过循环加递归来实现回溯：重复这样的过程，如果可以进栈则进一个元素，如果可以出栈则出一个元素。就这样一个一个地试探下去，当出栈元素个数达到n时就计数一次(这也是递归调用结束的条件)。

import java.math.BigDecimal;
import java.util.Scanner;

public class Main {
	public static Long jiecheng(int n) {
		if (n < 1)
			return 1L;
		return n * jiecheng(n - 1);
	}

	public static BigDecimal factorial(BigDecimal n) {
		BigDecimal bd1 = new BigDecimal(1);
		BigDecimal bd2 = new BigDecimal(2);
		BigDecimal result = bd1;
		while (n.compareTo(bd1) > 0) {
			result = result.multiply(n.multiply(n.subtract(bd1)));
			n = n.subtract(bd2);
		}
		return result;
	}

	public static void main(String[] args) {
		Scanner sc = new Scanner(System.in);
		BigDecimal n = sc.nextBigDecimal();
		BigDecimal sum = new BigDecimal(0);
		sum = factorial(n.multiply(new BigDecimal(2))).divide(
				(factorial(n.add(new BigDecimal(1)))).multiply(factorial(n)));
		 System.out.println(sum);
	}
}

 看数据结构，看到栈的时候，发现1个元素进栈，有1种出栈顺序；2个元素进栈，有2种出栈顺序；3个元素进栈，有5种出栈顺序，那么一个很自然地问题就是n个元素进栈，共有多少种出栈顺序？说来惭愧，以前学数据结构的时候竟然没有考虑过这个问题。最近在看动态规划，所以“子问题”这3个字一直在我脑中徘徊，于是解决这个问题的时候我也是用类似“子问题”的方法，说白了就是递推公式。

我们把n个元素的出栈个数的记为f(n), 那么对于1,2,3, 我们很容易得出：                    f(1) = 1     //即 1

                         f(2) = 2     //即 12、21

                         f(3) = 5     //即 123、132、213、321、231

然后我们来考虑f(4), 我们给4个元素编号为a,b,c,d, 那么考虑：元素a只可能出现在1号位置，2号位置，3号位置和4号位置(很容易理解，一共就4个位置，比如abcd,元素a就在1号位置)。

分析：

 1) 如果元素a在1号位置，那么只可能a进栈，马上出栈，此时还剩元素b、c、d等待操作，就是子问题f(3)；

 2) 如果元素a在2号位置，那么一定有一个元素比a先出栈，即有f(1)种可能顺序（只能是b），还剩c、d，即f(2)，根据乘法原理，一共的顺序个数为f(1) * f(2)；

 3) 如果元素a在3号位置，那么一定有两个元素比1先出栈，即有f(2)种可能顺序（只能是b、c），还剩d，即f(1)，根据乘法原理，一共的顺序个数为f(2) * f(1)；

 4) 如果元素a在4号位置，那么一定是a先进栈，最后出栈，那么元素b、c、d的出栈顺序即是此小问题的解，即f(3)；

 

结合所有情况，即f(4) = f(3) + f(2) * f(1) + f(1) * f(2) + f(3);

为了规整化，我们定义f(0) = 1；于是f(4)可以重新写为：

f(4) = f(0)*f(3) + f(1)*f(2) + f(2) * f(1) + f(3)*f(0)

 

然后我们推广到n，推广思路和n=4时完全一样，于是我们可以得到：

f(n) = f(0)*f(n-1) + f(1)*f(n-2) + ... + f(n-1)*f(0)

 

即

                  

 

但这只是一个递推公式（若编程实现，需要维护一个一维数组，时间复杂度为O(n^2)）。怎么把它转化为通项公式呢，复杂度仅为O(1)?

于是上网搜索一下，原来真的有这么一个公式：C(2n,n)/(n+1) （C(2n,n)表示2n里取n），并且有个名字叫Catalan数。附上wiki的链接，写得太详细了：http://en.wikipedia.org/wiki/Catalan_number