【2015.9.5】NOIP2012模拟赛题解by唐翔昊

T1 首先筛法暴力出1到30w中具有公因数最多的数量（算出只有180），然后我们惊喜地发现我们可以直接做，最小表示法枚举第k位，然后向前判断是否合法，处理处最小的能放进去的就可以了，同时判断不可能是看这一位能否放进去，然而数据中并没有impossible

T2 这题暴力的时间复杂度为O（n^2*2^k）,理论上极限数据即使给6s都跑不过,但是我们发现这是随机数据,于是我们可以水过而且很快【思考熊，具体就是状压二进制表示选择了哪些然后暴力验证

T3 原型POJ2323，递推，但是我们要写高精度，公式见代码

http://pan.baidu.com/s/1eQm8gIQ 内有题目和题解以及代码，以下是个人代码

T1

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#define N 300005
#define M 180
using namespace std;
int p[M+5],n,f=0;
char a[N];
bool v[27];

int main()
{
    freopen("arilang.in","r",stdin);
    freopen("arilang.out","w",stdout);
    scanf("%d\n",&n);memset(p,0,sizeof(p));
    for (int i=1;i<=n/2+1;i++)
        if (n%i==0) p[++p[0]]=i;
    for (int i=1;i<=n;i++)
    {
        memset(v,1,sizeof(v));
        for (int j=1;j<=p[0];j++)
        {
            if (p[j]<i)
            {
                         //printf("%d ",a[i-p[j]]-'a'+1);
                         if (i+p[j]>n)
                         {
                                      v[a[i+p[j]-n]-'a'+1]=0;
                                      }
                         v[a[i-p[j]]-'a'+1]=0;
                         }
            else break;
            }
        char l='a';int j;
        //for (int j=1;j<=26;j++) cout<<v[j];cout<<' '<<j<<endl;
        for (j=1;j<=26;j++)
            if (v[j]!=0) break;
            else ++l;
        a[i]=l;
        if (j==27) {f=1;break;}
        }
    if (f==1) printf("Impossible");
    else
    for (int i=1;i<=n;i++)
        printf("%c",a[i]);
    return 0;
}

T2

#include <iostream>
#include <cstdio>
#include <cstring>
#define N 35
#define M 25
#define inf 1000000007
using namespace std;
int g[N][N];
int v[N],q[40];
int ans,anss,K,n,m;

int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

bool bfs(int x)
{
     memset(v,0,sizeof(v));
     int h=0,t=1;q[t]=0;
     while (h<t)
     {
           int p=q[++h];
           for (int i=1;i<n;i++)
               if (!v[i]&&(g[p][i]&x))
               {
                                       v[i]=1;
                                       q[++t]=i;
                                       }
               }
     return v[1];
     }

int main()
{
    freopen("license.in","r",stdin);
	freopen("license.out","w",stdout);
    K=read(),n=read(),m=read();ans=inf;
    for (int i=1;i<=m;i++)
    {
        int x=read(),y=read(),w=read();
        g[x][y]=g[y][x]|=1<<w;
        }
    int x;
    for (int i=1;i<1<<K;i++)
    {
        x=0;
        for (int j=i;j;j>>=1)
        {
            x+=j&1;
            }
        if (x<ans&&bfs(i))
        {
                          ans=x;
                          anss=i;
                          }
        }
    printf("%d\n",ans);
    int k=0;
    if (ans!=inf)
    {
                 for (;anss;anss>>=1,k++)
                     if (anss&1) printf("%d ",k);
                     }
    else printf("Impossible"); 
    return 0;
}

T3

#include<iostream>  
#include<cstdio>
#include<string.h>
using namespace std;  
const int MAXT=100+9,MAXN=50+9,MAXM=1500+9,E=100000000;
int n,kk;
struct bign
{
	int a[9],l;
	void operator = (int r)
	{
		for(l=0;r;r/=E)
			a[l++]=r%E;
	}
	void operator += (const bign &r)
	{
		int i;
		for(i=0;i<r.l;++i)
			a[i]+=r.a[i];
		if(r.l>l)
			l=r.l;
		for(i=0;i<l;++i)
			if(a[i]>=E)
			{
				a[i]-=E;
				++a[i+1];
			}
		if(a[l])
			++l;
	}
}dp[MAXN][MAXM];
void print(bign x)
	{
		if (x.l)
		{
			printf("%d",x.a[x.l-1]);
			for(int i=x.l-2;i>=0;--i)
				printf("%08d",x.a[i]);
			putchar(10);
		}
		else puts("0");
	}
int main()  
{  
    freopen("inversion.in","r",stdin);
    freopen("inversion.out","w",stdout);
    /*dp[1][0]=1;  
    for(int i=1;i<51;i++)  
        for(int j=0;dp[i][j];j++)  
            for(int k=0;k<=i;k++)  
                dp[i+1][j+k]+=dp[i][j];*/
    dp[0][0]=1;dp[1][0]=1;dp[2][0]=1;dp[2][1]=1;
    for (int i=3;i<=50;i++)
    {
        dp[i][0]=1;
        for (int j=1;j<=i*(i-1)/2;j++)
        {
            dp[i][j]=0;
            for (int k=1;k<=i;k++)
                dp[i][j]+=dp[i-1][j-k+1];
            }
        }  
    for(int t=1;scanf("%d%d",&n,&kk),n;++t)
	{
		print(dp[n][kk]);
	}
}