UVa 125 - Numbering Paths

//\\__ hr1212 __//\\

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef pair<int,int> pii;
typedef map<int,int> mi;

#define si(a) scanf("%d",&a)
#define sii(a,b) scanf("%d %d",&a,&b)
#define siii(a,b,c) scanf("%d %d %d",&a,&b,&c)
#define pi(a) printf("%d\n",a)
#define nl printf("\n");
#define pb push_back
#define mp make_pair
#define all(c) (c).begin(),(c).end()
#define f(i,a,b) for(i=a;i<b;i++)
#define rf(i,a,b) for(i=a;i>=b;i--)
#define clr(x,a) memset(x,a,sizeof(x))
#define MAX 1000100
#define MOD 1000000007

int n,m,graph[50][50];

int main(){
    int r,k,i,c=0,x=0,y=0,j,t,l,z,x1=0,y1=0;
    ll ans=0;string p;

    t=-1;
    while(si(m)!=EOF){
        t++;
        clr(graph,0);
        n=0;
        f(i,0,m){
            sii(x,y);
            graph[x][y]=1;
            n=max(n,max(x,y));
        }

        n++;
        f(k,0,n){
            f(i,0,n)
                f(j,0,n)
                    if(graph[i][k] && graph[k][j])
                        graph[i][j]+=graph[i][k]*graph[k][j];
        }
        f(k,0,n){
            if(graph[k][k]){
                f(i,0,n){
                    f(j,0,n)
                        if(graph[i][k] && graph[k][j])
                            graph[i][j]=-1;
                }
            }
        }
        printf("matrix for city %d\n",t);
        f(i,0,n){
            f(j,0,n)
                if(j!=n-1)
                    printf("%d ",graph[i][j]);
                else
                    printf("%d",graph[i][j]);
            nl;
        }
    }

    return 0;
}