#include <ncbi.h>
#include "blastapp.h"

/* pcnt -- determine the "n" parameter for Poisson statistics */
void
pcnt(hlp)
	HitListPtr	hlp;
{
	HSPPtr	hp;
	unsigned	pluscnt = 0, minuscnt = 0;

	for (hp = hlp->hp; hp != NULL; hp = hp->next) {
		/* hits on opposite strands are counted separately */
		if (hp->frame >= 0)
			hp->n = ++pluscnt;
		else
			hp->n = ++minuscnt;
	}
}