(tarjan)洛谷 P2194 HXY烧情侣
思路：费用就是每个强连通分量里最小值的和，方案数就是每个强连通分量里最小值数量的乘积
众所周知，HXY 已经加入了 FFF 团。现在她要开始喜（sang）闻（xin）乐（bing）见（kuang）地烧情侣了。
这里有 nnn 座电影院，nnn 对情侣分别在每座电影院里，然后电影院里都有汽油，但是要使用它需要一定的费用。mmm 条单向通道连接相邻的两对情侣所在电影院。
HXY 有个绝技，如果她能从一个点开始烧，最后回到这个点，那么烧这条回路上的情侣的费用只需要该点的汽油费即可。并且每对情侣只需烧一遍，电影院可以重复去。然后她想花尽可能少的费用烧掉所有的情侣。
问：最少需要多少费用，并且当费用最少时的方案数是多少？由于方案数可能过大，所以请输出方案数对 109+710^9+7109+7 取模的结果。
（注：这里 HXY 每次可以从任何一个点开始走回路。就是说一个回路走完了，下一个开始位置可以任选。所以说不存在烧不了所有情侣的情况，即使图不连通，HXY 自行选择顶点进行烧情侣行动。且走过的道路可以重复走。）
输入格式
第一行一个正整数 nnn。
第二行 nnn 个正整数，表示每个点的汽油费 wiw_iwi。
第三行一个正整数 mmm。
接下来 mmm 行，每行两个正整数 xi,yix_i,y_ixi,yi，表示一条 xi→yix_i \to y_ixi→yi 的单向道路。
输出格式
输出一行两个整数，分别表示最小花费，和花费最小时的方案数。
输入输出样例
输入 #1
3 1 2 3 3 1 2 2 3 3 2
输出 #1
3 1
输入 #2
3 10 20 10 4 1 2 1 3 3 1 2 1
输出 #2
10 2
说明/提示
对于 30%30\%30% 的数据，1≤n,m≤201\le n,m \le 201≤n,m≤20；
对于另外 10%10\%10% 的数据，保证不存在回路；
对于 100%100\%100% 的数据，1≤n≤1051\le n \le 10^51≤n≤105，1≤m≤3×1051\le m \le 3\times 10^51≤m≤3×105，1≤wi≤1091\le w_i \le 10^91≤wi≤109。
#include
#include
#include
#include
#include
#include
#include
using namespace std;
#define ll long long
const int maxn=505000;
#define mod 1000000007
stackq;
setans;
vectora[maxn];
int cnt,n,m,head[maxn],dfn[maxn],low[maxn],co[maxn],tot,id,out[maxn],root;
struct Edge
{
int to,next;
}e[maxn];
void add(int u,int v)
{
e[cnt].to=v;
e[cnt].next=head[u];
head[u]=cnt++;
}
void init()
{
ans.clear();
for(int i=0;i=dfn[u]&&u!=root) ans.insert(u);
}
else if(!co[v])
low[u]=min(low[u],dfn[v]);
}
if(root==u&&ch>=2) ans.insert(u);
if(low[u]==dfn[u])
{
id++;
while(q.top()!=u)
{
int x=q.top();
q.pop();
co[x]=id;
}
co[u]=id;
q.pop();
}
}
int w[maxn];
int main()
{
scanf("%d",&n);
init();
for(int i=1;i<=n;i++) scanf("%d",&w[i]);
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
int u,v;
scanf("%d %d",&u,&v);
add(u,v);//add(v,u);
}
for(int i=1;i<=n;i++)
{
if(dfn[i]==0) {
root=i;
tarjan(i);
}
}
/*for(int i=1;i<=n;i++)
{
printf("co[%d]=%d\n",i,co[i]);
}*/
for(int i=1;i<=n;i++) a[co[i]].push_back(w[i]);
ll p1=0,p2=1;
for(int i=1;i<=id;i++)
{
ll p=0;
sort(a[i].begin(),a[i].end());
p1+=a[i][0];
for(int u:a[i])
{
if(u==a[i][0]) p++;
}
p2*=p;p2%=mod;
}
printf("%lld %lld\n",p1,p2);
}