【计算几何】【分类讨论】Gym - 101173C - Convex Contour
注意等边三角形的上顶点是卡不到边界上的。
于是整个凸包分成三部分:左边的连续的三角形、中间的、右边的连续的三角形。
套个计算几何板子求个三角形顶点到圆的切线、三角形顶点到正方形左上角距离啥的就行了,分类比较多。
#include<cstdio>
#include<cmath>
using namespace std;
const double PI=acos(-1.0);
int n;
char a[25];
struct Point{
double x,y;
double length(){
return sqrt(x*x+y*y);
}
};
typedef Point Vector;
Vector operator - (const Point &a,const Point &b){
return (Vector){a.x-b.x,a.y-b.y};
}
Vector operator + (const Vector &a,const Vector &b){
return (Vector){a.x+b.x,a.y+b.y};
}
Vector operator * (const double &K,const Vector &v){
return (Vector){K*v.x,K*v.y};
}
Vector Rotate(Vector A,double rad)
{
return (Vector){A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad)};
}
Vector unit(Vector A){
double l=A.length();
return (Vector){A.x/l,A.y/l};
}
double sqr(double x){
return x*x;
}
int main(){
scanf("%d%s",&n,a+1);
bool alltr=1;
for(int i=1;i<=n;++i){
if(a[i]!='T'){
alltr=0;
break;
}
}
if(alltr){
printf("%d\n",2*n+1);
}
else{
double ans=0;
int I,J;
for(int i=1;i<=n;++i){
if(a[i]!='T'){
if(i!=1){
if(a[i]=='S'){
ans+=((Point){0.5,sqrt(3.0)/2.0}-(Point){(double)(i-1),1.0}).length();
}
else if(a[i]=='C'){
Vector v=(Point){(double)i-0.5,0.5}-(Point){0.5,sqrt(3.0)/2.0};
ans+=sqrt(sqr(v.length())-0.5*0.5);
v=Rotate(v,asin(0.5/v.length()));
Point p=(Point){0.5,sqrt(3.0)/2.0}+ans*unit(v);
double xian=((Vector){(double)i-0.5,1.0}-p).length();
double jiao=acos((sqr(xian)-0.5*0.5*2.0)/(-2.0*0.5*0.5));
ans+=jiao*0.5;
}
}
else{
if(a[i]=='S'){
ans+=2.0;
}
else if(a[i]=='C'){
ans+=PI*0.5;
}
}
I=i;
break;
}
}
for(int j=n,i=1;j>=1;++i,--j){
if(a[j]!='T'){
if(j!=n){
if(a[j]=='S'){
ans+=((Point){0.5,sqrt(3.0)/2.0}-(Point){(double)(i-1),1.0}).length();
}
else if(a[j]=='C'){
Vector v=(Point){(double)i-0.5,0.5}-(Point){0.5,sqrt(3.0)/2.0};
double d=sqrt(sqr(v.length())-0.5*0.5);
ans+=d;
v=Rotate(v,asin(0.5/v.length()));
Point p=(Point){0.5,sqrt(3.0)/2.0}+d*unit(v);
double xian=((Vector){(double)i-0.5,1.0}-p).length();
double jiao=acos((sqr(xian)-0.5*0.5*2.0)/(-2.0*0.5*0.5));
ans+=jiao*0.5;
}
}
else{
if(a[j]=='S'){
ans+=2.0;
}
else if(a[j]=='C'){
ans+=PI*0.5;
}
}
J=j;
break;
}
}
if(I!=1 && J!=n){
ans+=((double)(I-1)+1.0);
ans+=((double)(n-J)+1.0);
ans+=(double)(J-I+1);
if(a[I]=='S' && a[J]=='S'){
ans+=(double)(J-I+1);
}
else if(a[I]=='C' && a[J]=='C'){
ans+=(double)(J-I);
}
else{
ans+=((double)(J-I)+0.5);
}
}
else if(I==1 && J==n){
ans+=(double)(n-1)*2.0;
}
else if(I==1 && J!=n){
ans+=((double)(n-J)+1.0);
if(a[J]=='S'){
ans+=((double)J-0.5)*2.0;
}
else{
ans+=(((double)J-0.5)*2.0-0.5);
} }
else{
ans+=((double)(I-1)+1.0);
if(a[I]=='S'){
ans+=((double)(n-I+1)-0.5)*2.0;
}
else{
ans+=(((double)(n-I+1)-0.5)*2.0-0.5);
}
}
printf("%.11f\n",ans);
}
return 0;
}
【计算几何】【分类讨论】Gym - 101173C - Convex Contour的更多相关文章
- dp+分类讨论 Gym 101128E
题目链接:http://codeforces.com/gym/101128 感觉这个人写的不错的(我只看了题目大意):http://blog.csdn.net/v5zsq/article/detail ...
- Bzoj4558:分类讨论 计算几何 组合数学
国际惯例的题面: 这题让我爆肝啦......这种计数显然容斥,正好不含任何坏点的我们不会算,但是我们能算至少含零个坏点的,至少含一个坏点的,至少含两个坏点的......所以最终的答案就是(至少含零个坏 ...
- Codeforces 460D Little Victor and Set --分类讨论+构造
题意:从区间[L,R]中选取不多于k个数,使这些数异或和尽量小,输出最小异或和以及选取的那些数. 解法:分类讨论. 设选取k个数. 1. k=4的时候如果区间长度>=4且L是偶数,那么可以构造四 ...
- BZOJ-1067 降雨量 线段树+分类讨论
这道B题,刚的不行,各种碎点及其容易忽略,受不鸟了直接 1067: [SCOI2007]降雨量 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 2859 ...
- UVaLive 6862 Triples (数学+分类讨论)
题意:给定一个n和m,问你x^j + y^j = z^j 的数量有多少个,其中0 <= x <= y <= z <= m, j = 2, 3, 4, ... n. 析:是一个数 ...
- 枚举(分类讨论):BZOJ 1177: [Apio2009]Oil
1177: [Apio2009]Oil Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 1477 Solved: 589[Submit] Descri ...
- Educational Codeforces Round 63 (Rated for Div. 2) D. Beautiful Array 分类讨论连续递推dp
题意:给出一个 数列 和一个x 可以对数列一个连续的部分 每个数乘以x 问该序列可以达到的最大连续序列和是多少 思路: 不是所有区间题目都是线段树!!!!!! 这题其实是一个很简单的dp 使用的是分 ...
- 【cf789B】Masha and geometric depression(分类讨论/暴力)
B. Masha and geometric depression 题意 在黑板上写数列,首项是b,公比是q,超过l时就停止不写.给定m个数,遇到后跳过不写.问一共写多少个数,如果无穷个输出inf. ...
- P2331 [SCOI2005]最大子矩阵 (动规:分类讨论状态)
题目链接:传送门 题目: 题目描述 这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩阵分值之和最大.注意:选出的k个子矩阵不能相互重叠. 输入输出格式 输入格式: 第一行为n,m,k( ...
随机推荐
- 如何免费上传4G以上大文件至百度云网盘
百度云网盘的容量高达2048G,因而如今使用百度云网盘的用户也越来越多, 但是百度云中如果要上传超过4G的大文件,必须要升级VIP才行,但这需要收费.那么,超过4G以上的大文件我们该怎样上传到百度云呢 ...
- javascript执行上下文学习一
原文: http://web.jobbole.com/84044/ http://blog.csdn.net/github_34514750/article/details/52901781 1.三种 ...
- arch中yaourt的安装和使用
yaourt-Yet AnOther User Repository Tool Yaourt是archlinux方便使用的关键部件之一,但没有被整合到系统安装中的工具.建议在装完系统重启之后,更新完p ...
- Unix/Linux Command Reference
- Perl6 Bailador框架(2):路径设置
use v6; use Bailador; =begin pod get表示是get发送 post表示是post发送 get/post 后面的 '/name' 表示是路径 => sub {} 是 ...
- Perl6 Bailador框架(1):开始
use v6; use Bailador; get '/' => sub { '<h1><center>Hello, World</center></h ...
- Python-字符串处理 str.format()
Python中内置的%操作符可用于格式化字符串操作,控制字符串的呈现格式.Python中还有其他的格式化字符串的方式,但%操作符的使用是最方便的. 另外python还有一个更强大的字符串处理函数 st ...
- 使用PTGui软件将全景图变成鱼眼图
把全景图变成鱼眼图.方法一部分是自己研究的,一部分是参考学妹街景合成鱼眼照片的方法. 需要使用的软件是PTGui.是个收费软件,价格还不便宜.操作一下,安装完后就可以开始合成鱼眼图了. 加载图像 打开 ...
- php 读写 csv文件
读取csv function input_csv($handle) { $out = array (); $n = 0; while ($data = fgetcsv($handle, 10000)) ...
- Shell三剑客之sed命令
Sed简介 Sed是Stream Editor(流编辑器)缩写,是操作.过滤和转换文本内容的强大工具,常用功能有增删改查. Sed命令执行流程 Sed语法格式 Sed [option] ‘[匹配][处 ...