【Uva 1543】Telescope
2017-10-04 18:44
309 查看
【Link】:
【Description】
给你一个圆和圆周上的n(3≤n≤40)个不同点。请选择其中的m(3≤m≤n)个,按照在圆 周上的顺序连成一个m边形,使得它的面积最大。
【Solution】
DP;
设f[i][j][k]表示在第i到第j个点之间一定选择了i和j的条件下选择了k个点组成的多边形的最大面积;
这样表示状态;在增加一个点的时候;
只要算i..j和新加的那个点组成的三角形面积,然后把这个三角形的面积加上去就好了;
逆推的话;
f[i][j][k] = f[i][l][k-1] + area(i,l,j)
l∈[i..j-1];
当k<3的时候,f[i][j][k]=0;
其中area(i,l,j)表示点i,l,j组成的三角形的面积;
面积可以用海伦公式搞;
其中两点之间的弦长为2∗sin(a2)
则用题中给的东西表示就为
2∗sin(|p[i]−p[j]|∗π)
这样就能求出三角形的三边了;
【NumberOf WA】
0
【Reviw】
状态表示得清晰一点,DP就不难写了;
【Code】
#include <bits/stdc++.h> using namespace std; #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define LL long long #define rep1(i,a,b) for (int i = a;i <= b;i++) #define rep2(i,a,b) for (int i = a;i >= b;i--) #define mp make_pair #define pb push_back #define fi first #define se second #define ms(x,y) memset(x,y,sizeof x) #define ri(x) scanf("%d",&x) #define rl(x) scanf("%lld",&x) #define rs(x) scanf("%s",x+1) #define rd(x) scanf("%lf",&x) #define oi(x) printf("%d",x) #define ol(x) printf("%lld",x) #define oc putchar(' ') #define os(x) printf(x) #define all(x) x.begin(),x.end() #define Open() freopen("F:\\rush.txt","r",stdin) #define Close() ios::sync_with_stdio(0) typedef pair<int,int> pii; typedef pair<LL,LL> pll; const int dx[9] = {0,1,-1,0,0,-1,-1,1,1}; const int dy[9] = {0,0,0,-1,1,-1,1,-1,1}; const double pi = acos(-1.0); const int N = 40; const int INF = 0x3f3f3f3f; int n,m; double p[N+10],f[N+10][N+10][N+10]; double dis(int i,int j){ double tp = abs(p[j]-p[i]); if (tp>0.5) tp = 1 - tp; return 2*sin(tp*pi); } double area(int i,int j,int k){ double a = dis(i,j),b = dis(i,k),c = dis(j,k); double temp = (a+b+c)/2.0; return sqrt(temp*(temp-a)*(temp-b)*(temp-c)); } double dfs(int l,int r,int num){ if (f[l][r][num]>=0) return f[l][r][num]; if (num < 3) return 0; double temp = 0; rep1(i,l,r-1) temp = max(temp,dfs(l,i,num-1) + area(l,i,r)); return f[l][r][num] = temp; } int main(){ //Open(); //Close(); while (~ri(n)){ ri(m); if (n == 0 && m == 0) break; rep1(i,1,n) rd(p[i]); rep1(i,1,n) rep1(j,1,n) rep1(k,1,n) f[i][j][k] = -1; double ans = 0; rep1(i,1,n) rep1(j,i+1,n) ans = max(ans,dfs(i,j,m)); printf("%.6f\n",ans); } return 0; }
相关文章推荐
- UVA 1543 Telescope
- 【UVA1543】Telescope
- uva 1543 - Telescope(dp+几何)
- Telescope UVA - 1543
- uva 1543 Telescope
- UVA 1543 - Telescope(dp+几何)
- UVA - 1543 Telescope——最优三角剖分
- UVA - 1543(dp)
- uva 1543 dp
- UVA - 1543 Telescope dp+几何
- uva1543
- UVA - 10023 - Square root (模拟手算开方)
- uva 1395 - Slim Span(生成树)
- UVA - 10591 Happy Number(hash)
- uva 10720 Graph Construction
- Uva10082 WERTYU -S.B.S.
- UVa 1585 Score 得分 (java版本)
- UVA242 背包问题&紫书dp习题9-5
- UVa 11021 麻球繁衍
- 例题6-15 给任务排序(Ordering Tasks, UVa 10305)