Pagini recente » Cod sursa (job #2725690) | Cod sursa (job #659481) | Cod sursa (job #3188680) | Cod sursa (job #2441495) | Cod sursa (job #80581)
Cod sursa(job #80581)
# include <stdio.h>
const long int NRBITI=30;
const long int MAXN=4096/NRBITI;
const long int MAXX=4096;
const long int MAXM=100000;
long int a[MAXX+4][MAXN+4],n,m;
long long int sol;
long int nbiti[MAXM+1];
void init(){
for (long int i=1;i<=(long int)1<<16;i++) nbiti[i]=1+nbiti[i&(i-1)];
}
long int nrbiti(long int a){return nbiti[a>>15]+nbiti[a%((long int)1<<16)];}
void add(long int loc,long int inf)
{
long int adaos=0;
if (inf%NRBITI==0) adaos=1;
a[loc][inf/NRBITI+1-adaos]|=(long int)1<<(inf%NRBITI);
}
void citire()
{
long int aa,bb;
FILE *f=fopen("triplete.in","r");
fscanf(f,"%ld%ld",&n,&m);
long int i;
for (i=1;i<=m;i++)
{
fscanf(f,"%ld%ld",&aa,&bb);
add(aa,bb);
add(bb,aa);
}
fclose(f);
}
void calculeaza()
{
FILE *f=fopen("triplete.in","r");
fscanf(f,"%ld%ld",&n,&m);
long int i,aa,bb,j;
for (i=1;i<=m;i++)
{
fscanf(f,"%ld%ld",&aa,&bb);
for (j=1;j<=n/NRBITI+1;j++)
sol+=nrbiti(a[aa][j]&a[bb][j]);
}
FILE *g=fopen("triplete.out","w");
fprintf(g,"%lld\n",sol/3);
fcloseall();
}
int main()
{
citire();
calculeaza();
return 0;
}