Pagini recente » Cod sursa (job #1484188) | Cod sursa (job #944363) | Cod sursa (job #1829606) | Cod sursa (job #2907998) | Cod sursa (job #2490802)
#include <fstream>
#include <vector>
#include <algorithm>
#include <deque>
#include <bitset>
using namespace std;
ifstream fin("biconex.in");
ofstream fout("biconex.out");
int n, i, t[100001], f[100001], low[100001], m, x, y, componente, j, p;
bitset<100002> pun;
vector<int> v[100001], punctedearticulatie;
deque<pair<int, int> > stiva;
deque<int> sol[100001];
void dfs(int nod, int niv){
f[nod] = niv;
/// low[nod] = cel mai mic nivel la care pot ajunge dintr-un fiu al lui nod pe alta cale decat trecand iar prin tata
/// cu siguranta, din nod pot ajunge in el insusi, deci low[nod] este macar egal cu niv
low[nod] = niv;
int fii = 0; /// numarul de componente conexe care pornesc din nod
for(int i = 0;i<v[nod].size();i++){
int vecin = v[nod][i];
if(f[vecin] == 0){
t[vecin] = nod;
fii++;
stiva.push_back(make_pair(nod, vecin));
dfs(vecin, niv+1);
low[nod] = min(low[nod], low[vecin]);
if(low[vecin] >= niv){ /// daca din unul dintre fii nu se poate ajunge mai sus de nod
if(niv != 1 && pun[nod] == 0)
punctedearticulatie.push_back(nod), pun[nod] = 1; /// atunci nod este punct de articulatie
/// daca este punct de articulatie, inseamna ca delimiteaza o componenta conexa;
/// asa ca sterg toate muchiile adaugate in stiva pana la aceasta muchie
/// ele formand componenta biconexa
componente++;
y = stiva.back().second;
x = stiva.back().first;
if(stiva.back().first == nod){
sol[componente].push_back(stiva.back().second);
sol[componente].push_back(stiva.back().first);
stiva.pop_back();
}
else{
do{
sol[componente].push_back(stiva.back().second);
stiva.pop_back();
}
while(!stiva.empty() && stiva.back().first != nod);
stiva.pop_back();
sol[componente].push_back(nod);
sol[componente].push_back(vecin);
}
}
}
else
if(t[nod] != vecin)
low[nod] = min(low[nod], f[vecin]);
}
if(fii >= 2 && nod == 1 && pun[nod] == 0)
punctedearticulatie.push_back(nod), pun[nod] = 1;
}
int main(){
fin>>n>>m;
for(i=1;i<=m;i++){
fin>>x>>y;
v[x].push_back(y);
v[y].push_back(x);
}
dfs(1, 1);
fout<<componente<<"\n";
for(i=1;i<=componente;i++){
for(j=0;j<sol[i].size();j++)
fout<<sol[i][j]<<" ";
fout<<"\n";
}
return 0;
}