abstract:

  1. V const & a 加速

F. Mirror

题意

链接

问题:

有n个人在y=0的平面上(及xoz平面)。z=0平面上有一面镜子(边平行于坐标轴)。z=a平面上有q个点(保证a大于所有人的z坐标)。 所有人面朝镜子,且在镜子和q个点之间(即每个人的z坐标保证0<z<a)。

问对于某个点,让所有人能够通过镜子看到那个点的镜子的最小面积。


题解

三维几何+镜像+凸包

首先考虑镜面,我们可以通过(初中科学的)镜面反射原理,关于z=0做出z=a的对称平面z=-a。问题就变成了n个人看z=-a上的某个点。(下图绿点是人,红点是询问点)

然后观察,镜子的高和宽是独立的。 于是我们分别求它们的最大值即可。

求高比较简单,我们朝-x方向看yoz平面。通过把每个人跟点Q的像Q‘相连,我们发现离Q’z轴距离最近的人对应着镜子的下边界,最远的人对应着上边界,通过维护所有人z坐标的max_z&min_z以及相似三角形可以直接求出两个边界,复杂度为O(1)。

我们用同样的方法,朝-y方向看xoz平面。 通过把每个人跟点Q的像Q‘相连,我们发现,左右边界并不对应着最左边与最右边的人。而且随着询问点的变化,对应着左右边界的人也在变化。

如果我们暴力的找对应左右边界的人,复杂度为n*q ,不可行。

我们发现,对应着左右边界的人虽然随着询问点变化,但他们都在凸包上(如下图)。

更进一步,如果我们将询问点按照x坐标排序,随着询问的x坐标增加,左右边界的人在凸包上的变化是顺时针旋转的。(考虑你从左到右观察一个正前方的凸包)

于是我们就能够通过一个nlogn的凸包预处理然后O(1)地回答每个询问,复杂度为O(q+nlogn)

剩下的是实现”从左到右看凸包时凸包左右边界的顺时针更新“。

首先是写out,in函数(右手法则,向外转就是逆时针),用来逆时针、顺时针遍历凸包上的点。因为极角排序凸包存的点是逆时针的(极角排序的那个角是与y轴的逆时针夹角。),所以out就是++。

先找到凸包的下上顶点,作为初始的左右边界的对应点。

然后根据x坐标从小到大枚举询问点Q。

对于每个Q,不断顺时针更新左边界对应的人,直到他与Q的连线在他凸包上顺时针的下一个人与Q的连线的外面(直观上显然正确)。 右边界同理。

某些编辑器比如codeforces不能混用iostream与stdio


代码

#include<cmath>
#include<algorithm>
#include<iostream>
#include<stdio.h>
#include<map>
#include<string.h>
#include<queue>
#include<stack>
using namespace std;
#define debug(x) cerr<<#x<<" = "<<(x)<<endl
#define rep(i,j,k) for(int i = (int)j;i <= (int)k;i ++)
#define FAST_IO ios_base::sync_with_stdio(false); cin.tie(nullptr)
//#define double long long
typedef long long ll;
typedef double db;
const int maxn = 2e5 + 5;
const db eps = 1e-7;
int n, q, a; long double ans[maxn];
bool eq(double a, double b) { return abs(a - b) <= eps; }
struct V {
double x, y;
V(double a = 0.0, double b = 0.0) :x(a), y(b) {}
void sc() { scanf("%lf%lf", &x, &y); }
double operator |(V const &o)const {
return x * o.y - o.x * y;
}
bool operator <(V const &o)const {
if (eq(x, o.x))return y < o.y;
return x < o.x;
}
V operator -(V const &o)const { return V(x - o.x, y - o.y); }
void pr() { printf("%lf %lf\n", x, y); }
}st[maxn];
pair<V, int> Q[maxn];
bool cmpr(V const &a, V const &b) {
V v1 = a - st[0], v2 = b - st[0];
return (v1 | v2) < -eps;
} vector<V> ch;
void getCH() {
sort(st + 1, st + n, cmpr);
ch.push_back(st[0]);
rep(i, 1, n-1) {
while (ch.size() > 1 && (st[i] - ch.back() | ch.back() - ch[ch.size() - 2]) < eps)ch.pop_back();
ch.push_back(st[i]);
}
}
int out(int x) { return x ? x - 1 : ch.size() - 1; }
int in(int x) { return x + 1 == (int)ch.size() ? 0 : x + 1; } double getx(double xs, double z, double xq) {
return xs + z / (z + a) * (xq - xs);
} int main() {
//FAST_IO;
int t;
cin >> t;
while (t--) { cin >> n >> a; rep(i, 0, n - 1)st[i].sc();
db zmn = st[0].y, zmx = st[0].y;
rep(i, 0, n - 1) {
zmn = min(zmn, st[i].y);
zmx = max(zmx, st[i].y); }
rep(i, 0, n - 1)st[i].y = a - st[i].y;
sort(st, st + n);
ch.clear();
getCH(); cin >> q;
ll qx = 0, qy = 0;
rep(i, 1, q) {
Q[i].first.sc();
Q[i].second = i;
}
sort(Q + 1, Q + 1 + q); int lp = 0, rp = 0;
rep(i, 0, ch.size() - 1)ch[i].y = a - ch[i].y;
while (ch[in(rp)].y < ch[rp].y)rp = in(rp);
while (ch[out(lp)].y > ch[lp].y)lp = out(lp); rep(i, 1, q) {
while (true) {
int ni = in(lp);
if (getx(ch[ni].x, ch[ni].y, Q[i].first.x) < getx(ch[lp].x, ch[lp].y, Q[i].first.x))lp = ni;
else break;
}
while (true) {
int ni = in(rp);
if (getx(ch[ni].x, ch[ni].y, Q[i].first.x) > getx(ch[rp].x, ch[rp].y, Q[i].first.x))rp = ni;
else break;
}
double x = abs(getx(ch[rp].x, ch[rp].y, Q[i].first.x) - getx(ch[lp].x, ch[lp].y, Q[i].first.x));
double h = getx(0, zmx, Q[i].first.y) - getx(0, zmn, Q[i].first.y);
ans[Q[i].second] = (long double)x * h;
//printf("%.20lf\n", x*h);
}
rep(i, 1, q)printf("%.20lf\n", (double)ans[i]);
}
cin >> n;
} /*
1
3 3
-2 1
7 2
3 1
3
2 5
-2 4
8 10
*/

心路历程

当有两个以上的bug时你就炸了
wa0:FAST_IO codeforces以用就wa wa1:输出rep(i,1,q) not rep(i,1,n)//最后才发现
wa2:凸包板子里面下标从0开始。

Problem F Fair Chocolate-Cutting

题意

来源:ICPC Asia Regional Contest, Yokohama, 2018–12–09

让你求所有平分凸多边形的线段中最长和最短的那两条。

输出它们的长度。

题解

证明题+硬核平面几何

我们有结论:(目前不会证明)

某合法(平分多边形的)线段的端点在顶点上时一定是局部最长或最短的。

除此之外,当某合法线段垂直于任意“外角”平分线时,也是局部最短的

于是问题就转化成

1.对于每个顶点,找到过它的合法线段。

2.对于每个角(跳过一些线段的“相邻”两条线段的延长线交出的角),找到垂直它的平分线的合法线段。注意特判掉角不存在的情况,即两边平行。

对于1,只要从该顶点开始不断逆时针扫顶点,直到刚刚超过1/2的总面积。然后把上一个三角形分割一下。

对于2,直接用平面几何知识以及三角函数求解

代码

#include<complex>
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <map>
#include<algorithm>
#include<iomanip>
typedef long long ll;
typedef double db;
typedef long double ld; using namespace std;
typedef complex<double> point;
#define rep(i,j,k) for(int i = (int)j;i <= (int)k;i ++)
#define F(i,N)for(int i = 0;i < (int)(N);i ++)
#define Decimal fixed<<setprecision(20)
#define PB push_back
#define EB emplace_back
#define X real();
#define Y imag();
#define curr(PP,i) PP[i]
#define next(PP,i) PP[(i+1)%PP.size()]
#define diff(PP,i) (next(PP,i)-curr(PP,i))
const double eps = 1e-7;
const double PI = acos(-1.0);
const double INF = 1e18;
int n;
vector<point> V;
double sumS,nowS,mx,mn;
double cross(const point& a, const point& b) {
return imag(conj(a)*b);
}
double dot(const point& a, const point& b) {
return real(conj(a)*b);
}
double area(const vector<point>& p) {
double A = 0;
F(i, p.size()) A += cross(curr(p, i), next(p, i));
return A / 2.;
}
double triArea(int a, int b, int c) {
a %= n, b %= n, c %= n;
if (a == b || b == c || c == a)return 0;
vector<point> tmp = { V[a],V[b],V[c] };
return area(tmp);
}
double distancePP(point const &a, point const &b) {
return abs(a - b);
}
//seg
struct segment :public vector<point> {
segment(const point &a, const point &b) {
push_back(a); push_back(b);
}
};
point crosspoint(segment const &l,segment const &m) {
double A = cross(l[1] - l[0], m[1] - m[0]);
double B = cross(l[1] - l[0], l[1] - m[0]);
if (abs(A) < eps&&abs(B) < eps)return m[0];//same line
if (abs(A) < eps)return point(INF, INF);//parallel
return m[0] + B / A*(m[1] - m[0]);
}
double angle(const point& a, const point& b) {
auto tmp = abs(arg(a) - arg(b));
return min(tmp, 2 * PI - tmp);
}
double angle(const segment &s1, const segment &s2) {
return angle(s1[1] - s1[0], s2[1] - s2[0]);
}
point projection(const segment &l, const point &p) {
double t = dot(p - l[0], l[0] - l[1]) / norm(l[0] - l[1]);
return l[0] + t * (l[0] - l[1]);
}
double distanceLP(const segment &l, const point &p) {
return abs(p - projection(l, p));
}
bool intersectSP(const segment &s,const point &p) {//point on seg
return abs(s[0] - p) + abs(s[1] - p) - abs(s[1] - s[0]) < eps;
}
double seica(int a,int b) {
a %= n; b %= n;
vector<point> tmp;
int i = a;
while (1) {
tmp.EB(V[i]);
if (i == b)break;
i = (i + 1) % n;
}
if (tmp.size() < 3)return 0;
else return area(tmp);
}
void checkV(int p,int q) {
p %= n; q %= n;
double remain = (sumS / 2) - nowS;
double rate = remain / triArea(p, q, q + 1);
point tmp = V[q] + rate * (V[(q + 1) % n] - V[q]);
mn = min(mn, distancePP(V[p], tmp));
mx = max(mx, distancePP(V[p], tmp));
}
void checkSeg(int p, int q) {
p %= n, q %= n;
nowS -= triArea(p, p + 1, q);
segment a(V[(p + 1) % n], V[p]);
segment b(V[(q) % n], V[(q + 1) % n]);
point crossp = crosspoint(a, b);
if (crossp == point(INF, INF)) {
double sqS = area({ V[p],V[(p + 1) % n],V[q],V[(q + 1) % n] });
if (nowS + sqS < sumS / 2 + eps)return;
bool flag = 0;
if (intersectSP(b, projection(b, a[0])) || intersectSP(b, projection(b, a[1])))flag = 1;
if (intersectSP(a, projection(a, b[0])) || intersectSP(a, projection(a, b[1])))flag = 1;
if (!flag)return;
mn = min(mn, distanceLP(b, a[0]));
return;
}
double tmpS=0;
if (distancePP(crossp, a[0]) > distancePP(crossp, a[1]) + eps) {
tmpS = seica(q + 1, p);
swap(a[0], a[1]);
swap(b[0], b[1]);
}
else {
tmpS= seica(p + 1, q);
}
double theta = angle(a, b);
double T = distancePP(crossp, a[0])*distancePP(crossp, b[0])*sin(theta) / 2;
double leg = (2 * (T - tmpS) + sumS) / sin(theta);
bool flag = 0;
if (leg > distancePP(a[1], crossp)*distancePP(a[1], crossp) + eps)flag = 1;
if (leg + eps < distancePP(a[0], crossp)*distancePP(a[0], crossp))flag = 1;
if (leg > distancePP(b[1], crossp)*distancePP(b[1], crossp) + eps)flag = 1;
if (leg + eps < distancePP(b[0], crossp)*distancePP(b[0], crossp))flag = 1;
if (flag)return;
mn = min(mn, 2*sqrt(leg)*sin(theta / 2));
}
int main() {
mx = 0;
mn = 1e18;
cin >> n;
F(i, n) {
db a, b;
cin >> a >> b;
V.EB(a, b);
}
int q = 1;
sumS = area(V);
F(i, n) {
while (nowS + triArea(i, q, q + 1) < sumS / 2 + eps) {
nowS += triArea(i, q, q + 1);
q++;
}
checkV(i,q);
checkSeg(i,q);
}
cout << Decimal << mn << endl;
cout << Decimal << mx << endl;
cin >> n;
return 0;
} /*
4
0 0
10 0
10 10
0 10 3
0 0
6 0
3 10
*/

心路历程


关键是不会证明。
复数大法好啊!!!(其实还是慢5倍??下面那题)

B. Ali and Wi-Fi

题意

来源ACM International Collegiate Programming Contest, Tishreen Collegiate Programming Contest

100个圆,每个圆代表wifi覆盖区域,每个wifi有一个网速。你最多可以连入m个WiFi,网速为这些wifi的总和。

问你最大的网速是多少。


题解

想法题+圆模板

我们知道n个圆最多可以把平面分成n^2-n+2个区域(递推进阶题233)。所以可以考虑遍历每个区域,计算它被哪些圆。。。

我们发现维护一块区域很难,而要取最大的m个更难考虑。

于是考虑维护圆与圆的交点。我们发现每个区域都是圆交出来的,i.e.交点都在区域的边界上。

于是我们可以枚举交点来代替枚举区域,显然交点个数也是n^2级别的。(考虑圆内含的情况,所以要把每个圆心也枚举一遍)

这样枚举的另一个好处是可以很简单地处理m:

枚举每个交点,计算有几个圆包含它。然后直接在这些圆中取最大的m个。

复杂度为n^2nnlogn。

如果我们用优先队列来维护最大的m个,那么复杂度里可以去掉一个n。


代码

#include<cmath>
#include<algorithm>
#include<iostream>
#include<stdio.h>
#include<map>
#include<string.h>
#include<queue>
#include<stack>
#include<iomanip>
#include<complex>
using namespace std;
#define debug(x) cerr<<#x<<" = "<<(x)<<endl
#define rep(i,j,k) for(int i = (int)j;i <= (int)k;i ++)
#define F(i,N)for(int i = 0;i < (int)(N);i ++)
#define FAST_IO ios_base::sync_with_stdio(false); cin.tie(nullptr)
#define EB emplace_back
#define Decimal fixed<<setprecision(20)
#define fi first
#define se second
#define INF 1000000000 typedef long long ll;
typedef double db;
typedef long double ld;
#define double long double
//typedef complex<double> point;
#define X real()
#define Y imag()
const int maxn = 2e5 + 5;
const db eps = 1e-10;
const db PI = acos(-1.);
int dcmp(ld x)
{
if (fabs(x) < eps) return 0;
return x < 0 ? -1 : 1;
}
struct V {
ld x, y;
V() {}
void sc() { scanf("%lf%lf", &x, &y); }
V(double a, double b) : x(a), y(b) { }
V operator+(const V &o)const { return V(x + o.x, y + o.y); }
V operator-(const V &o)const { return V(x - o.x, y - o.y); }
double L() { return sqrt(x * x + y * y); }
V N() {
double l = L();
return V(x / l, y / l);
}
V rot(double th) { return V(x * cos(th) - y * sin(th), x * sin(th) + y * cos(th)); }
V operator*(double z) { return V(x * z, y * z); }
double operator*(const V &o) const { return x * o.x + y * o.y; }
double operator|(const V &o) const { return x * o.y - o.x * y; }
double operator==(const V &o) const { return dcmp(x - o.x) == 0 && dcmp(y - o.y) == 0; }
void pr() { printf("%lf %lf\n", x, y); }
};
vector<V> p;
typedef V point;
V conj(V p) { return V(p.x, -p.y); }
double arg(V p) { return atan2(p.y, p.x); }
double norm(V p) { return p.L()*p.L(); }
double abs(V p) { return p.L(); } double cross(const point& a, const point& b) {
return imag(conj(a)*b);
}
double dot(const point& a, const point& b) {
return real(conj(a)*b);
}
struct circle
{
point c;
ld r, v;
circle() {}
circle(point c, ld r, ld v) :c(c), r(r), v(v) {} inline point pt(double a)
{
return point(c.x + cos(a)*r, c.y + sin(a)*r);
}
};
int getCircleCircleIntersection(circle C1, circle C2, point &t1, point &t2)
{
ld d = abs(C1.c - C2.c);
if (dcmp(d) == 0)
{
if (dcmp(C1.r - C2.r) == 0) return -1;// same cir
return 0;// include
}
if (dcmp(C1.r + C2.r - d) < 0) return 0;// disjoint
if (dcmp(fabs(C1.r - C2.r) - d) > 0) return 0;// include
ld a = arg(C2.c - C1.c);
ld da = acos((C1.r*C1.r + d * d - C2.r*C2.r) / (2 * C1.r*d));
point p1 = C1.pt(a - da), p2 = C1.pt(a + da);
t1 = p1;
if (p1 == p2) return 1;
t2 = p2;
return 2;
}
vector<circle> A;
circle aa[105];
vector<point> P;
point pp[10105];
priority_queue<int > Q;
int ans, sum;
int n, m, tot;
int main() {
FAST_IO;
int t; cin >> t;
while (t--) { tot = ans = 0;
cin >> n >> m;
F(i, n) { cin >> aa[i].c.x >> aa[i].c.y >> aa[i].r >> aa[i].v; }
F(i, n) {
rep(j, i + 1, n - 1) {
point t1, t2;
int num = getCircleCircleIntersection(aa[i], aa[j], t1, t2);
if (num == 1) {
pp[tot++] = t1;
}
else if (num == 2) {
pp[tot++] = t1;
pp[tot++] = t2;
}
}
}
F(i, n)pp[tot++] = aa[i].c;
F(j, tot) {
F(i, n) {
if (dcmp(abs(pp[j] - aa[i].c) - aa[i].r) <= 0) {
Q.push(-aa[i].v);
if (Q.size() > m)Q.pop();
}
}
sum = 0;
while (!Q.empty()) {
sum -= Q.top();
Q.pop();
}
ans = max(ans, sum);
}
cout << ans << endl;
}
cin >> t;
return 0;
}

心路历程

没初始化。
瞎debug发现long double int求和会爆精度。
priority_Q太久没用,忘记pop是从头上弹出的了orz。
vector 和对于1e6的数据其实只比数组慢10ms。

几何入门合集 gym101968 problem F. Mirror + gym102082 Problem F Fair Chocolate-Cutting + gym101915 problem B. Ali and Wi-Fi的更多相关文章

  1. JMeter入门合集

    JMeter从入门到精通 http://blog.csdn.net/lihengxin/article/details/4325918 jmeter入门教程- Jmeter教程及技巧汇总 http:/ ...

  2. UVa 12304 (6个二维几何问题合集) 2D Geometry 110 in 1!

    这个题能1A纯属运气,要是WA掉,可真不知道该怎么去调了. 题意: 这是完全独立的6个子问题.代码中是根据字符串的长度来区分问题编号的. 给出三角形三点坐标,求外接圆圆心和半径. 给出三角形三点坐标, ...

  3. Sorl 4.10 入门合集

    Sorl4.10 + Tomcat 7.0  win7环境下的安装 1.首先是到apache官网下载sorl 4.10 ,解压 2.进入路径\solr-4.10.4\example\webapps,拷 ...

  4. 转:android studio入门合集

    http://blog.csdn.net/column/details/zsl-androidstudio.html

  5. 深度学习优质学习项目大放送!-AI Studio精选开源项目合集推荐

    近期 在AI Studio上发现了不少优质的开源深度学习项目,从深度学习入门到进阶,涵盖了CV.NLP.生成对抗网络.强化学习多个研究方向,还有最新的动态图,都以NoteBook的方式直接开源出来,并 ...

  6. react-native 入门资源合集

    # 了解react-native React Native enables you to build world-class application experiences on native pla ...

  7. dotnet 从入门到放弃的 500 篇文章合集

    本文是记录我从入门到放弃写的博客 博客包括 C#.WPF.UWP.dotnet core .git 和 VisualStudio 和一些算法,所有博客使用 docx 保存 下载:dotnet 从入门到 ...

  8. [转载]HTML5开发入门经典教程和案例合集(含视频教程)

    http://www.iteye.com/topic/1132555 HTML5作为下一代网页语言,对Web开发者而言,是一门必修课.本文档收集了多个HTML5经典技术文档(HTML5入门资料.经典) ...

  9. Web 前端从入门菜鸟到实践老司机所需要的资料与指南合集

    http://web.jobbole.com/89188/ 2016 – 对于未来五年内Web发展的7个预测 2015 – 我的前端之路:从命令式到响应式,以及组件化与工程化的变革 怎么成为一名优秀的 ...

随机推荐

  1. C# 自动Ping服务

    using Newtonsoft.Json; using Newtonsoft.Json.Linq; using System; using System.Collections.Generic; u ...

  2. cocos2d-x实战 C++卷 学习笔记--第6章 场景与层

    前言: 一个场景(Scene)是由多个层(Layer)组成,而且层的个数要至少是1,不能为0. 场景切换 场景切换相关函数 1)void  runWithScene(Scene*  scene) 该函 ...

  3. 一步步学习Python-django开发-Mac下搭建Python-Django环境

    首先确定学习那一个版本,以之前对Python的接触,Python2.x 的版本仍然占据很大的市场.但作为初学者我还是决定学习新的版本,相信未来会逐渐转向3.x版本的. 其实搭建Python不是一件很难 ...

  4. 逻辑性最强的React Native环境搭建与调试

    React Native(以下简称RN),已经“火”了好一段时间了,网上的资料相对也很丰富,只是一直迟迟没有发布1.0,不过出身豪门(Facebook)的RN和国内顶级互联网公司对于RN的实践与应用, ...

  5. Chrome调试工具developer tool技巧

    Chrome这个浏览器赞的不能再赞了,给前端的开发调试工作带来了极大的效率提升. Chrome的简洁.快速吸引了无数人,它的启动速度.页面解析速度都很快,同时得益于Google V8的快速,Javas ...

  6. Asp.Net Core API网关Ocelot

    首先,让我们简单了解下什么是API网关? API网关是一个服务器,是系统的唯一入口.从面向对象设计的角度看,它与外观模式类似.API网关封装了系统内部架构,为每个客户端提供一个定制的API.它可能还具 ...

  7. [LeetCode] My Calendar II 我的日历之二

    Implement a MyCalendarTwo class to store your events. A new event can be added if adding the event w ...

  8. Docker学习笔记【三】安装Redis

    项目中使用到Redis,平常都是别人搭建的,今天试着在Google Cloud Platform 上搭建一个学习环境. 1.使用 docker pull redis 从docker hub中下载镜像 ...

  9. 【Idea】好的插件集合,持续更新

    UploadJar,用于配合Nexus上传jar包,方便上传 Key Promoter X,用于显示快捷键,学习快捷键非常实用 lombok,getter/setter使用注解,而不需要写 自动生成g ...

  10. think in java 读书笔记

    第三章 操作符 3.1 更简单的打印语句  原:System.out.println("打印"); 简单:print("我是更简单的");  =======需要 ...