UVA 1543 Telescope

原创 2016年08月29日 12:27:35

题目链接:http://acm.hust.edu.cn/vjudge/problem/41494


题意:一个圆上有n个点,选择其中的m个点按照顺序连成m边形,求最大的面积。


思路:dp[i][j][k],表示从i号点开始作为起点的j个点中选择k个点连成图形的最大面积,且区间内的左右边界点必须选择。

dp[i][j][k] = max( dp[i][l][k-1] + area(i,l,j) ) 枚举起点,长度,和选择点的个数k,以及第k-1个点选择的位置进行转移。

还有就是给定三个点,计算他们围成的面积,可以根据角度计算出各边长,然后用海伦公式计算面积。


#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod 100000007

const int maxn = 50;
const double pi = acos(-1.0);
double x[maxn];
double dp[maxn][maxn][maxn];
double area[maxn][maxn][maxn];

double cal( int a,int b,int c )
{
    double a1 = 2 * sin( min( fabs(x[a]-x[b]) , 1 - fabs(x[a]-x[b]) )*pi );
    double a2 = 2 * sin( min( fabs(x[a]-x[c]) , 1 - fabs(x[a]-x[c]) )*pi );
    double a3 = 2 * sin( min( fabs(x[c]-x[b]) , 1 - fabs(x[c]-x[b]) )*pi );
    double p = ( a1 + a2 + a3 ) / 2;
    return sqrt( p * (p-a1) * (p-a2) * (p-a3) );
}
int n,m;

void init()
{
    rep(i,1,n)
        rep(j,1,n)
            rep(k,1,n)
            if ( i != j && i != k && k != j )
            area[i][j][k] = area[i][k][j] = area[j][i][k] = area[j][k][i] = area[k][i][j] = area[k][j][i] = cal( i,j,k );
}

int main()
{
    while( cin>>n>>m )
    {
        if ( n + m == 0 ) break;
        rep(i,1,n) scanf("%lf",&x[i]);
        Clean(dp,0);
        init();
        rep(l,3,n) //长度
            rep(i,1,n) //起点
            {
                int uplim = min( m , l ); //此区间最多选择的点数
                rep( k , 3 , uplim )
                    Rrep(mid,l-1,2) //枚举第k-1个点的位置
                    {
                        if ( mid < k-1 ) break; //此区间长度不够k-1,放不下k-1个点
                        dp[i][l][k] = max( dp[i][l][k] , dp[i][mid][k-1] + area[i][ (i+mid-2)%n + 1 ][ (i+l-2)%n + 1 ] ) ;
                    }
            }
        double ans = 0;
        rep(i,1,n)
            rep(j,m,n)
            ans = max( ans , dp[i][j][m] );
        printf("%0.6f\n",ans);
    }
    return 0;
}





版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

uva 1543 Telescope

终于找到一点自信,没有看别人题解ac这道lrj紫书上的题目,一旦想到用海伦公式其他的就很简单了,我发现用递归写dp,真的比循环写dp要容易些。主要思路就是dfs函数中的循环语句。其实这道题目也可以看做...

uva1543

dp[i][j][k] 从第i个点到第j个点之间取k个点去的面积的最大值 #include #include #include #include using namespace std; const...

POJ 3675 Telescope

多边形与圆的相交面积问题 TelescopeTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 1012 Accepted: 266De...

POJ 3675 Telescope

POJ 3675 Telescope 题目链接:POJ 3675 题意: 圆和多边形的公共面积。 分析: 将多边形分成n个三角形和圆的公共面积。 实现1: /* ...

poj 3675 Telescope(三角剖分求多边形和圆交面积)

Telescope Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 2271   Accepted: 673 ...

PKU 3675 Telescope

Telescope Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 1178   Acce...

Codeforces Gym 100800L Telescope (二维树状数组 或 暴力水过)

题目大意: 貌似这个题就是题目意思不太好懂...其实题目本身还是比较简单的 题意是给出一张照片(最多1000*1000的十六进制矩阵) 然后这是一个星体的照片, 并且告诉的是每个位置(x, y)的值是...

HDU 1543:Avoid The Lakes 【深搜】

Avoid The Lakes Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 131072/65536K (Java/Other) ...

刷题打卡1 poj1543 ,poj1006,poj1013

刷题打卡1 poj1543 ,poj1006,poj1013

POJ 1543 Perfect Cubes(我的水题之路——四重暴力水)

Perfect Cubes Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 10707   ...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)