【计算几何】【分类讨论】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( ...
随机推荐
- Kill windows和linux 进程
Windows
- SpringBoot工程目录配置
Spring Boot建议的目录结果如下: root package结构:com.example.myproject com +- example +- myproject +- Applicat ...
- Java案例之随机验证码功能实现
实现的功能比较简单,就是随机产生了四个字符然后输出.效果图如下,下面我会详细说一下实现这个功能用到了那些知识点,并且会把 这些知识点详细的介绍出来.哈哈 ,大神勿喷,对于初学Java的人帮助应该蛮大的 ...
- perl HTML::LinkExtor模块(2)
use LWP::Simple; use HTML::LinkExtor; $html_code = get("https://tieba.baidu.com/p/4929234512&qu ...
- 9.quartus_warning_altera_reserved_tck
编译的时候没有注意,整个工程都可以在板子上跑起来.但是做Powerplay的时候,出现了这个Critical Warning:. Critical Warning: The following clo ...
- python科学计算整理
网站: http://bokeh.pydata.org/gallery.html
- java===java基础学习(8)---静态域与静态方法
静态域:如果将域定义为static,每个类中只有一个这样的域.而每一个对象对于所有的实例域却都有自己的一份拷贝.例如,加入需要给每一个雇员赋予唯一的标识码.这里给的Employee类添加一个实例域id ...
- python实战===用python对比两张图片的不同
from PIL import Image from PIL import ImageChops def compare_images(path_one, path_two, diff_save_lo ...
- App推广
一行业常见名词解释 开发商:也叫CP,即Content Provider内容提供商的英文首字母缩写. 发行商(运营商):即代理CP开发出来的产品. 渠道:拥有用户,能够进行流量分发的公司,即可成为渠道 ...
- Leetcode 之Length of Last Word(37)
扫描每个WORD的长度并记录即可. int lengthOfLast(const char *s) { //扫描统计每个word的长度 ; while (*s) { if (*s++ != ' ')/ ...