1. 程式人生 > >動態規劃-凸多邊形

動態規劃-凸多邊形

時間限制:10000ms 單點時限:1000ms 記憶體限制:256MB

描述

給定一個凸多邊形的N個頂點。你需要在凸多邊形內找到M個點,使得這M個點也圍成一個凸多邊形,並且圍成的面積儘可能大。

輸入

第一行包含兩個整數N和M,意義如前文所述。

接下來N行,每行兩個整數Ai和Bi,表示按照逆時針順序排列的凸多邊形頂點座標。

對於30%的資料,滿足N<=5

對於100%的資料,滿足N<=100

對於100%的資料,滿足3<=M<N, |Ai|,|Bi|<=10000

輸出

輸出新凸多邊形最大的面積,保留兩位小數。

樣例輸入
4 3
0 0
1 0
1 1
0 1
樣例輸出
0.50
import java.util.*;

public class Main {
    public static void main(String[] args) {

        Scanner in = new Scanner(System.in);

        while (in.hasNext()) {
            int n = in.nextInt(),m = in.nextInt();
            int[] x = new int[n];
            int[] y = new int[n];
            for (int i = 0; i < n; i++) {
                x[i] = in.nextInt();
                y[i] = in.nextInt();
            }

            double[][][] dp = new double[n][n][n];
            double[][][] s = new double[n][n][n];
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                    for(int k=0;k<n;k++) {
                        dp[i][j][k] = 0;
                        s[i][j][k]=0.5*Math.abs(x[i]*y[j]+x[j]*y[k]+x[k]*y[i]-x[i]*y[k]-x[j]*y[i]-x[k]*y[j]);
                    }

            for(int k=3;k<=m;k++)
                for(int i=0;i<n;i++)
                    for(int j=(i+1)%n;j!=i;j=(j+1)%n)
                        for(int u=(i+1)%n;u!=j;u=(u+1)%n)
                            dp[i][j][k]=Math.max(dp[i][j][k],dp[i][u][k-1]+s[i][u][j]);


            double max=0;
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                    max=Math.max(max, dp[i][j][m]);
            System.out.println(String.format("%.2f",max));

        }

    }
}