1. 程式人生 > >最優三角劃分(動態規劃) By ACReaper

最優三角劃分(動態規劃) By ACReaper

1.有狀態轉移方程為f(i,j) = max(f(i,k) + f(k,j) + w(i,j,k));f(i,j)表示從vi到vj的簡單通路(我們預設vj指向vi),這樣才能構成凸多邊形,同時表示在這個凸多邊形中權和的最大值。

    2.由3邊的凸多邊形也就是三角形開始遞推上去,構建整個動態規劃。
#include <stdio.h>
#define MAXN 1000
int f[MAXN][MAXN];
int w[MAXN][1];
int n;
int ww(int i,int j,int k){
    return w[i][0] + w[j][0] + w[k][0];
}

int main(){
    while(scanf("%d",&n) != EOF){
        
        for(int i =1 ; i <= n; i++){//輸入各個頂點權值 
            scanf("%d",&w[i][0]);
        }
        for(int i = 1; i <= n; i++){//無法構成凸邊形,說白點就是邊為1條 
            f[i][i + 1] = 0;
        } 
        //從通路長度為2的開始遞推上去,通路為1即為一條邊,通路為2為三角形 
        for(int d = 2; d <= n - 1; d++){//因為最後一個即為n邊形,只需要n - 1唱的的通路,加上最後一個頂點指向第一個構成迴路,也就是凸多變形,tips:d = 2是三角形
            for(int i = 1; i <= n - d; i++){//列舉起點 
                int j = d + i;//終點 
                f[i][j] = f[i + 1][j] + ww(i,i + 1, j); 
                for(int k = i + 1; k < j; k++){
                    int temp = f[i][k] + f[k][j] + ww(i,j,k);
                    f[i][j] = f[i][j] < temp?f[i][j] : temp;
                }
            }
            
        } 
        printf("%d\n",f[1][n]);
    }
    return 0;
}