bzoj1179

其实是最大团问题，但是直接dfs不保证不tle

题解再次来自lsj（懒得敲):

对于一个强连通分量, 一定是整个走或者不走, 所以tarjan缩点然后跑dijkstra或spfa 

spfa

 1 #include<bits/stdc++.h>
 2 #define clr(a,x) memset(a,x,sizeof(a))
 3 #define rep(i,l,r) for(int i=l;i<r;i++)
 4 #define Rep(i,a) rep(i,0,e[a].size())
 5 #define REP(i,a) rep(i,0,e1[a].size())
 6 typedef long long ll;
 7 using namespace std;
 8 int read()
 9 {
10     char c=getchar();
11     int ans=0,f=1;
12     while(!isdigit(c)){
13         if(c=='-') f=-1;
14         c=getchar();
15     }
16     while(isdigit(c)){
17         ans=ans*10+c-'0';
18         c=getchar();
19     }
20     return ans*f;
21 }
22 const int maxn=500005,inf=0x7fffffff;
23 bool b[maxn],p[maxn];
24 int dfstime,cnt,s,k,n,m,w[maxn],low[maxn],pre[maxn],scc[maxn],d[maxn],v[maxn];
25 stack<int>S;
26 queue<int>Q;
27 vector<int>e[maxn];
28 vector<int>e1[maxn];
29 void dfs(int a)
30 {
31     low[a]=pre[a]=++dfstime;
32     S.push(a);
33     Rep(i,a){
34         if(!pre[e[a][i]]){
35             dfs(e[a][i]);
36             low[a]=min(low[a],low[e[a][i]]);
37         }else if(!scc[e[a][i]]){
38             low[a]=min(low[a],pre[e[a][i]]);
39         }
40     }
41     if(low[a]==pre[a]){
42         ++cnt;
43         int x=0;
44         while(x!=a){
45             x=S.top();
46             S.pop();
47             scc[x]=cnt;
48             w[cnt]+=v[x];
49             b[cnt]|=p[x];
50         }
51     }
52 }
53 void spfa()
54 {    
55     clr(d,0);
56     d[scc[s]]=w[scc[s]];
57     Q.push(scc[s]);
58     while(!Q.empty()){
59         int now=Q.front();
60         Q.pop();
61         REP(i,now){
62             if(d[e1[now][i]]<d[now]+w[e1[now][i]]){
63                 d[e1[now][i]]=d[now]+w[e1[now][i]];
64                 Q.push(e1[now][i]);
65             }
66         }
67     }
68 }
69 int main()
70 {
71     n=read(),m=read(),cnt=dfstime=0,clr(low,0),clr(pre,0),clr(scc,0),clr(b,0),clr(p,0);
72     rep(i,0,m){
73         int from=read(),to=read();
74         e[from].push_back(to);
75     }
76     rep(i,1,n+1) v[i]=read();
77     s=read(),k=read();
78     while(k--){
79         int t=read();
80         p[t]=1;
81     }
82     rep(i,1,n+1){
83         if(!scc[i]) dfs(i);
84     }
85 
86     rep(i,1,n+1){
87         Rep(j,i){
88             if(scc[i]!=scc[e[i][j]])
89             e1[scc[i]].push_back(scc[e[i][j]]);
90         }
91     }
92     spfa();
93     int ans=0;
94     rep(i,1,cnt+1){
95         if(b[i]) ans=max(ans,d[i]);
96     }
97     printf("%d\n",ans);
98     return 0;
99 }
View Code

dijkstra

  1 #include<bits/stdc++.h>
  2 #define clr(a,x) memset(a,x,sizeof(a))
  3 #define rep(i,l,r) for(int i=l;i<r;i++)
  4 #define Rep(i,a) rep(i,0,e[a].size())
  5 #define REP(i,a) rep(i,0,e1[a].size())
  6 typedef long long ll;
  7 using namespace std;
  8 int read()
  9 {
 10     char c=getchar();
 11     int ans=0,f=1;
 12     while(!isdigit(c)){
 13         if(c=='-') f=-1;
 14         c=getchar();
 15     }
 16     while(isdigit(c)){
 17         ans=ans*10+c-'0';
 18         c=getchar();
 19     }
 20     return ans*f;
 21 }
 22 struct node{
 23     int d,num;
 24     inline bool operator <(const node&A)const{
 25         return d>A.d;
 26     }
 27 };
 28 const int maxn=500005,inf=0x7fffffff;
 29 bool b[maxn],p[maxn];
 30 int dfstime,cnt,s,k,n,m,w[maxn],low[maxn],pre[maxn],scc[maxn],d[maxn],v[maxn];
 31 stack<int>S;
 32 vector<int>e[maxn];
 33 vector<int>e1[maxn];
 34 priority_queue<node>Q;
 35 void dfs(int a)
 36 {
 37     low[a]=pre[a]=++dfstime;
 38     S.push(a);
 39     Rep(i,a){
 40         if(!pre[e[a][i]]){
 41             dfs(e[a][i]);
 42             low[a]=min(low[a],low[e[a][i]]);
 43         }else if(!scc[e[a][i]]){
 44             low[a]=min(low[a],pre[e[a][i]]);
 45         }
 46     }
 47     if(low[a]==pre[a]){
 48         ++cnt;
 49         int x=0;
 50         while(x!=a){
 51             x=S.top();
 52             S.pop();
 53             scc[x]=cnt;
 54             w[cnt]+=v[x];
 55             b[cnt]|=p[x];
 56         }
 57     }
 58 }
 59 void dijkstra()
 60 {    
 61     clr(d,0);
 62     d[scc[s]]=w[scc[s]];
 63     node start;
 64     start.num=scc[s],start.d=w[scc[s]];
 65     Q.push(start);
 66     while(!Q.empty()){
 67         node now=Q.top();
 68         Q.pop();
 69         if(now.d==d[now.num]){
 70             REP(i,now.num){
 71                 if(d[now.num]+w[e1[now.num][i]]>d[e1[now.num][i]]){
 72                     d[e1[now.num][i]]=d[now.num]+w[e1[now.num][i]];
 73                     node next;
 74                     next.d=d[e1[now.num][i]];
 75                     next.num=e1[now.num][i];
 76                     Q.push(next);
 77                 }
 78             }
 79         }
 80     }
 81 }
 82 int main()
 83 {
 84     n=read(),m=read(),cnt=dfstime=0,clr(low,0),clr(pre,0),clr(scc,0),clr(b,0),clr(p,0);
 85     rep(i,0,m){
 86         int from=read(),to=read();
 87         e[from].push_back(to);
 88     }
 89     rep(i,1,n+1) v[i]=read();
 90     s=read(),k=read();
 91     while(k--){
 92         int t=read();
 93         p[t]=1;
 94     }
 95     rep(i,1,n+1){
 96         if(!scc[i]) dfs(i);
 97     }
 98     rep(i,1,n+1){
 99         Rep(j,i){
100             if(scc[i]!=scc[e[i][j]])
101             e1[scc[i]].push_back(scc[e[i][j]]);
102         }
103     }
104     dijkstra();
105     int ans=0;
106     rep(i,1,n+1){
107         if(b[i]) ans=max(ans,d[i]);
108     }
109     printf("%d\n",ans);
110     return 0;
111 }
View Code

1179: [Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 1921  Solved: 763
[Submit][Status][Discuss]

Description

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号

Output

输出一个整数，表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。

Sample Input

6 7
1 2
2 3
3 5
2 4
4 1
2 6
6 5
10
12
8
16
1 5
1 4
4
3
5
6

Sample Output

47

HINT

 

50%的输入保证N, M<=3000。所有的输入保证N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

 

Source

 
[Submit][Status][Discuss]
posted @ 2015-07-27 23:47  ChenThree  阅读(173)  评论(1)  编辑  收藏  举报