#include <cstdio>
#include <algorithm>

#define DIM1 4100
#define DIM2 66000
#define BASE 32
using namespace std;

int N, M, X, Y, sol;
int V[DIM2], W[DIM2];
int A[DIM1][DIM1/BASE];

inline int ub(int X){
    int nr = 0;

    while(X != 0){
        nr ++;
        X = (X & (X - 1));
    }
    return nr;
}

int main(){

    freopen("triplete.in" ,"r", stdin );
    freopen("triplete.out","w", stdout);

    scanf("%d %d", &N, &M);

    for(int i = 1; i <= M; i ++){
        scanf("%d %d", &V[i], &W[i]);

        if(V[i] > W[i])
            swap(V[i], W[i]);

        A[V[i]][W[i] / BASE] += (1 << (W[i] % BASE));
    }

    for(int i = 1; i <= M; i ++){
        for(int j = 0; j <= N / BASE; j ++)
            sol += ub( A[ V[i] ][j] & A[ W[i] ][j] );
    }

    printf("%d\n", sol);

    fclose(stdin );
    fclose(stdout);

    return 0;
}