UVA 160

/***
Md. Namzul Hasan
Shahjalal University of Science & Technology,Sylhet.
hasan08sust@gmail.com
***/
#include<iostream>
#include<cstdio>
#include<stack>
#include<queue>
#include<map>
#include<vector>
#include<algorithm>
#include<cstring>
#include<sstream>
#include<cmath>
#include <iomanip>
using namespace std ;
typedef long long ll ;
typedef int in ;
typedef unsigned long long ull ;
const double pi = 2*acos(0) ;
#define maxi 40000
#define pf printf
#define sc scanf
#define pb push_back
#define MEM(x,y) (memset((x),(y),sizeof(x)))
#define MIN(x,y) ((x) < (y) ? (x) : (y))
#define MAX(x,y) ((x) > (y) ? (x) : (y))
#define load(array,size)  for(int i=0 ; i<size ; i++) cin>>array[i]  ;
#define new_line  pf("\n")
#define clear_data(array) memset(array,0,sizeof(array))
#define highest_int 2147483647

int isPrime(int n){
   if(n==2){
    return 1;
   }
   if(n%2==0){
    return 0;
   }

   for(int i=3 ; i<sqrt(n) ; i+=2){
    if(n%i==0){
        return 0;
    }
   }

   return 1;
}


int main()
{
   int n,d,p;
   int data[10001];
   vector<int>index;

   while(sc("%d",&n)){
        memset(data,0,sizeof(data));
    if(n==0)
        return 0;
    for(int i=n ; i>=2 ; i--){
        d = i;
        p=2 ;
        while(1){
           if(d%p==0){
            if(isPrime(p)){

                if(data[p] ==0){
                    index.pb(p);
                }
                data[p]++;
            }
            d = d/p;
            p=2;
           }else{
             p++;
           }
            if(d==1)
                break;
        }
    }

    sort(index.begin(),index.end());

    cout << setw(3) << n << "! =";
    for(int i=0 ;i<index.size();i++){
         if (i % 15 == 0 && i != 0)
             cout << "\n      ";
             cout << setw(3) << data[index[i]];
    }
    pf("\n");

    index.clear();
   }

    return 0;
}