题目
Description
给定一张无重边、自环、割点的平面图,你需要回答 Q 次询问,每次询问会给出一个简单环,你需要回答在由这个简单环围成的多边形内部(包括边界上)的点有多少个。
保证图中每条线段不会通过除这条线段端点以外的其他点。
请仔细阅读题目读入输出的格式。
Input
第一行两个正整数 n, m 表示平面图有 n 个点和 m 条边。
接下来 m 行,每行两个正整数 a, b 表示平面图上的一条边 (a, b)。
接下来 n 行,每行两个整数 xi , yi 表示第 i 个点的坐标。
接下来一行一个整数 Q 表示有 Q 组询问。
接下来 Q 行每行表示一个询问。
每一行第一个整数 K 表示询问的环的大小,
接下来 K 个整数,按照顺时针或者逆时针的顺序给出环上的点。
保证读入数据一定有至少一个对应的 AC 自动机。
Output
Q 行,每行输出一个整数表示第 i 个询问的答案。
Sample Input
5 8
3 4
5 2
5 1
1 4
2 4
5 3
1 3
2 3
1 5
4 -5
4 -3
5 2
-5 -2
4
4 5 2 3 1
3 3 1 4
3 5 2 3
4 2 3 1 4
Sample Output
4
3
3
4
Data Constraint
思路
考虑一个模型:每个点有一个单位的流量从这个点流出,流量沿着边传送,最后所有点流出的流量汇集到某个点。
那么每条边在某个方向上可能有一个流量,除了汇点其他所有点的流出流量-流入流量=1。
一个环当中的点的个数就是流出环的流量总和减去流入环的流量总和。
为了避免出现询问最外面的环的情况(在绝大多数数据中我加了这种情况),我们在这个图的最外面加一个点并向图
连边,以这个新加的点作为汇点。标程的做法是在x坐标最小的点左边连接这个点。
构造这个流量可以在图上随便搞一个以汇点为根的生成树,每个点往父亲的流量大小是这个点的子树大小。
注意到环的总长不超过3e5 ,所以每次绕环走一圈是可以做到的。
代码
#include <bits/stdc++.h>
#define ll long long
#define db double
using namespace std;
const int N = 1e5 + 10;
int n, m;
struct P{
int x,y;
P(int _x=0, int _y=0) {x=_x,y=_y;}
P operator-(P b){return P(x-b.x,y-b.y);}
P operator+(P b){return P(x+b.x,y+b.y);}
ll operator*(P b){return (ll)x*b.y-(ll)y*b.x;}
P operator*(ll r){return P(x*r,y*r);}
} p[N];
vector<pair<db,int> > e[N];
vector<int>sum[N];
map<int,int> loc[N];
int sz[N],vis[N];
int fa[N];
void bfs() {
static int Q[N], h = 0, t = 0;
Q[++t] = n + 1;
vis[n+1] = 1;
while(h < t){
int x = Q[++h];
for(auto y:e[x]){
if(!vis[y.second]){
Q[++t]=y.second;
fa[y.second] = x;
vis[y.second]=1;
}
}
}
assert(t == n + 1);
for(int i = n + 1; i; i--) {
int x = Q[i];
if(x == n + 1) continue;
sz[x]++;
sz[fa[x]] += sz[x];
sum[x].resize(e[x].size());
for(auto y:e[x]){
if(fa[y.second] == x){
sum[x][loc[x][y.second]]=-sz[y.second];
}
}
sum[x][loc[x][fa[x]]]=sz[x];
for(int i=1; i<sum[x].size(); i++)
sum[x][i]+=sum[x][i-1];
}
}
int t[N];
void wr(P a){
printf("%d %d\n",a.x,a.y);
}
void wr(vector<int> a){
for(int x:a)printf("%d ",x);
printf("\n");
}
int main() {
freopen("graph.in","r",stdin);
freopen("graph.out","w",stdout);
cin >> n >> m;
int mix=1e9;
for(int i = 1; i <= m; i++) {
int x,y;scanf("%d %d",&x,&y);
e[x].push_back(make_pair(0,y));
e[y].push_back(make_pair(0,x));
}
p[n+1]=P(-1e9,-1e9);
for(int i = 1; i <= n; i++) {
scanf("%d %d",&p[i].x,&p[i].y);
mix=min(mix,p[i].x);
}
for(int i = 1; i <= n; i++) if(p[i].x == mix) {
e[i].push_back(make_pair(0,n+1));
e[n+1].push_back(make_pair(0,i));
break;
}
for(int i = 1; i <= n + 1; i++) {
for(auto &y:e[i]){
y.first = atan2(p[y.second].y - p[i].y, p[y.second].x - p[i].x);
}
sort(e[i].begin(), e[i].end());
for(int j = 0; j < e[i].size(); j++) {
loc[i][e[i][j].second] = j;
}
}
bfs();
int q = 0; cin >> q;
for(int o = 1; o <= q; o++) {
int k; scanf("%d", &k);
int ans = 0;
for(int i = 1; i <= k; i++) scanf("%d",&t[i]);
ll u=0;
for(int i=2; i<k; i++)
u+=(p[t[i]]-p[t[1]])*(p[t[i+1]]-p[t[1]]);
if(u<0)reverse(t+1,t+1+k);
t[0]=t[k],t[k+1]=t[1];
for(int i=1; i<=k; i++){
int l=loc[t[i]][t[i-1]],r=loc[t[i]][t[i+1]];
if(l<=r){
assert(l!=r);
ans+=sum[t[i]][r-1]-sum[t[i]][l];
}else{
ans+=sum[t[i]].back()-sum[t[i]][l];
if(r>0) ans+=sum[t[i]][r-1];
}
}
printf("%d\n",ans);
}
}