Cod sursa(job #138030)

Utilizator ciprianfFarcasanu Alexandru Ciprian ciprianf Data 17 februarie 2008 19:55:51
Problema Factoriale Scor 20
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.41 kb
#include <stdio.h>
#include <math.h>
#include <string.h>
int prim[200],v[200],frecv[200],k=0,soli=1,aux[200];
int a[20000];
void mul(int A[200], int B[200])
{
      int i, j, t, C[100];
      memset(C, 0, sizeof(C));
      for (i = 1; i <= A[0]; i++)
      {
              for (t=0, j=1; j <= B[0] || t; j++, t/=10)
                      C[i+j-1]=(t+=C[i+j-1]+A[i]*B[j])%10;
              if (i + j - 2 > C[0]) C[0] = i + j - 2;
      }
      memcpy(A, C, sizeof(C));
}

int main(){
	int i,j,x,p,n,ok=1,t;
	long long unsigned co;
	FILE*f=fopen("factoriale.in","r");
	FILE*g=fopen("factoriale.out","w");
	a[1]=1;
	a[0]=1;
	for(i=2;i<=100;i++)
		if(!v[i]){
			prim[++k]=i;
			for(j=i+i;j<=100;j+=i)
				v[j]=1;
		}
	fscanf(f,"%d%d",&n,&p);
	for (i=1;i<=n;i++){
		fscanf(f,"%d",&x);
		for(j=1;prim[j]<=x && j<=k;j++){
			co=prim[j];
			while(x/co){
				frecv[prim[j]]+=x/co;
				co=co*prim[j];
			}
		}
	}
	for(i=2;i<=100;i++)
		if(frecv[i] && frecv[i]%p){
			aux[0]=1;
			aux[1]=1;
			for(t=1;t<=p-(frecv[i]%p);t++)
				for(j=1;j<=aux[0];j++){
					aux[j]*=i;
					if(aux[j]>9) { aux[j+1]*=i; aux[j+1]+=aux[j]/10; aux[j]%=10; aux[0]++; j++;}
				}
			for(j=1;j<=aux[0];j++)
				if(aux[j]>9) { aux[j+1]+=aux[j]/10; aux[j]%=10; aux[0]++;}
			while(!aux[aux[0]]) aux[0]--;
			mul(a,aux);
			for(j=1;j<=120;j++) aux[i]=0;
		}
	for(i=a[0];i>=1;i--)
		fprintf(g,"%d",a[i]);
	return 0;	
}