Pagini recente » Statistici Mocioi Alexandru (pletutza) | Istoria paginii problema/photo | Rating Pop Cristian (Cristi3956) | Cod sursa (job #731071) | Cod sursa (job #134056)
Cod sursa(job #134056)
#include <stdio.h>
#include <vector>
#include <string.h>
#include <queue>
using namespace std;
#define pb push_back
#define sz size()
#define NMAX 410
#define INFI 0x3f3f3f3f
int n;
int cost[NMAX][NMAX];
short cap[NMAX][NMAX], aux[NMAX][NMAX];
int sink, source, until;
int _min;
vector<int> g[NMAX];
void read()
{
scanf("%d", &n);
until = 2*n;
int i, j;
for(i = 1; i <= n; ++i)
for(j = 1; j <= n; ++j)
scanf("%d", &cost[i][n+j]), cost[n+j][i] = -cost[i][n+j], cap[i][n+j] = 1;
source = until+3;
sink = until+4;
for(i = 1; i <= n; ++i)
cap[source][i] = 1, cap[i+n][sink] = 1;
}
int bf()
{
queue<int> q;
int i, x, t[NMAX], dist[NMAX];
memset(t, 0, sizeof(t));
memset(dist, INFI, sizeof(dist));
q.push(source);
dist[source] = 0;
while(!q.empty())
{
x = q.front(), q.pop();
for(i = 1; i <= sink; ++i)
if(cap[x][i] && dist[i] > dist[x] + cost[x][i])
{
dist[i] = dist[x] + cost[x][i];
t[i] = x;
q.push(i);
}
}
if(dist[sink] == 1061109567)
return -1;
for(i = sink; i != source; i = t[i])
--cap[t[i]][i], ++cap[i][t[i]];
return dist[sink];
}
int flux()
{
int aux, c = 0;
while(1)
{
aux = bf();
//printf("%d\n", aux);
if(aux == -1) break;
}
int i, j;
for(i = 1; i <= n; ++i)
for(j = 1; j <= n; ++j)
c += cost[i][j+n] * (!cap[i][j+n] && cap[j+n][i]);//, printf("%d %d %d\n", i, j+n, cap[i][j+n]);
return c;
}
int mate[NMAX];
int form(int x, int y)
{
int i, j;
memset(cap, 0, sizeof(cap));
for(i = 1; i <= n; ++i)
if(i != x)
cap[source][i] = 1;
for(i = 1; i <= n; ++i)
if(i+n != y)
cap[i+n][sink] = 1;
for(i = 1; i <= n; ++i)
if(i != x)
for(j = 1; j <= n; ++j)
if(j+n != y)
cap[i][j+n] = 1;
}
int main()
{
freopen("paznici2.in", "r", stdin);
freopen("paznici2.out", "w", stdout);
read();
_min = flux();
memcpy(aux, cap, sizeof(cap));
printf("%d\n", _min);
int i, j;
int aux1, aux2, z;
for(i = 1; i <= n; ++i)
for(j = 1; j <= n; ++j)
if(cap[j+n][i])
{
mate[i] = j+n;
mate[j+n] = i;
break;
}
for(j = 1; j <= n; ++j)
{
for(i = 1; i <= n; ++i)
{
if(cap[j+n][i] && !cap[i][j+n])
g[j].pb(i);
else
{
form(i, j+n);
if(_min == flux() + cost[i][j+n])
{
g[j].pb(i);
}
}
}
}
vector<int> :: iterator it;
for(i = 1; i <= n; ++i)
{
printf("%d ", g[i].sz);
for(it = g[i].begin(); it != g[i].end(); ++it) printf("%d ", *it);
printf("\n");
}
return 0;
}