离散数学中的warshall算法应用

#include<stdio.h>
#include<string.h>
int a[110][110];

void warshell(int n)
{

    int k, i, j;
    for(i=1; i<=n; i++)
    {
        for(j=1; j<=n; j++)
        {
            if(a[j][i])
            {
                for(k=1; k<=n; k++)
                {
                    a[j][k] = a[j][k] + a[i][k];
                    if(a[j][k] >= 1)
                        a[j][k] = 1;
                }
            }
        }
    }
}
int main()
{
    int m, n, x, y, i, j, ans, flag[110];
    while(~scanf("%d %d", &n, &m))
    {
        ans = 0;
        memset(a, 0, sizeof(a));
        memset(flag, 1, sizeof(flag));
        while(m--)
        {
            scanf("%d %d", &x, &y);
            a[x][y] = 1;
        }
        warshell(n);
        for(j=1; j<=n; j++)
        {

            for(i=1; i<=n; i++)
            {
                if(i == j)
                {
                    if(a[i][j])
                    {
                        flag[j] = 0;
                        break;
                    }
                }
                else if(!a[i][j])
                {
                    if(!a[j][i])
                    {
                        flag[j] = 0;
                        break;
                    }
                }
            }
            if(flag[j])
            {
                ans++;
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}