4541: [Hnoi2016]矿区
学习了一下平面图剖分的姿势,orz cbh
每次只要随便选择一条边,然后不停尽量向左转就行
#include <bits/stdc++.h>
#define N 1300000
#define M 5000013
#define LL long long
#define pb push_back
using namespace std;
LL n, m, k;
struct point
{
LL x, y;
} S[N];
vector <LL> bi[N];
vector <LL> re[N], vis[N], s1[N], s2[N];
vector <LL> bt[N], ca[N], cb[N];
LL tot_are;
LL siz[N], sum1[N], sum2[N], tot[N];
LL nwc;
LL comp(LL a, LL b)
{
return atan2(S[a].y - S[nwc].y, S[a].x - S[nwc].x) > atan2(S[b].y - S[nwc].y, S[b].x - S[nwc].x);
}
LL getsiz(point a, point b, point c)
{
return (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y);
} namespace addr
{
LL hs[M]; LL dt[M];
LL get(LL a, LL b)
{
return 1ll * a * + b;
}
LL & find(LL a, LL b)
{
LL p = get(a, b); LL q = p % M;
while (hs[q] && hs[q] != p) q = (q + ) % M;
hs[q] = p; return dt[q];
}
}
LL tvis[N], fa[N];
void dfs(LL t)
{
//cout << t << "\n"; tvis[t] = ;
sum1[t] = siz[t] * siz[t];
sum2[t] = siz[t] * ;
for (LL i = ; i < bt[t].size(); ++ i)
if (!tvis[bt[t][i]])
{
dfs(bt[t][i]); fa[bt[t][i]] = t;
sum1[t] += sum1[bt[t][i]];
sum2[t] += sum2[bt[t][i]];
s1[ca[t][i]][addr :: find(ca[t][i], cb[t][i])] = sum1[bt[t][i]];
s2[ca[t][i]][addr :: find(ca[t][i], cb[t][i])] = sum2[bt[t][i]];
s1[cb[t][i]][addr :: find(cb[t][i], ca[t][i])] = -sum1[bt[t][i]];
s2[cb[t][i]][addr :: find(cb[t][i], ca[t][i])] = -sum2[bt[t][i]];
}
}
LL read()
{
LL x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
LL c = ;
int main()
{
//freopen("mine2.in", "r", stdin);
n = read(); m = read(); k = read();
for (LL i = ; i <= n; ++ i) S[i].x = read(), S[i].y = read();
for (LL i = ; i <= m; ++ i)
{
LL a, b;
a = read(); b = read();
bi[a].push_back(b);
bi[b].push_back(a);
}
for (LL i = ; i <= n; ++ i) re[i].resize(bi[i].size()), s1[i] = s2[i] = vis[i] = re[i];
for (LL i = ; i <= n; ++ i) nwc = i, sort(bi[i].begin(), bi[i].end(), comp);
for (LL i = ; i <= n; ++ i)
for (LL j = ; j < bi[i].size(); ++ j)
addr :: find(i, bi[i][j]) = j;
for (LL i = ; i <= n; ++ i)
for (LL j = ; j < bi[i].size(); ++ j)
re[i][j] = addr :: find(bi[i][j], i); for (LL i = ; i <= n; ++ i)
for (LL j = ; j < bi[i].size(); ++ j)
if (!vis[i][j])
{
tot_are ++;
for (LL k = i, p = j; !vis[k][p]; )
{
siz[tot_are] += getsiz(S[i], S[k], S[bi[k][p]]);
vis[k][p] = tot_are;
LL np = (re[k][p] + ) % bi[bi[k][p]].size();
k = bi[k][p];
p = np;
}
if (siz[tot_are] < ) c = tot_are;
}
for (LL i = ; i <= n; ++ i)
for (LL j = ; j < bi[i].size(); ++ j)
bt[vis[i][j]].push_back(vis[bi[i][j]][re[i][j]]),
ca[vis[i][j]].push_back(i),
cb[vis[i][j]].push_back(bi[i][j]);
dfs(c);
for (LL i = , last = ; i <= k; ++ i)
{
LL d, ns1 = , ns2 = ;
d = read(); d = (d + last) % n + ;
LL fs, ls, nw;
fs = read(); fs = (fs + last) % n + ; ls = fs;
for (LL j = ; j <= d; ++ j, ls = nw)
{
nw = read(); nw = (nw + last) % n + ;
ns1 += s1[ls][addr :: find(ls, nw)];
ns2 += s2[ls][addr :: find(ls, nw)];
}
ns1 += s1[ls][addr :: find(ls, fs)];
ns2 += s2[ls][addr :: find(ls, fs)];
LL g = __gcd(ns1, ns2);
cout << ns1 / g << " " << ns2 / g << "\n";
last = ns1 / g;
}
}
放在class里的东西还会爆栈QAQ,以后不敢用了
4541: [Hnoi2016]矿区的更多相关文章
- BZOJ 4541: [Hnoi2016]矿区 平面图转对偶图+DFS树
4541: [Hnoi2016]矿区 Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 433 Solved: 182[Submit][Status][ ...
- ●BZOJ 4541 [Hnoi2016]矿区
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=4541 题解: 平面图的对偶图,dfs树 平面图的对偶图的求法: 把所有双向边拆为两条互为反向 ...
- bzoj 4541: [Hnoi2016]矿区【平面图转对偶图+生成树】
首先平面图转对偶图,大概思路是每条边存正反,每个点存出边按极角排序,然后找每条边在它到达点的出边中极角排序的下一个,这样一定是这条边所属最小多边形的临边,然后根据next边找出所有多边形,用三角剖分计 ...
- [HNOI2016]矿区
[HNOI2016]矿区 平面图转对偶图 方法: 1.分成正反两个单向边,每个边属于一个面 2.每个点按照极角序sort出边 3.枚举每一个边,这个边的nxt就是反边的前一个(这样找到的是面的边逆时针 ...
- 【LG3249】[HNOI2016]矿区
[LG3249][HNOI2016]矿区 题面 洛谷 题解 先平面图转对偶图, 建好了对偶图之后随意拿出一个生成树,以无边界的范围为根. 无边界的范围很好求,用叉积算出有向面积时,算出来是负数的就是无 ...
- BZOJ4541 [Hnoi2016]矿区
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...
- 【bzoj4541】 Hnoi2016—矿区
http://www.lydsy.com/JudgeOnline/problem.php?id=4541 (题目链接) 题意 给出一个平面图,若干询问,每次询问一个凸多边形内小多边形面积的平方和与面积 ...
- BZOJ4541 HNOI2016矿区(平面图转对偶图)
考虑先将平面图转化为对偶图.具体地,将无向边拆成两条有向边.每次考虑找到包围一个区域的所有边.对当前考虑的边,找到该边的反向边在该边终点的出边集中,按极角序排序的后继,这条后继边也是包围该区域的边.这 ...
- [BZOJ4541][HNOI2016]矿区(平面图转对偶图)
https://www.cnblogs.com/ljh2000-jump/p/6423399.html #include<cmath> #include<vector> #in ...
随机推荐
- [LeetCode] Longest Palindrome 最长回文串
Given a string which consists of lowercase or uppercase letters, find the length of the longest pali ...
- StringIO和BytesIO
1. StringIO 很多时候,数据读写不一定是文件,也可以在内存中读写. StringIO顾名思义就是在内存中读写str. 要把str写入StringIO,我们需要先创建一个StringIO,然后 ...
- 《UNIX环境高级编程》笔记——2.标准和实现
随着UNIX各种衍生版本不断发展壮大,标准化工作就十分必要.其实干啥事都是这样,玩的人多了,必须进行标准化. 一.UNIX标准 1.1 ISO C(ANSI C) ANSI:Amerocan Nato ...
- 推公式hdu2298
(0,0)点到(x,y)以v射箭 问 最小的角度 不能射到 -1 x=v*cos(a)*t; y=v*sin(a)*t-0.5*g*t*t; 消去t; y=x*tan(a)-0.5*g*x*x/(v* ...
- 敏捷组织中PMO应遵循的准则
敏捷改变了人们的工作方式,不仅仅是开发部门,而且还包括其它的部门,例如HR.财务以及PMO等.在大多数组织中,PMO是一个控制体.它指导项目团队的规范.模板以及流程.目前,大多数的IT组织都敏捷化了. ...
- IRandomAccessStream, IBuffer, Stream, byte[] 之间相互转换
/* * 用于实现 IRandomAccessStream, IBuffer, Stream, byte[] 之间相互转换的帮助类 */ using System;using System.IO;us ...
- C#面向对象设计模式纵横谈——4.Builder 生成器模式(创建型模式)
动机 (Motivation) 在软件系统中,有时候面临着“一个复杂对象”的创建工作,其通常由各个部分的子对象用一定的算法构成:由于需求的变化,这个复杂对象的各个部分经常面临着剧烈的变化,但是它们组合 ...
- BZOJ 1221: [HNOI2001] 软件开发
1221: [HNOI2001] 软件开发 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1428 Solved: 791[Submit][Stat ...
- Zabbix监控redis status
概述 zabbix采用Trapper方式监控redis status 原理 redis-cli info命令得到redis服务器的统计信息,脚本对信息分两部分处理: (1)# Keyspace部分为Z ...
- UEFI模式安装Win10和Linux双系统
最近心血来潮,想装一个Linux.Windows双系统,吸取上次安装的经验,这次一定都要使用UEFI模式启动,需要注意的是必须是支持此种启动模式的系统(一般解压之后都有efi文件夹不需要刻录),这次遇 ...