思路:费用就是每个强连通分量里最小值的和,方案数就是每个强连通分量里最小值数量的乘积
众所周知,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 <cstdio>
#include<iostream>
#include <cstring>
#include <algorithm>
#include <stack>
#include<set>
#include<vector>
using namespace std;
#define ll long long
const int maxn=505000;
#define mod 1000000007
stack<int>q;
set<int>ans;
vector<int>a[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<maxn;i++) a[i].clear();
memset(out,0,sizeof(out));
memset(co,0,sizeof(co));
memset(head,-1,sizeof(head));
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
while(!q.empty()) q.pop();
cnt=0;id=0,tot=0;
}
void tarjan(int u)
{
q.push(u);
dfn[u]=low[u]=++tot;int ch=0;
for(int i=head[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
if(!dfn[v])
{
ch++;
tarjan(v);
low[u]=min(low[u],low[v]);
if(low[v]>=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);
}