#include <stdio.h>
#include <math.h>
#include <stdlib.h>

main()
{
      int s[41][41];
      int i,j,ix,jx,nx,ns,nstop;
      long iseed,m,mx,mm;
      double w,e,p,pl;
      FILE *fp;

      if ((fp=fopen("moment.data","w+"))==NULL) {
           fprintf(stderr,"cannot open moment.data\n");
           exit(-2);
      }
      scanf("%d %d %d %f %d",&ix,&jx,&nstop,&w,&iseed);
      srand48(iseed);
      mx = (ix-1)*(jx-1);
      for (i=0; i<=ix; i++) {
        for (j=0; j<=jx; j++) {
            s[i][j]=1;
        }
      }
      for (ns=0; ns < nstop; ns++) {   
          for (m=0; m<mx; m++) {
              i=(int)((ix-1)*drand48())+1;
              j=(int)((jx-1)*drand48())+1;
              e=-s[i][j]*(s[i][j-1]+s[i][j+1]+s[i-1][j]+s[i+1][j]);
              p=exp((e-4.0)/w);
	      if (drand48() < p) s[i][j]=-s[i][j];
              if ((i==ix-1)||(i==1)||(j==jx-1)||(j==1)) {
	         s[ix][j] =s[1][j];
                 s[0][j]  =s[ix-1][j];
                 s[i][jx] =s[i][1];
                 s[i][0]  =s[i][jx-1];
                 s[ix][jx]=s[1][1];
                 s[ix][0] =s[0][jx-1];
                 s[0][0]  =s[ix-1][jx-1];
                 s[0][jx] =s[ix-1][1];
              }
          }
          mm = 0;
          for (i=1; i<ix; i++) {
             for (j=1; j<jx; j++) {
                 mm = mm+s[i][j];
             }
          }
          pl = mm/(float)mx;
          fprintf(fp,"%d %f\n",ns,pl);
       }

       for (i=0; i<=ix; i++) {
           for (j=0; j<=jx; j++) {
               if (s[i][j]==1) printf("%d %d\n",i,j);
           }
       }       
       fclose(fp);
}