#include <bits/stdc++.h> using namespace std; struct big { int l,num[10000]; void read() { string s; cin>>s; l=s.size(); for(int i=0;i<l;i++) num[l-i]=s[i]-'0'; } void print() { for(int i=l;i>=1;i--) cout<<num[i]; } void conv(string s) { l=s.size(); for(int i=0;i<l;i++) num[l-i]=s[i]-'0'; } }; big operator+(big a,big b) { big c; c.l=max(a.l,b.l); int t=0; for(int i=1;i<=c.l;i++) { if(i<=a.l) t+=a.num[i]; if(i<=b.l) t+=b.num[i]; c.num[i]=t%10; t/=10; } if(t>0) c.num[++c.l]=t; return c; } big operator*(big a,big b) { big c; c.l=a.l+b.l-1; for(int i=1;i<=a.l;i++) for(int j=1;j<=b.l;j++) c.num[i+j-1]+=a.num[i]*b.num[j]; for(int i=1;i<c.l;i++) { c.num[i+1]+=c.num[i]/10; c.num[i]%=10; } while(c.num[c.l]>=10) { c.num[c.l+1]+=c.num[c.l]/10; c.num[c.l]%=10; c.l++; } return c; } big operator*(big a,int b) { big c; c.l=a.l; int t=0; for(int i=1;i<=c.l;i++) { t+=a.num[i]*b; c.num[i]=t%10; t/=10; } while(t>0) { c[++c.l]=t%10; t/=10; } return c; } int main() { big a; a.conv("1"); int n; cin>>n; for(int i=2;i<=n;i++) a=a*i; cout<<n<<"!="; a.print(); return 0; }