#include <pam.h>



#include <stdio.h>

#include <stdlib.h>







int main(int argc, char **argv) {

        FILE * fin = stdin;



        struct pam inpam;

        tuple * tuplerow;

        unsigned int row;



        unsigned int i;

        unsigned int * csc = (unsigned int *) malloc(inpam.depth*sizeof(unsigned int));





        pm_init(argv[0], 0); // *** ???

        pnm_readpaminit(fin, &inpam, PAM_STRUCT_SIZE(tuple_type));

        tuplerow = pnm_allocpamrow(&inpam);



        for(row=0; row<inpam.height; row++) {

                unsigned int column;

                pnm_readpamrow(&inpam, tuplerow);

                for(column=0; column<inpam.width; ++column) {

                        unsigned int plane;

                        unsigned int sum = 0;

                        for(plane=0; plane<inpam.depth; ++plane) sum += tuplerow[column][plane];

                        for(plane=0; plane<inpam.depth; ++plane) {

                                printf("%d/%d\t",tuplerow[column][plane],sum);

                                unsigned int val = tuplerow[column][plane];

                                if(4*val > sum && 2*val < sum) {

                                } else  csc[plane]++;

                        } printf("\n");

                } printf("\n");

        }



        pnm_freepamrow(tuplerow);



        for(i=0; i<inpam.depth; i++)

                printf("%d\t%d\n", i,csc[i]);



        return EXIT_SUCCESS;

}