LJOJ 1500:
题目：http://www.docin.com/p-601990756.html
Sol:贪心，从叶子结点往上加入无法传递了，就需要建设。 Dfs返回的是到达叶子节点最多所要的能量，如果大于最大能量就需要建设放大器.
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 const int Maxn=100100; 6 using namespace std; 7 int head[Maxn],vis[Maxn],cnt,Ans,x,v,w,n,Power; 8 struct EDGE{int to,next,w;}edge[Maxn<<2]; 9 inline void Add(int u,int v,int w) {edge[cnt].to=v;edge[cnt].next=head[u];edge[cnt].w=w;head[u]=cnt++;} 10 inline int Max(int x,int y) {return x>y?x:y;} 11 int Dfs(int u) 12 { 13 int Res=0; vis[u]=true; 14 for (int i=head[u];i!=-1;i=edge[i].next) 15 if (!vis[edge[i].to]) 16 { 17 int Tmp=Dfs(edge[i].to); 18 if (Tmp+edge[i].w>=Power) {if (edge[i].w>Res) Res=edge[i].w; Ans++;} 19 else if (Tmp+edge[i].w>Res) Res=Tmp+edge[i].w; 20 } 21 return Res; 22 } 23 int main() 24 { 25 scanf("%d",&n); int Mx=0; 26 cnt=0; memset(head,-1,sizeof(head)); Ans=0; 27 for (int i=1;i<=n;i++) 28 { 29 scanf("%d",&x); 30 for (int j=1;j<=x;j++) 31 { 32 scanf("%d%d",&v,&w); 33 Add(i,v,w),Add(v,i,w); Mx=Max(Mx,w); 34 } 35 } 36 scanf("%d",&Power); 37 if (Mx>=Power) {puts("No solution."); return 0;} 38 Dfs(1); 39 printf("%d ",Ans); 40 return 0; 41 }
LJOJ 1501:
题目描述 出于最高安全性考虑，司令部采用了特殊的安全操作系统，该系统采用一个特殊的文件系统。在这个文件系统中所有磁盘空间都被分成了相同尺寸的N块，用整数1到N标识。每个文件占用磁盘上任意区域的一块或多块存储区，未被文件占用的存储块被认为是可使用的。如果文件存储在磁盘上自然连续的存储块中，则能被以最快的速度读出。 因为磁盘是匀速转动的，所以存取上面不同的存储块需要的时间也不同。读取磁盘开头处的存储块比读取磁盘结尾处的存储块快。根据以上现象，我们事先将文件按其存取频率的大小用整数1到K标识。按文件在磁盘上的最佳存储方法，1号文件将占用1,2,…,S1的存储块，2号文件将占用S1+1,S1+2,...,S1+S2的存储块，以此类推（Si是被第i个文件占用的存储块的个数）。为了将文件以最佳形式存储在磁盘上，需要执行存储块移动操作。一个存储块移动操作包括从磁盘上读取一个被占用的存储块至内存并将它写入其它空的存储块，然后宣称前一个存储块被释放，后一个存储块被占用。 本程序的目的是通过执行最少次数的存储块移动操作，将文件按最佳方式存储到磁盘上，注意同一个文件的存储块在移动之后其相对次序不可改变。 输入 每个磁盘说明的第一行包含两个用空格隔开的整数N和K，1<=K<=N<=100000，接下来的K行每行说明一个文件，对第i个文件的说明是这样的：首先以整数Si开头，表示第i个文件的存储块数量，1<=Si<=N-K，然后跟Si个整数，每个整数之间用空格隔开，表示该文件按自然顺序在磁盘上占用的存储块的标识。所有这些数都介于1和N之间，包括1和N。 一个磁盘说明中所有存储块的标识都是不同的，并且该盘至少有一个空的存储块。 输出 对于每一个磁盘说明，只需输出一行句子"We need M move operations"，M表示将文件按最佳方式存储到磁盘上所需进行的最少存储块移动操作次数。如果文件已按最佳方式存储，仅需输出“No optimization needed.”。 样例输入 20 3 4 2 3 11 12 1 7 3 18 5 10 样例输出 We need 9 move operations.
Sol:环需要有多余的一块进行交换，不然会覆盖
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 using namespace std; 5 int cnt,mb[100001],n,k,p,x,v[100001]; 6 int work(int x) 7 { 8 if (mb[x]==x) return 0; 9 int p=x,t=x,tot=0;v[p]=1; 10 while (1) 11 { 12 p=mb[p]; 13 tot++; 14 if (p!=t && p!=0 && v[p]) return tot; 15 v[p]=1; 16 if (p==0) return tot-1; 17 if (p==t) return tot+1; 18 19 } 20 } 21 int main() 22 { 23 scanf("%d%d",&n,&k); 24 for (int i=1;i<=k;i++) 25 { 26 scanf("%d",&p); 27 for (int i=1;i<=p;i++) 28 { 29 scanf("%d",&x); 30 mb[++cnt]=x; 31 } 32 } 33 int ans=0; 34 for (int i=1;i<=cnt;i++) if (!v[i]) ans+=work(i); 35 if (ans==0) printf("No optimization needed."); 36 else printf("We need %d move operations. ",ans); 37 }
LJOJ 1502:
题目:删去一条边使得最短路最长
Sol:Dij后找到最短路，依次删去最短路上的边进行最短路即可
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 #include <vector> 6 #include <queue> 7 #define Pa pair<int,int> 8 #define mp make_pair 9 using namespace std; 10 const int Inf=0x3f3f3f3f; 11 const int Maxn=1010; 12 const int Maxm=1000100; 13 priority_queue<Pa,vector<Pa>,greater<Pa> > Q; 14 struct EDGE{int to,next,w;}edge[Maxm<<2]; 15 int head[Maxn],Dis[Maxn],vis[Maxn],n,m,u,v,w,cnt,Pre[Maxn],Ans,Ban,Ban_Edge[Maxn]; 16 inline void Add(int u,int v,int w) {edge[cnt].to=v;edge[cnt].next=head[u];edge[cnt].w=w; head[u]=cnt++;} 17 int Dij(int Flag) 18 { 19 memset(vis,false,sizeof(vis)); 20 for (int i=1;i<=n;i++) Dis[i]=Inf; 21 Dis[1]=0; Q.push(mp(0,1)); 22 while (!Q.empty()) 23 { 24 int u=Q.top().second; Q.pop(); 25 if (vis[u]) continue; vis[u]=true; 26 for (int i=head[u];i!=-1;i=edge[i].next) 27 if (!vis[edge[i].to] && Dis[edge[i].to]>Dis[u]+edge[i].w && Ban!=i && Ban!=(i^1)) 28 { 29 Dis[edge[i].to]=Dis[u]+edge[i].w; 30 Q.push(mp(Dis[edge[i].to],edge[i].to)); 31 if (Flag) Pre[edge[i].to]=u,Ban_Edge[edge[i].to]=i; 32 } 33 } 34 return Dis[n]; 35 } 36 int main() 37 { 38 // freopen("c.in","r",stdin); 39 scanf("%d%d",&n,&m); memset(head,-1,sizeof(head)); 40 for (int i=1;i<=m;i++) scanf("%d%d%d",&u,&v,&w),Add(u,v,w),Add(v,u,w); 41 Ban=Inf; 42 Dij(true); Ans=0; 43 u=n; 44 while (true) 45 { 46 Ban=Ban_Edge[u]; 47 int Tmp=Dij(false); 48 if (Tmp>Ans) Ans=Tmp; 49 u=Pre[u]; 50 if (u==1) break; 51 } 52 printf("%d ",Ans); 53 return 0; 54 55 }
LJOJ 1503
题目:经典网络流:Pig
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 #define Inf 1<<29 6 using namespace std; 7 const int Maxn=110; 8 struct node 9 { 10 int to,next,w; 11 }edge[4000000]; 12 int cnt,Level[Maxn<<8],pigs[Maxn<<8],Q[4000000],head[Maxn<<8],n,m,Pre[Maxn<<8],limit; 13 inline void Add(int u,int v,int w) 14 { 15 edge[cnt].to=v; 16 edge[cnt].next=head[u]; 17 edge[cnt].w=w; 18 head[u]=cnt++; 19 edge[cnt].to=u; 20 edge[cnt].next=head[v]; 21 edge[cnt].w=0;//写成w 22 head[v]=cnt++; 23 } 24 25 inline int Min(int a,int b){return a>b?b:a;} 26 bool Bfs() 27 { 28 memset(Level,-1,sizeof(Level)); 29 int l=1,r=1; 30 Q[1]=0;Level[0]=0; 31 while (l<=r) 32 { 33 int u=Q[l++]; 34 for (int i=head[u];i!=-1;i=edge[i].next) 35 if (Level[edge[i].to]==-1 && edge[i].w>0) 36 { 37 Level[edge[i].to]=Level[u]+1; 38 Q[++r]=edge[i].to; 39 } 40 } 41 if (Level[n+1]>0) return true; 42 return false; 43 } 44 45 int Find(int u,int low) 46 { 47 if (u==n+1) return low; 48 int tmp=0; 49 for (int i=head[u];i!=-1;i=edge[i].next) 50 { 51 int v=edge[i].to; 52 if (edge[i].w>0 53 && Level[v]==Level[u]+1 54 && (tmp=Find(v,Min(low,edge[i].w)))) 55 { 56 edge[i].w-=tmp,edge[i^1].w+=tmp; 57 return tmp; 58 } 59 } 60 Level[u]=-1; 61 return 0; 62 } 63 64 int main() 65 { 66 scanf("%d%d",&m,&n); 67 memset(Pre,0,sizeof(Pre)); 68 memset(head,-1,sizeof(head)); 69 cnt=0; 70 for (int i=1;i<=m;i++) scanf("%d",&pigs[i]); 71 for (int i=1;i<=n;i++) 72 { 73 int tn; 74 scanf("%d",&tn); 75 int sum=0; 76 for (int j=1;j<=tn;j++) 77 { 78 int x; 79 scanf("%d",&x); 80 if (Pre[x]==0) 81 sum+=pigs[x]; else 82 Add(Pre[x],i,Inf); 83 Pre[x]=i; 84 } 85 if (sum!=0) Add(0,i,sum); 86 scanf("%d",&limit); 87 Add(i,n+1,limit); 88 } 89 int ans=0,temp; 90 while (Bfs()) 91 while (temp=Find(0,0x7fffffff)) 92 ans+=temp; 93 printf("%d ",ans); 94 return 0; 95 } 96
LJOJ 1504
题目: Codevs 3637
Sol:显然不可能形成三角形，假设每点左边的边小于右边的边，会矛盾，然后用Prim就可以了.
1 #include<cstdio> 2 #include<cmath> 3 #include<algorithm> 4 using namespace std; 5 6 const int MAXN=5010; 7 const int INF=1e8; 8 9 struct point{ 10 int x,y; 11 point(){} 12 point(int x,int y):x(x),y(y){} 13 }; 14 15 int n; 16 point p[MAXN]; 17 double ans=0,d[MAXN]; 18 bool vis[MAXN]={0}; 19 20 double dist(point a,point b){ 21 return sqrt((double)(a.x-b.x)*(double)(a.x-b.x)+(double)(a.y-b.y)*(double)(a.y-b.y)); 22 } 23 24 int main(){ 25 scanf("%d",&n); 26 for(int i=1;i<=n;i++){ 27 scanf("%d%d",&p[i].x,&p[i].y); 28 d[i]=INF; 29 } 30 d[1]=0; 31 for(int k,i=1;i<=n;i++){ 32 double mi=INF; 33 for(int j=1;j<=n;j++){ 34 if(!vis[j]&&mi>d[j]){ 35 mi=d[j]; 36 k=j; 37 } 38 } 39 if(mi==INF)break; 40 vis[k]=1; 41 ans+=mi; 42 for(int j=1;j<=n;j++){ 43 if(!vis[j]){ 44 d[j]=min(d[j],dist(p[j],p[k])); 45 } 46 } 47 } 48 printf("%.2lf ",ans); 49 return 0; 50 }
LJOJ 1505
1 题目描述 2 在这个繁忙的社会中，我们往往不再去选择最短的道路，而是选择最快的路线。开车时每条道路的限速成为最关键的问题。不幸的是，有一些限速的标志丢失了，因此你无法得知应该开多快。一种可以辩解的解决方案是，按照原来的速度行驶。你的任务是计算两地间的最快路线。 你将获得一现代化城市的道路交通信息。为了使问题简化，地图只包括路口和道路。每条道路是有向的，只连接了两条道路，并且最多只有一块限速标志，位于路的起点。两地A和B，最多只有一条道路从A连接到B。你可以假设加速能够在瞬间完成并且不会有交通堵塞等情况会影响你。当然，你的车速不能超过当前的速度限制。 3 输入 4 输入文件SPEED.IN的第一行是3个整数 N,M和D,2<=N<=150,表示道路的数目,用0..N-1标记。M是道路的总数，D表示你的目的地。接下来的M行，每行描述一条道路，每行有4个整数A (0 <= A < N), B (0 <= B < N), V (0 <= V <= 500) and L (1<= L <= 500),这条路是从A到B的，速度限制是V，长度为L。如果V是0，表示这条路的限速未知。如果V不为0，则经过该路的时间T=L/V。否则T = L / Vold，Vold是你到达该路口前的速度。开始时你位于0点，并且速度为70。 5 输出 6 输出文件SPEED.OUT仅一行整数，表示从0到D经过的城市。输出的顺序必须按照你经过这些城市的顺序，以0开始，以D结束。仅有一条最快路线。 7 样例输入 8 6 15 1 9 0 1 25 68 10 0 2 30 50 11 0 5 0 101 12 1 2 70 77 13 1 3 35 42 14 2 0 0 22 15 2 1 40 86 16 2 3 0 23 17 2 4 45 40 18 3 1 64 14 19 3 5 0 23 20 4 1 95 8 21 5 1 0 84 22 5 2 90 64 23 5 3 36 40 24 样例输出 25 0 5 2 3 1 26 提示 27 The required time is in this case 2,628 units. 28 29 来源
Sol:Bfs记录上次的速度
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 #include <vector> 6 #include <queue> 7 #define fi first 8 #define se second 9 #define Pa pair<int,int> 10 #define mp make_pair 11 #define pb push_back 12 using namespace std; 13 int n,m,T,u,v,Speed,Len,vis[200][600]; 14 double f[200][600]; 15 vector<pair<int,pair<int,int> > > E[200]; 16 vector<int> Ret; 17 Pa Pre[200][600]; 18 int main() 19 { 20 scanf("%d%d%d",&n,&m,&T); 21 for (int i=1;i<=m;i++) 22 { 23 scanf("%d%d%d%d",&u,&v,&Speed,&Len); 24 E[u].pb(mp(v,mp(Speed,Len))); 25 } 26 queue<Pa> Q; 27 memset(f,0x7f,sizeof(f)); 28 memset(vis,false,sizeof(vis)); 29 f[0][70]=0,vis[0][70]=true; 30 Q.push(mp(0,70)); 31 while(!Q.empty()) 32 { 33 Pa u=Q.front(); vis[u.fi][u.se]=0; Q.pop(); 34 for (int i=0;i<E[u.fi].size();i++) 35 { 36 int v=E[u.fi][i].fi,Race=u.se; 37 if (E[u.fi][i].se.fi!=0) Race=E[u.fi][i].se.fi; 38 if (f[u.fi][u.se]+E[u.fi][i].se.se*1.0/Race<f[v][Race]) 39 { 40 f[v][Race]=f[u.fi][u.se]+E[u.fi][i].se.se*1.0/Race; 41 Pre[v][Race]=u; 42 if (!vis[v][Race]) 43 { 44 vis[v][Race]=true; 45 Q.push(mp(v,Race)); 46 } 47 } 48 } 49 } 50 double Ans=1e30; Pa Opt; 51 for (int i=0;i<=500;i++) if (f[T][i]<Ans) Ans=f[T][i],Opt=mp(T,i); 52 while (Opt.fi!=0) Ret.pb(Opt.fi),Opt=Pre[Opt.fi][Opt.se]; 53 printf("0 "); 54 for (int i=Ret.size()-1;i>0;i--) printf("%d ",Ret[i]); 55 printf("%d ",Ret[0]); 56 }
LJOJ 1506
1 题目描述 2 你的任务是把一些人分成两组，使得： •每个人都被分到其中一组； •每个组都至少有一个人； •一组中的每个人都认识其他同组成员； •两组的成员人数尽量接近。 这个问题可能有多个解决方案，你只要输出任意一个即可，或者输出这样的分组法不存在。 3 4 输入 5 为了简单起见，所有的人都用一个整数标记，每个人号码不同，从1到N。 输入文件的第一行包括一个整数N（2≤N≤100），N就是需要分组的总人数；接下来的N行对应了这N个人，按每个人号码的升序排列，每一行给出了一串号码Aij (1≤Aij≤N，Aij≠i)，代表了第i个人所认识的人的号码，号码之间用空格隔开，并以一个“0”结束。 6 7 输出 8 如果分组方法不存在，就输出信息“No solution”（输出时无需加引号）至输出文件； 否则输出两组人数的差的绝对值 9 10 样例输入 11 5 12 2 3 5 0 13 1 4 5 3 0 14 1 2 5 0 15 1 2 3 0 16 4 3 2 1 0 17 样例输出 18 1
Sol:两个不认识的人肯定不能在一组里，那么就可以相互染色，然后可以分成两块。
1 #include <cstdio> 2 const int Inf=0x3f3f3f3f; 3 int n,Cnt,Ans; 4 bool G[110][110],Tmp[110][110],F[110][110]; 5 int vis[110],c[110],a[110],b[110],x; 6 inline int Abs(int x) {return x<0?-x:x;} 7 bool Dfs(int u,int Color) 8 { 9 vis[u]=true; c[u]=Color; 10 for (int i=1;i<=n;i++) 11 { 12 if (!G[u][i]) continue; 13 if (c[i]==Color) return false; 14 else if (c[i]==0 && !Dfs(i,-Color)) return false; 15 } 16 return true; 17 } 18 bool Work() 19 { 20 for (int i=1;i<=n;i++) 21 { 22 if (vis[i]) continue; 23 if (!Dfs(i,i)) return false; 24 Cnt++; 25 for (int j=1;j<=n;j++) 26 { 27 if (Abs(c[j])!=i) continue; 28 a[Cnt]+=(c[j]==i); 29 b[Cnt]+=(c[j]==-i); 30 } 31 } 32 return true; 33 } 34 int main() 35 { 36 scanf("%d",&n); 37 for (int i=1;i<=n;i++) 38 { 39 while (true) 40 { 41 scanf("%d",&x); if (x==0) break; 42 Tmp[i][x]=true; 43 } 44 } 45 for (int i=1;i<=n;i++) for (int j=i+1;j<=n;j++) if (!(Tmp[i][j]&&Tmp[j][i])) G[i][j]=G[j][i]=true; 46 if (!Work()) {puts("No solution"); return 0;} 47 F[0][0]=1; 48 for (int Op=1;Op<=Cnt;Op++) 49 { 50 for (int i=n;i>=0;i--) for (int j=n;j>=0;j--) 51 { 52 if (i>=a[Op] && j>=b[Op] && F[i-a[Op]][j-b[Op]]) F[i][j]=true; 53 if (i>=b[Op] && j>=a[Op] && F[i-b[Op]][j-a[Op]]) F[i][j]=true; 54 } 55 } 56 Ans=10000; 57 for (int i=0;i<=n;i++) 58 { 59 int j=n-i; 60 if (Abs(i-j)<Ans && F[i][j]) Ans=Abs(i-j); 61 } 62 printf("%d ",Ans); 63 return 0; 64 }
LJOJ 1507
题目:RQNOJ P319
Sol: 可以计算出移动之后的Ans减少了Size[u]又增加加了Tot-Size[u]，Dfs找到不能移动为止。
1 #include <cstdio> 2 #include <cstring> 3 #define LL long long 4 const LL Maxn=100100; 5 struct EDGE{LL to,next;}edge[Maxn<<2]; 6 LL head[Maxn],Sum[Maxn],Dep[Maxn],Ans,n,l,r,w[Maxn],cnt,Tot; 7 inline void Add(LL u,LL v) 8 {edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt++;} 9 void Dfs(LL u,LL fa) 10 { 11 Sum[u]=w[u]; 12 for (LL i=head[u];i!=-1;i=edge[i].next) 13 { 14 if (edge[i].to==fa) continue; 15 Dep[edge[i].to]=Dep[u]+1; 16 Dfs(edge[i].to,u); 17 Sum[u]+=Sum[edge[i].to]; 18 } 19 } 20 void Find(LL u,LL fa) 21 { 22 for (LL i=head[u];i!=-1;i=edge[i].next) 23 { 24 if (edge[i].to==fa) continue; 25 if (Tot-2*Sum[edge[i].to]<0) 26 { 27 Ans+=Tot-2*Sum[edge[i].to]; 28 Find(edge[i].to,u); 29 } 30 } 31 } 32 int main() 33 { 34 scanf("%lld",&n); memset(head,-1,sizeof(head)); 35 for (LL i=1;i<=n;i++) 36 { 37 scanf("%lld%lld%lld",&w[i],&l,&r); 38 if (l) Add(i,l),Add(l,i); 39 if (r) Add(i,r),Add(r,i); 40 Tot+=w[i]; 41 } 42 Dep[1]=1; 43 Dfs(1,0); 44 for (LL i=1;i<=n;i++) Ans+=(Dep[i]-Dep[1])*w[i]; 45 Find(1,0); 46 printf("%lld ",Ans); 47 return 0; 48 }
LJOJ 1508
题目：POJ 3694
Sol:求割点当且仅当 Low[v]>=Dfn[u]
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include <algorithm> 5 using namespace std; 6 const int Maxn=1010; 7 const int Maxm=2001000; 8 struct EDGE{int to,next;}edge[Maxm]; 9 int head[Maxn],Dfn[Maxn],Low[Maxn],Root,Son,u,v,n,cnt,Stamp; 10 bool Cut[Maxn]; 11 inline int Min(int x,int y) {return x>y?y:x;} 12 inline void Add(int u,int v) {edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt++;} 13 void Tarjan(int u) 14 { 15 Dfn[u]=Low[u]=++Stamp; 16 for (int i=head[u];i!=-1;i=edge[i].next) 17 { 18 int v=edge[i].to; 19 if (!Dfn[v]) 20 { 21 Tarjan(v); 22 if (u==Root) Son++; 23 else 24 { 25 Low[u]=Min(Low[u],Low[v]); 26 if (Dfn[u]<=Low[v]) Cut[u]=true; 27 } 28 } else 29 Low[u]=Min(Low[u],Dfn[v]); 30 } 31 } 32 int main() 33 { 34 while (scanf("%d",&n)!=EOF && n!=0) 35 { 36 memset(head,-1,sizeof(head)); 37 memset(Dfn,0,sizeof(Dfn)); 38 memset(Low,0,sizeof(Low)); 39 memset(Cut,false,sizeof(Cut)); Stamp=0; cnt=0; 40 while (scanf("%d",&u)!=EOF && u!=0) 41 while (getchar()!=' ') 42 { 43 scanf("%d",&v); 44 Add(u,v),Add(v,u); 45 } 46 Root=1,Son=0; 47 Tarjan(1); 48 int Ans=0; 49 if (Son>1) Ans=1; 50 for (int i=1;i<=n;i++) if (Cut[i]) Ans++; 51 printf("%d ",Ans); 52 } 53 return 0; 54 }
LJOJ 1509
题目:POJ 3177
Sol：Tarjan缩点之后找到桥的个数/2
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 using namespace std; 6 const int Maxn=1100; 7 const int Maxm=10010; 8 struct EDGE{int to,next;}edge[Maxm<<2]; 9 int head[Maxn],Dfn[Maxn],Low[Maxn],vis[Maxn],n,m,u,v,cnt,Top,Stack[Maxn],d[Maxn],Scc,Stamp,Belong[Maxn],Ans; 10 inline void Add(int u,int v) {edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt++;} 11 inline int Min(int x,int y) {return x>y?y:x;} 12 void Tarjan(int u,int fa) 13 { 14 Dfn[u]=Low[u]=++Stamp; Stack[++Top]=u; vis[u]=true; 15 for (int i=head[u];i!=-1;i=edge[i].next) 16 { 17 int v=edge[i].to; 18 if (v==fa) continue; 19 if (!Dfn[v] && !vis[v]) 20 { 21 Tarjan(v,u); 22 Low[u]=Min(Low[u],Low[v]); 23 } else if (vis[v]) Low[u]=Min(Low[u],Dfn[v]); 24 } 25 if (Dfn[u]==Low[u]) 26 { 27 Scc++; 28 while (Stack[Top]!=u && Top>=1) Belong[Stack[Top--]]=Scc; 29 Belong[Stack[Top--]]=Scc; 30 } 31 } 32 int main() 33 { 34 scanf("%d%d",&n,&m); 35 memset(head,-1,sizeof(head)); 36 for (int i=1;i<=m;i++) 37 { 38 scanf("%d%d",&u,&v); 39 Add(u,v),Add(v,u); 40 } 41 Tarjan(1,0); 42 for (int u=1;u<=n;u++) 43 { 44 for (int i=head[u];i!=-1;i=edge[i].next) 45 if (Belong[u]!=Belong[edge[i].to]) d[Belong[edge[i].to]]++; 46 } 47 48 for (int i=1;i<=Scc;i++) if (d[i]==1) Ans++; 49 printf("%d ",(Ans+1)>>1); 50 return 0; 51 }
LJOJ 1510
题目: BZOJ 1596
Sol: 树形DP
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 using namespace std; 6 const int Maxn=20010; 7 const int Inf=10000000; 8 struct Node{int to,next;}edge[Maxn<<2]; 9 int F[Maxn][4],cnt,n,u,v,head[Maxn]; 10 bool vis[Maxn]; 11 inline int Max(int x,int y) {return x>y?x:y;} 12 inline int Min(int x,int y) {return x>y?y:x;} 13 inline int Min3(int x,int y,int z) {return Min(x,Min(y,z));} 14 inline void Add(int u,int v) 15 {edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt++;} 16 void Dp(int u,int fa) 17 { 18 F[u][0]=1; F[u][2]=0; int Mx=-Inf; 19 for (int i=head[u];i!=-1;i=edge[i].next) 20 { 21 if (fa==edge[i].to) continue; 22 Dp(edge[i].to,u); 23 F[u][0]+=Min3(F[edge[i].to][0],F[edge[i].to][1],F[edge[i].to][2]); 24 F[u][1]+=Min(F[edge[i].to][0],F[edge[i].to][1]); 25 Mx=Max(Mx,Min(F[edge[i].to][0],F[edge[i].to][1])-F[edge[i].to][0]); 26 F[u][2]+=F[edge[i].to][1]; 27 } 28 F[u][1]-=Mx; 29 } 30 int main() 31 { 32 scanf("%d",&n); 33 memset(head,-1,sizeof(head)); 34 for (int i=1;i<n;i++) scanf("%d%d",&u,&v),Add(u,v),Add(v,u); 35 36 Dp(1,0); 37 printf("%d ",Min(F[1][0],F[1][1])); 38 return 0; 39 }
POJ 1274
二分图匹配
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 const int Maxn=210; 6 int n,m,G[Maxn][Maxn],Link[Maxn],vis[Maxn],num,tmp; 7 bool Dfs(int u) 8 { 9 for (int i=1;i<=m;i++) 10 if (G[u][i] && !vis[i]) 11 { 12 vis[i]=true; 13 if (Link[i]==-1 || Dfs(Link[i])) 14 { 15 Link[i]=u; 16 return true; 17 } 18 } 19 return false; 20 } 21 int Hungary() 22 { 23 int Ret=0; 24 memset(Link,-1,sizeof(Link)); 25 for (int i=1;i<=n;i++) 26 { 27 for (int j=1;j<=m;j++) vis[j]=false; 28 if (Dfs(i)) Ret++; 29 } 30 return Ret; 31 } 32 int main() 33 { 34 while (scanf("%d%d",&n,&m)!=EOF) 35 { 36 memset(G,false,sizeof(G)); 37 for (int i=1;i<=n;i++) 38 { 39 scanf("%d",&num); 40 for (int j=1;j<=num;j++) 41 { 42 scanf("%d",&tmp); 43 G[i][tmp]=true; 44 } 45 } 46 printf("%d ",Hungary()); 47 } 48 return 0; 49 }
BZOJ 3436
差分约束系统
a比b至少多c -> Add(a,b,-c)
a比b至多多c -> b比a至少多-c Add(b,a,c);
a比b相等 ->Add(a,b,0)
加入满足所有要求则无负环
1 #include <cstdio> 2 #include <cstring> 3 const int Inf=0x3f3f3f3f; 4 const int Maxn=10010; 5 struct EDGE{int to,next,w;}edge[Maxn<<1]; 6 int head[Maxn<<1],n,m,vis[Maxn],Type,x,y,z,Dis[Maxn],flag,cnt; 7 inline void Add(int u,int v,int w) 8 {edge[cnt].to=v;edge[cnt].next=head[u];edge[cnt].w=w;head[u]=cnt++;} 9 void Spfa(int u) 10 { 11 vis[u]=true; 12 for (int i=head[u];i!=-1;i=edge[i].next) 13 if (Dis[edge[i].to]>Dis[u]+edge[i].w) 14 { 15 if (vis[edge[i].to]) 16 { 17 flag=true; 18 break; 19 } 20 Dis[edge[i].to]=Dis[u]+edge[i].w; 21 Spfa(edge[i].to); 22 } 23 vis[u]=false; 24 } 25 26 int main() 27 { 28 29 // freopen("c.in","r",stdin); 30 scanf("%d%d",&n,&m); 31 memset(head,-1,sizeof(head)); 32 for (int i=1;i<=m;i++) 33 { 34 scanf("%d",&Type); 35 if (Type==1) 36 { 37 scanf("%d%d%d",&x,&y,&z); 38 Add(x,y,-z); 39 } 40 if (Type==2) 41 { 42 scanf("%d%d%d",&x,&y,&z); 43 Add(y,x,z); 44 } 45 if (Type==3) 46 { 47 scanf("%d%d",&x,&y); 48 Add(x,y,0); 49 } 50 } 51 for (int i=1;i<=n;i++) 52 { 53 Dis[i]=0; 54 Spfa(i); 55 } 56 if (flag) puts("No"); else puts("Yes"); 57 return 0; 58 }
BZOJ 2068
拓扑排序，把所有入度为0的节点u加入对列，那么a[u]就是可以选择。再把a[a[u]]为0的加入
剩下的还有环就是环内总的节点数/2就可以了
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <algorithm> 5 using namespace std; 6 const int Maxn=1001000; 7 int a[Maxn],In[Maxn],Q[Maxn],n,Ans; 8 bool vis[Maxn]; 9 void Work() 10 { 11 int l=1,r=0; 12 for (int i=1;i<=n;i++) if (In[i]==0) Q[++r]=i,vis[i]=true; 13 while (l<=r) 14 { 15 int u=Q[l++]; 16 if (!vis[a[u]]) 17 { 18 Ans++,vis[a[u]]=true,In[a[a[u]]]--; 19 if (!vis[a[a[u]]] && In[a[a[u]]]==0) 20 { 21 vis[a[a[u]]]=true; 22 Q[++r]=a[a[u]]; 23 } 24 } 25 } 26 for (int i=1;i<=n;i++) 27 if (!vis[i]) 28 { 29 int t=i,Ret=0; 30 while (!vis[t]) Ret++,vis[t]=true,t=a[t]; 31 Ans+=(Ret>>1); 32 } 33 } 34 int main() 35 { 36 37 scanf("%d",&n); 38 for (int i=1;i<=n;i++) 39 { 40 scanf("%d",&a[i]); 41 In[a[i]]++; 42 } 43 Work(); 44 printf("%d ",Ans); 45 return 0; 46 }
BZOJ 1734 Floyd求最小环，就是多了一层对于已经得知道最短路，然后从后面的值开始Floyd
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 5 using namespace std; 6 const int N=105; 7 const int Inf=9999999; 8 9 int G[N][N],F[N][N],u,v,w; 10 int Road[N][N],Path[N]; 11 12 int m,n,cnt,Ans; 13 14 void Record(int s,int t) 15 { 16 if (Road[s][t]) 17 { 18 Record(s,Road[s][t]); 19 Record(Road[s][t],t); 20 } 21 else Path[cnt++]=t; 22 } 23 24 void Floyd() 25 { 26 Ans=Inf; 27 for (int k=1;k<=n;k++) 28 { 29 for (int i=1;i<k;i++) 30 for (int j=i+1;j<k;j++) 31 if (Ans>F[i][j]+G[i][k]+G[k][j]) 32 { 33 Ans=F[i][j]+G[i][k]+G[k][j]; 34 cnt=0; 35 Path[cnt++]=i; 36 Record(i,j); 37 Path[cnt++]=k; 38 } 39 for (int i=1;i<=n;i++) 40 for (int j=1;j<=n;j++) 41 if (F[i][j]>F[i][k]+F[k][j]) 42 { 43 F[i][j]=F[i][k]+F[k][j]; 44 Road[i][j]=k; 45 } 46 } 47 } 48 int main() 49 { 50 while (scanf("%d%d",&n,&m)!=EOF) 51 { 52 for (int i=1;i<=n;i++) 53 for (int j=1;j<=n;j++) F[i][j]=Inf,Road[i][j]=0; 54 for (int i=1;i<=m;i++) 55 { 56 scanf("%d%d%d",&u,&v,&w); 57 if (w<F[u][v]) F[u][v]=F[v][u]=w; 58 } 59 for (int i=1;i<=n;i++) 60 for (int j=1;j<=n;j++) G[i][j]=F[i][j]; 61 Floyd(); 62 if (Ans==Inf) puts("No solution."); 63 else 64 { 65 printf("%d",Path[0]); 66 for (int i=1;i<cnt;i++) 67 printf(" %d",Path[i]); 68 puts(""); 69 } 70 } 71 return 0; 72 }