loj2587铁人两项

无向图，图中选出定点三元组（a,b,c），a->b->c的路径没有重复边。问方案有多少？

————————————————————————————————————————————

首先求出圆方树，方点权值为连接的圆点数量，圆点权值为-1

这时，枚举a,c点，b点的方案数为a,c路径上的点权和。

枚举a,c点然后计算点权和明显超时。于是我们枚举b点，计算通过它的方案数。

所以神搜后有三种可能，b点的各个子分支与其它点，b点的向上分支与其它点，如果b点为a或c点的情况

————————————————————————————————————————————

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 const int maxn = 2e5 + 10, maxm = 8e5 + 10;
 4 int n, m, cnt;
 5 struct edge {
 6     int u, v, nxt;
 7 } e[maxm], ee[maxm];
 8 int head[maxn], js, headd[maxn], jss;
 9 void addage(int u, int v) {
10     e[++js].u = u;
11     e[js].v = v;
12     e[js].nxt = head[u];
13     head[u] = js;
14 }
15 void addagee(int u, int v) {
16     ee[++jss].u = u;
17     ee[jss].v = v;
18     ee[jss].nxt = headd[u];
19     headd[u] = jss;
20 }
21 int dfn[maxn], low[maxn], tim, sta[maxn], top, pw[maxn];
22 int sum, siz[maxn];
23 bool vis[maxn];
24 void tarjan(int u, int fa) {
25     dfn[u] = low[u] = ++tim;
26     sta[++top] = u;
27     for (int i = head[u]; i; i = e[i].nxt) {
28         int v = e[i].v;
29         if (!dfn[v]) {
30             tarjan(v, u);
31             low[u] = min(low[v], low[u]);
32             if (low[v] >= dfn[u]) {
33                 cnt++;
34                 int tp = 0;
35                 while (tp != v) {
36                     tp = sta[top--];
37                     addagee(cnt, tp);
38                     addagee(tp, cnt);
39                     pw[cnt]++;
40                     pw[tp] = -1;
41                 }
42                 addagee(cnt, u);
43                 addagee(u, cnt);
44                 pw[cnt]++;
45                 pw[u] = -1;
46             }
47         } else if (v != fa)
48             low[u] = min(low[u], dfn[v]);
49     }
50 }
51 long long ans;
52 void dfs1(int u, int fa) {
53     vis[u] = 1;
54     siz[u] = (u <= n);
55     for (int i = headd[u]; i; i = ee[i].nxt) {
56         int v = ee[i].v;
57         if (v == fa)
58             continue;
59         dfs1(v, u);
60         siz[u] += siz[v];
61     }
62 }
63 void dfs2(int u, int fa) {
64     for (int i = headd[u]; i; i = ee[i].nxt) {
65         int v = ee[i].v;
66         if (v == fa)
67             continue;
68         dfs2(v, u);
69         ans += (long long)pw[u] * (sum - siz[v]) * siz[v];  // u的每一个子分支与其它点的
70     }
71     ans += (long long)pw[u] * (sum - siz[u]) * siz[u];  // u向上的分支与其它的点
72     if (u <= n)
73         ans += (long long)pw[u] * (sum - 1);  // u作为端点
74 }
75 int main() {
76     scanf("%d%d", &n, &m);
77     for (int u, v, i = 0; i < m; ++i) {
78         scanf("%d%d", &u, &v);
79         addage(u, v);
80         addage(v, u);
81     }
82     cnt = n;
83     for (int i = 1; i <= n; ++i)
84         if (!dfn[i])
85             tarjan(i, 0);
86     for (int i = 1; i <= n; ++i)
87         if (!vis[i]) {
88             dfs1(i, 0);
89             sum = siz[i];
90             dfs2(i, 0);
91         }
92     cout << ans;
93
94     return 0;
95 }