【m*sqrt(m)暴力三元环】hdu 6184

xiaoxiao2021-02-28  41

一道赛场上读错题意的题。。。

传送门: http://acm.hdu.edu.cn/showproblem.php?pid=6184

题意: 给一张图,问能组成多少个,以同一条边构成的两个三元环(就是类似四边形,然后连一条对角线)。

思路: 听说叉姐的camp的第一场有m*sqrt(m)得方法。没去看camp血亏啊orz。

暴力枚举每一条边(也就是枚举两个端点x,y,先枚举x再枚举y,枚举x的时候可以记录一下x可到达的z点)。 对于第三个端点,如果后枚举的端点y的边数要小于sqrt(m),则枚举它。 否则枚举前一个端点x的第三个端点。

这样做的极端情况。 对于第一种枚举端点y的第三个端点z。由于前面枚举x的时候,已经顺带的把x和y的关系已经记录下来了,且因为满足小于sqrt(m),所以时间复杂度为m * sqrt(m)。

对于第二种枚举端点x的第三个端点z,最差的情况下为完全图。题目给的边数是2*10^5。sqrt一下最多也就100+个点。也是sqrt(m) * sqrt(m) * sqrt(m) 再上多个set的查询log(m)。 且这种情况数较少。所以复杂度差不多能约为m * sqrt(m)。如果第一种情况也用set来找的话,会给卡常数。

#include <bits/stdc++.h> #define MAXN 100005 #define ll long long using namespace std; vector<int> e[MAXN]; bool mark[MAXN]; set<ll> st; int link[MAXN]; ll cal(ll a) { return a * (a - 1) / 2; } void init() { memset(mark, false, sizeof(mark)); memset(link, 0, sizeof(link)); for (int i = 0; i < MAXN; i++) { e[i].clear(); } st.clear(); } void addEdge(int u, int v) { e[u].push_back(v); } int main() { int n, m, u, v; while (~scanf("%d %d", &n, &m)) { init(); int b = sqrt(m); for (int i = 0; i < m; i++) { scanf("%d %d", &u, &v); addEdge(u, v); addEdge(v, u); st.insert((ll)u * (n + 1) + v); st.insert((ll)v * (n + 1) + u); } ll ans = 0; for (int i = 1; i <= n; i++) { int x = i; mark[i] = true; for (int j = 0; j < e[x].size(); j++) { link[e[x][j]] = x; } for (int j = 0; j < e[x].size(); j++) { //暴力枚举每一条边 int y = e[x][j]; ll sum = 0; if (mark[y]) { continue; } if (e[y].size() <= b) { //选择一个边少的进行枚举sqrt(m) for (int k = 0; k < e[y].size(); k++) { //如果枚举的是y点的边,判断x是否有连边 int z = e[y][k]; if (link[z] == x) { //这边点集少,只需要跑m*sqrt(m),如果用set的话还要加上一个log(m)导致tle sum++; } } } else { for (int k = 0; k < e[x].size(); k++) { //如果枚举的是x点的边,判断是否在set里 int z = e[x][k]; if (st.find((ll)z * (n + 1) + y) != st.end()) { sum++; } } } ans += cal(sum); } } printf("%lld\n", ans); } }
转载请注明原文地址: https://www.6miu.com/read-78671.html

最新回复(0)