Go: package main

import (
	"sort"
	"bufio"
	"os"
	"fmt"
	"strconv"
)

// Define disjoint-set forests

type vertexNo int

type vertex struct {
	number vertexNo
	rank int
	parent *vertex
}

func makeSet(number vertexNo) vertex {
	return vertex {
		number: number,
		rank: 0,
		parent: nil,
	}
}

func union(x, y *vertex) *vertex {
	return link(findSet(x), findSet(y))
}

func link(x, y *vertex) *vertex {
	if x.rank > y.rank {
		y.parent = x
		return x
	} else {
		x.parent = y
		if x.rank == y.rank {
			y.rank += 1
		}
		return y
	}
}

func findSet(v *vertex) *vertex {
	if v.parent == nil {
		return v
	}
	v.parent = findSet(v.parent)
	return v.parent
}

// Function to calculate minimum spanning tree

type matrix [][]int

type edge struct {
	one vertexNo
	two vertexNo
	weight int
}

const(
	MAX_VERTICES = 100
)

func newEdgesFromMatrix(matrix matrix) []edge {
	edges := make([]edge, 0, MAX_VERTICES * MAX_VERTICES)

	for i := 0; i < len(matrix); i++ {
		row := matrix[i]
		for j := i; j < len(row); j++ {
			if row[j] >= 0 {
				edge := edge{one: vertexNo(i+1), two: vertexNo(j+1), weight: row[j]}
				edges = append(edges, edge)
			}
		}
	}

	return edges
}

func solve(matrix matrix) int {
	nVertices := len(matrix)

	// initialize vertices
	vertices := make([]vertex, nVertices)
	for i := 0; i < nVertices; i++ {
		vertices[i] = makeSet(vertexNo(i+1))
	}

	// prepare edges sorted in the ascending order of weight
	edges := newEdgesFromMatrix(matrix)
	sort.Slice(edges, func(i, j int) bool {
		return edges[i].weight < edges[j].weight
	})

	unionCount := 0
	weightSum := 0
	for unionCount < nVertices-1 {
		for i := 0; i < len(edges); i++ {
			edge := edges[i]
			v1 := &vertices[edge.one-1]
			v2 := &vertices[edge.two-1]
			rep1 := findSet(v1)
			rep2 := findSet(v2)
			if rep1.number == rep2.number {
				continue
			}
			union(rep1, rep2)
			edges = append(edges[0:i], edges[i+1:]...)
			weightSum += edge.weight
			unionCount++
			break
		}
	}

	return weightSum
}

// main

func main() {
	sc := bufio.NewScanner(os.Stdin)
	sc.Split(bufio.ScanWords)

	n, err := nextInt(sc)
	if err != nil {
		panic(err)
	}

	matrix := parseMatrix(n, sc)

	edges := newEdgesFromMatrix(matrix)
	sort.Slice(edges, func(i, j int) bool {
		return edges[i].weight < edges[j].weight
	})

	fmt.Println(solve(matrix))
}

func parseMatrix(n int, sc *bufio.Scanner) [][]int {
	matrix := make([][]int, n)
	for i := 0; i < n; i++ {
		matrix[i] = make([]int, n)
		for j := 0; j < n; j++ {
			x, err := nextInt(sc)
			if err != nil {
				panic(err)
			}
			matrix[i][j] = x
		}
	}
	return matrix
}

func nextText(sc *bufio.Scanner) string {
	sc.Scan()
	return sc.Text()
}

func nextInt(sc *bufio.Scanner) (int, error) {
	return strconv.Atoi(nextText(sc))
}

C++: #include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;

#define INF 20000000
typedef struct Edge edge;

struct Edge{
  int src,dst,weight;
  Edge(int start,int end,int cost) : src(start),dst(end),weight(cost) {}
};

typedef vector<edge> Edges;

bool operator < (const edge &e,const edge &f){
  return (e.weight != f.weight) ? (e.weight > f.weight) : (e.src != f.src) ? (e.src < f.src) : (e.dst < f.dst);
}

class UnionFind{
public:
  int par[1001];
  int rank[1001];

  UnionFind(){
	for(int i = 0;i < 1001;i++){
	  par[i] = i;
	  rank[i] =1;
	}
  }

  int root(int);
  bool findset(int,int);
  void unionset(int,int);
};

void kruskal(Edges);

int main(int argc,char* argv[]){
  int V;
  Edges es;
  int dist;

  cin >> V;
  for(int i = 0;i < V;i++){
	for(int j = 0;j < V;j++){
	  cin >> dist;
	  if(dist != -1){
		edge e(i,j,dist);
		es.push_back(e);
	  }
	}
  }

  kruskal(es);
  return 0;
}

void kruskal(Edges es){
  priority_queue<edge> que;
  int cost = 0;
  UnionFind uf;
  for(int i = 0;i < (int)es.size();i++)
	que.push(es[i]);

  while(!que.empty()){
	edge e = que.top();que.pop();

	if(!uf.findset(e.src,e.dst)){
	  uf.unionset(e.src,e.dst);
	  cost += e.weight;
	}
  }

  printf("%d\n",cost);
}


int UnionFind::root(int x){
  return (par[x] == x) ? x : root(par[x]);
}

bool UnionFind::findset(int x,int y){
  return root(x) == root(y);
}

void UnionFind::unionset(int x,int y){
  x = root(x);
  y = root(y);

  if(x == y)
	return;

  if(rank[y] < rank[x])
	swap(x,y);
  par[y] = x;
  rank[x] += rank[y];
}

C++: #include <iostream>
#include <vector>

using namespace std;
const int INF = 1001001001;

int prim(vector< vector<int> > cost)
{
    vector<int> V(1, 0), T(cost.size());
    for(int i=0;i<cost.size();++i) T[i] = i;
    T[0] = INF;
    int mini, minj, min_cost, ans=0;
    while(V.size() < cost.size())
    {
        min_cost = INF;
        for(int i : V){
            for(int j : T)
            {
                if(j == INF) continue;
                if(cost[i][j] < min_cost)
                {
                    mini = i;
                    minj = j;
                    min_cost = cost[i][j];
                }
            }
        }
        V.push_back(minj);
        T[minj] = INF;
        cost[mini][minj] = cost[minj][mini] = INF;
        ans += min_cost;
    }
    return ans;
}

int main()
{
    int n;
    cin >> n;
    vector<vector<int> > cost(n, vector<int>(n));

    for(int i=0;i<n;++i)
    {
        for(int j=0;j<n;++j)
        {
            cin >> cost[i][j];
            if(cost[i][j] == -1) cost[i][j] = INF;
        }
    }

    cout << prim(cost) << endl;
}

C++: #include <stdio.h>

int n,M[110][110];

int minimum(){
  int color[110],d[110],p[110];
  int i=0,s=0;
  int mincost,v,u;
  while(i<=n){color[i]=0;d[i]=5000;p[i]=-1;i++;}
  d[0]=0;
  while(1){
    mincost = 5000;
    for(i=0;i<n;i++){
      if(color[i] != 2 && d[i]<mincost){
	mincost = d[i];
	u=i;
      }
    }
    if(mincost == 5000){
      break;
    }
    color[u] = 2;
    for(v=0;v<n;v++){
      if(color[v] != 2 && M[u][v]<d[v] && M[u][v]!=-1){
	d[v] = M[u][v];
	p[v] = u;
	color[v] = 1;
      }
    }
  }
  for(i=0;i<n;i++){
    if(d[i]!=5000) s+=d[i];
  }
  return s;
}
int main(void){
  int i,j;
  scanf("%d",&n);
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      scanf("%d",&M[i][j]);
    }
  }
  printf("%d\n",minimum());
  return 0;
}

C++: #include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;

int x[10000][10000];
int MIN[10000];
int MINX = 0;
int MINZ = 0;

int main() {
	memset(MIN, 127, sizeof(MIN));
	int n, a;
	cin >> n;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			cin >> a;
			x[i][j] = a;
		}
	}
	for(int h=0;h<n;h++){
		memset(MIN, 127, sizeof(MIN));
		MIN[h] = 0;
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				for (int k = 0; k < n; k++) {
					if (x[j][k] != -1) {
						MIN[k] = min(MIN[k], x[j][k]);
					}
				}
			}
		}
		for (int i = 0; i < n; i++) {
			MINX += MIN[i];
		}
		MINZ=min(MINZ,MINX);
	}
	cout << MINZ << endl;
	return 0;
}