1 条题解

  • 0
    @ 2026-7-4 22:44:48

    #include <cstdio>
    const int M = 1205;
    const int MOD = 1e9+7;
    #define int long long
    int read()
    {
    	int x=0,f=1;char c;
    	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
    	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
    	return x*f;
    }
    int n,m,a[M],g[M][M],f[M][M],C[M][M];
    void add(int &x,int y) {x=(x+y)%MOD;}
    signed main()
    {
    	n=read();m=n<<1;
    	for(int i=1;i<=n;i++) a[read()]=1;
    	for(int i=0;i<=m;i++)
    	{
    		C[i][0]=1;
    		for(int j=1;j<=i;j++)
    			C[i][j]=(C[i-1][j-1]+C[i-1][j])%MOD;
    	}
    	g[0][0]=1;
    	for(int i=1;i<=n;i++)
    		for(int j=0;j<=i;j++)
    		{
    			g[i][j]=g[i-1][j];
    			if(j>0) add(g[i][j],2*j*g[i-1][j-1]);
    			if(j>1) add(g[i][j],j*(j-1)*g[i-1][j-2]);
    		}
    	f[m+1][0]=1;
    	for(int i=m,c0=0,c1=0;i>=1;i--)
    	{
    		if(a[i])
    		{
    			c1++;
    			for(int j=c0;j<=c1-1;j++)
    			{
    				add(f[i][j],f[i+1][j]);
    				for(int k=1;k<=c1-j;k++)
    					add(f[i][j+k],f[i+1][j]*C[c1-j-1][k-1]
    					%MOD*g[k-1][k-1]%MOD*(k+1));
    			}
    		}
    		else
    		{
    			c0++;
    			for(int j=c0;j<=c1;j++)
    				add(f[i][j],f[i+1][j]*(j-c0+1));
    		}
    	}
    	int ans=f[1][n],inv2=(MOD+1)/2;
    	for(int i=1;i<=n;i++) ans=ans*inv2%MOD;
    	printf("%lld\n",ans);
    }
    
    
    • 1

    信息

    ID
    8399
    时间
    4000ms
    内存
    512MiB
    难度
    10
    标签
    递交数
    2
    已通过
    1
    上传者