题意：给出平面上\(n\)个点，要求选出\(k\)个点，使得这些点形成一个凸包，且凸包内部没有点，求最大面积。无解输出\(0\)。

题解：枚举凸包最左的点\(p\)，删除所有在\(p\)左边的点，然后把\(p\)定为原点。将所有点按极角排序，相邻两个点之间连边，那么会形成一个星状多边形，合法的凸包一定在这个多边形内部。

先考虑求出这张图的visibility graph，显然合法的凸包所有边都是visibility graph上的边。求法大概是逆时针枚举所有点，对于每个点维护一个队列维护未来可能加入的边，实际上是对于每个点\(i\)，维护所有满足\(ij\)在visibility graph上，并且当前还没找到\(k(k>i)\)使得\(jk\)在visibility graph上的\(j\)，详见代码。复杂度\(O(E)\)。

考虑在visibility graph上DP。顺时针枚举所有点，设\(f_{i,j,k}\)表示最后一条选取的边为\(i,j\)，选了\(k\)条边的最大面积。转移时可以枚举一个\(l\)，如果\(i,j\)和\(l,i\)这两条边可以同时存在（不会使得凸包不满足凸性）则可以转移到\(f_{l,i,k+1}\)。朴素DP复杂度\(O(n^3k)\)，可以对于每个点将转移出去和进来的边分别排序后（其实根据求visibility graph的过程，这些边是已经排好序的）双指针+前缀和优化，复杂度\(O(n^2k)\)。

总复杂度\(O(n^3k)\)。

#include<bits/stdc++.h>
using namespace std;
const int N = 210;
typedef long long ll;
typedef double db;
#define pb push_backint gi() {int x = 0, o = 1;char ch = getchar();while((ch < '0' || ch > '9') && ch != '-') {ch = getchar();}if(ch == '-') {o = -1, ch = getchar();}while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0', ch = getchar();}return x * o;
}struct point {int x, y;db k;point(int x = 0, int y = 0): x(x), y(y) {k = atan2(y, x);}point operator-(const point &A) const {return point(x - A.x, y - A.y);}ll operator%(const point &A) const {return 1ll * x * A.y - 1ll * y * A.x;}bool operator<(const point &A) const {return k < A.k;}
} a[N], p[N];int n, m, tt;
ll f[N][N][55], mx[55], ans = 0;
queue<int> q[N];
vector<int> E[N], G[N];void add(int x, int y) {while(!q[x].empty() && (p[q[x].front()] - p[x]) % (p[y] - p[x]) < 0) {add(q[x].front(), y), q[x].pop();}G[x].pb(y), E[y].pb(x), q[y].push(x);
}int main() {
#ifndef ONLINE_JUDGEfreopen("a.in", "r", stdin);freopen("a.out", "w", stdout);
#endifcin >> n >> m;for(int i = 1; i <= n; i++) {a[i].x = gi(), a[i].y = gi();}for(int s = 1; s <= n; s++) {tt = 0;for(int i = 1; i <= n; i++) if(a[i].x > a[s].x || (a[i].x == a[s].x && a[i].y > a[s].y)) {p[++tt] = a[i] - a[s];}sort(p + 1, p + tt + 1);for(int i = 1; i <= tt; i++) {E[i].clear(), G[i].clear();while(!q[i].empty()) {q[i].pop();}}for(int i = 1; i < tt; i++) {add(i, i + 1);}memset(f, 0xc0, sizeof(f));for(int i = tt; i; i--) {memset(mx, 0xc0, sizeof(mx));reverse(E[i].begin(), E[i].end());int cur = G[i].size() - 1;for(auto j : E[i]) {f[i][j][1] = p[j] % p[i];while(~cur && (p[j] - p[i]) % (p[G[i][cur]] - p[i]) < 0) {for(int k = 1; k < m; k++) {mx[k] = max(mx[k], f[G[i][cur]][i][k]);}--cur;}for(int k = 1; k < m; k++) {f[i][j][k + 1] = mx[k] + p[j] % p[i];}}}for(int i = 1; i <= tt; i++)for(auto j : E[i]) {ans = max(ans, f[i][j][m - 2]);}}printf("%.2lf\n", 1.0 * ans / 2);return 0;
}