import java.util.Random;

public class Gillespie {

    private static Random random = new Random();

    public static double uniform() {
        return random.nextDouble();
    }

    public static int poisson(double lambda) {
        // using algorithm given by Knuth
        // see http://en.wikipedia.org/wiki/Poisson_distribution
        int k = 0;
        double p = 1.0;
        double L = Math.exp(-lambda);
        do {
            k++;
            p *= uniform();
        } while (p >= L);
        return k-1;
    }

    public static double exp(double lambda) {
        return -Math.log(1 - Math.random()) / lambda;
    }


    public static void main(String[] args) {
	int n                 = (args.length - 1)/3;           // n species
	int[] N               = new int[n];                    // population size
	double[] avgBirthRate = new double[n];                 // mean birth rate
	double[] avgDeathRate = new double[n];                 // mean death rate
	int M                 = Integer.parseInt(args[n*3]);   // number of iterations
	int totN              = 0;                             // total population

	for (int i=0;i<n;i++){
	    N[i]            = Integer.parseInt(args[3*i]);
	    avgBirthRate[i] = Double.parseDouble(args[3*i+1]);
	    avgDeathRate[i] = Double.parseDouble(args[3*i+2]);
	    totN            = totN + N[i];
	}

	/*
	System.out.println(M);	
	for (int i=0;i<n;i++)
	    System.out.println(N[i] +" "+ avgBirthRate[i] +" "+ avgDeathRate[i]);
	*/

	double[] birthRate  = new double[n];
	double[] deathRate  = new double[n];
	double[] pBirth     = new double[n];
	double[] pDeath     = new double[n];
	double lambda       = 0.0;
	double t            = 0.0;
	double p            = 0.0;
	double time         = 0.0;
	double lwb          = 0.0;
	
	for (int iter=0;iter<M && totN > 0; iter++){
	    lambda = 0.0;

	    for (int i=0;i<n;i++){
		birthRate[i] = avgBirthRate[i] * N[i];
		deathRate[i] = avgDeathRate[i] * N[i];
		lambda    = lambda + birthRate[i] + deathRate[i];
	    }

	    for (int i=0;i<n;i++){
		pBirth[i] = birthRate[i]/lambda;
		pDeath[i] = deathRate[i]/lambda;
	    }

	    t    = exp(lambda);
	    p    = Math.random();
	    lwb  = 0.0;
	    time = time + t;

	    for (int i=0;i<n;i++){
		if (lwb <= p && p < pBirth[i]+lwb){N[i]++;totN++;break;}
		lwb = lwb + pBirth[i];
		if (lwb <= p && p < pDeath[i]+lwb){N[i]--;totN--;break;}
		lwb = lwb + pDeath[i];
	    }

	    System.out.printf("%.5f",time);
	    for (int i=0;i<n;i++) System.out.printf(" %d",N[i]);
	    System.out.println(" "+ totN +" "+ iter);
	}
    }
}