时光真疯狂, 我一路执迷于匆忙.

那么,你是想放手一搏,还是等到年华老去,心中充满遗憾,孤独地迈向黄泉路?...

【NOIP2013模拟】粉刷匠

Description

给出n个球,其中有C1个球是颜色1的,有C2个球是颜色2的,有C3个球是颜色3的……
有Ck个球是颜色k的。求相邻两个球颜色不同的排列方案。
k<=15,ci<=6,数据组数<=2000

Solution

这种题一般很难有直接的通式,可以考虑Dp。
设Fi,j表示,前i中颜色已经放完了,其中产生了j个不合法的空的方案数。(这里的空指的是每两个球中间的位置,包括开头结尾,不合法即指两边的球的颜色相同)
我们再枚举一个p和q,表示这种颜色要放到p个合法空和q个不合法空。
那么新的j呢?
因为我们放入了q个不合法空,所以j会减去q,然后,总共有C[i+1]个球,在p+q个空中各放了一个后,剩下的无论怎么放都会产生一个不合法空,所以加上C[i+1]-p-q
即:F[i,j]->F[i+1,j+C[i+1]-p-q]
然后如何转移,首先,我们要把C[i+1]个球放到p+q个空中,每个空都必须要放,所以有Cp+q1C[i+1]1种方法放,然后,在现有的j个空中选q个,乘上Cqj,设总共有sum个空,再乘上Cpsumj
然后这道题就解决了。
Ps:由于比赛时特别方,我疯狂滴优化常数,打得非常的渣渣。

Code

#include<cstdio>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define mo 1000000007
#define N 105
using namespace std;
typedef long long ll;
int n,ty,x,c[N][N],f[20][N];
int min(int x,int y) {if (x<y) return x;else return y;}
int get() {
    char ch;
    while (!isdigit(ch=getchar()));int o=ch-48;
    while (isdigit(ch=getchar())) o=o*10+ch-48;
    return o;
}
int main() {
    c[0][0]=1;
    fo(i,1,95) {
        c[i][0]=1;
        fo(j,1,i) {
            c[i][j]=c[i-1][j]+c[i-1][j-1];
            c[i][j]-=c[i][j]/mo*mo;
        }
    } 
    for(ty=get();ty;ty--) {
        n=get();f[0][0]=1;int sum=0;
        fo(i,1,n) {
            x=get();fo(j,0,sum+x+1) f[i][j]=0;
            fo(j,0,sum) 
                fo(k,0,min(j,x))
                    fo(l,0,min(x-k,sum+1-j)) {
                        if (!(k+l)) continue;
                        int ans=f[i-1][j];int y=j+x-2*k-l;
                        ans=(ll)ans*c[x-1][k+l-1]-(ll)ans*c[x-1][k+l-1]/mo*mo;
                        ans=(ll)ans*c[j][k]-(ll)ans*c[j][k]/mo*mo;
                        ans=(ll)ans*c[sum+1-j][l]-(ll)ans*c[sum+1-j][l]/mo*mo;
                        f[i][y]+=ans;f[i][y]-=f[i][y]/mo*mo;
                    }sum+=x;
        }
        printf("%d\n",f[n][0]);
    }
}
阅读更多
版权声明:既然是蒟蒻写的文,那么各位大爷就将就着看吧~跑的比西方记者慢多了233 https://blog.csdn.net/alan_cty/article/details/51540488
个人分类: 其他dp 数论
想对作者说点什么? 我来说一句

bzoj1296【SCOI2009】粉刷

DP

AaronGZK AaronGZK

2016-02-27 20:03:42

阅读数:1098

没有更多推荐了,返回首页

不良信息举报

【NOIP2013模拟】粉刷匠

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭