POJ 3597 Polygon Division 多边形剖分

题目链接：

http://poj.org/problem?id=3597

Polygon Division

Time Limit: 2000MS
Memory Limit: 131072K #### 问题描述 > Given a regular polygon, there are numerous ways to divide it into several triangles and/or quadrangles by adding some diagonals that do not properly intersect each other. For example, Figure 4 shows all ten different divisions of a regular pentagon into triangles and quadrangles. > > Figure 4: Divisions of a regular pentagon into triangles and quadrangles > > Given n, the number of sides of the polygon, compute the number of such divisions. #### 输入 > The input contains multiple test cases. Each test case consists of a single integer n (3 ≤ n ≤ 5000) on a separate line. The input ends where EOF is met. #### 输出 > For each test case, print the answer modulo 264 on a separate line. ####样例输入 > 3 > 4 > 5 > 6 > 7 > 8 > 9 > 10

样例输出

1
3
10
38
154
654
2871
12925

题意

把正凸多边形剖成四边形和三角形的所有方案。

题解

[port]

代码

include

include

include

include

include

include

include

include

include

include

include

include

include

include

include

using namespace std;
#define X first
#define Y second
#define mkp make_pair
#define lson (o<<1)
#define rson ((o<<1)|1)
#define mid (l+(r-l)/2)
#define sz() size()
#define pb(v) push_back(v)
#define all(o) (o).begin(),(o).end()
#define clr(a,v) memset(a,v,sizeof(a))
#define bug(a) cout<<#a<<" = "<<a<<endl
#define rep(i,a,b) for(int i=a;i<(b);i++)
#define scf scanf
#define prf printf

typedef long long LL;
typedef vector VI;
typedef pair<int,int> PII;
typedef vector<pair<int,int> > VPII;
typedef unsigned long long ULL;

const int INF=0x3f3f3f3f;
const LL INFL=0x3f3f3f3f3f3f3f3fLL;
const double eps=1e-8;
const double PI = acos(-1.0);

//start----------------------------------------------------------------------

const int maxn=5555;

ULL f[maxn],g[maxn];

void pre(){
clr(g,0),clr(f,0);
g[1]=g[2]=f[1]=f[2]=1;
for(int i=3;i<maxn;i++){
g[i]=0;
for(int j=2;j<i;j++) g[i]+=f[j]f[i-j+1];
f[i]=g[i];
for(int j=2;j<=i-2;j++) f[i]+=f[j]g[i-j+1];
}
}

int main() {
pre();
int tc,kase=0;
int n;

while(scf("%d",&n)==1){
if(n<3) puts("0");
prf("%llu\n",f[n]);
}
return 0;
}

//end-----------------------------------------------------------------------