唉我还是too naive

  这是个整数划分题……

  我想的DP方式是f[i][j][k]表示前 i 个数拼出 j 用了 k 个数的方案数……

  转移当然是比较直观……

  但是只能得30分QAQ

  正确的DP姿势:http://blog.csdn.net/Vmurder/article/details/42551603

分析:

数据范围不大,我们可以写整数划分。

f[i][j]表示将i划分成j个互不相同的正整数,且最大不超过n 的划分方案数。


这里说一下这道题的整数划分。


我们不妨先来反向思考一下。---------------------------------------------------------------------

首先考虑f[i][j](下图每一列都代表一个数,高度就是数值)

【BZOJ】【3612】【HEOI 2014】平衡

它可以在底下添加一行,进行转移

分为两种情况:

Ⅰ. 转移过后最小数不为1

【BZOJ】【3612】【HEOI 2014】平衡

Ⅱ. 转移过后最小数为1

【BZOJ】【3612】【HEOI 2014】平衡

然后f[i][j]可以向这两个方向转移。

而我们还需要保证最大数不能大于n,那么如下图

【BZOJ】【3612】【HEOI 2014】平衡

在最后加入一层,使得当前所有整数都+1时,发现出现了一个101,而不妨当作n是100

那么显然我们可以很方便地清掉这个数带来的情况。

即当前是f[i][j],那么最后一列是大于n的情况显然只能是有一个整数n+1

不看这个n+1,情况数是f[i-(n+1)][j-1] ,我们把这个情况集删掉就好了。



正向考虑:-----------------------------------------------------------------------------------

首先不妨把刚才的图片按顺序记作图1、2、3、4。

我们把f[i][j](图1)这么多方案分成两种情况:



Ⅰ. 最小的数不为1:

好说。 直接由f[i-j][j]在底下加一行得到。就是图2。

此时原来划分出来的整数不同,新的这些整数显然依然不同。

 

Ⅱ. 最小的数为1:

那么显然它可以由f[i-j][j-1]转移得到,

即在f[i-j][j-1]代表图形下面整体+1,最后加上一个整数1,即图3。

注意此时f[i-j][j-1]代表的所有图形整数都不同(性质/定义),那么新加1后所有整数依然不同,且均>=2

这个时候再来个整数1,依然满足所有整数不同。

 

而这两种情况显然互补,即这两种情况的转移包含了f[i][j]的所有情况(两个命题“最小数是1”,“最小数不是1”,显然包含全部情况),也就是说转移完成。


但是我们注意到还需要让最大数不能超过n,

所以有了图4。

也就是我们要减去最大数超过n的情况,方法前文图下有说明。


这道题难点解决了。

现在说一下其它细节:


f[i][j]算出来后直接暴力枚举两边的权值,及用点个数(不要忘了中心支点)

然后check。end。

 1 /**************************************************************
 2     Problem: 3612
 3     User: Tunix
 4     Language: C++
 5     Result: Accepted
 6     Time:3260 ms
 7     Memory:15336 kb
 8 ****************************************************************/
 9  
10 //Huce #2 B 
11 #include<vector>
12 #include<cstdio>
13 #include<cstdlib>
14 #include<cstring>
15 #include<iostream>
16 #include<algorithm>
17 #define rep(i,n) for(int i=0;i<n;++i)
18 #define F(i,j,n) for(int i=j;i<=n;++i)
19 #define D(i,j,n) for(int i=j;i>=n;--i)
20 using namespace std;
21  
22 int getint(){
23     int v=0,sign=1; char ch=getchar();
24     while(ch<'0'||ch>'9') {if (ch=='-') sign=-1; ch=getchar();}
25     while(ch>='0'&&ch<='9') {v=v*10+ch-'0'; ch=getchar();}
26     return v*sign;
27 }
28 typedef long long LL;
29 const int N=100010,INF=~0u>>2;
30 /*******************tamplate********************/
31 LL f[N][12];
32 int n,k,P,sum[12],mx[N][12];
33  
34 int main(){
35 #ifndef ONLINE_JUDGE
36     freopen("B.in","r",stdin);
37     freopen("B.out","w",stdout);
38 #endif
39     int T=getint();
40     f[0][0]=1;
41     while(T--){
42         n=getint(); k=getint(); P=getint();
43         if (k==1){puts("1");continue;}
44         int w=n*(k-1);
45         F(i,1,w)
46             F(j,1,k-1){
47                 f[i][j]=i>=j ? (f[i-j][j]+f[i-j][j-1])%P : 0;
48                 f[i][j]=i>=n+1 ? (f[i][j]-f[i-(n+1)][j-1]+P)%P : f[i][j];
49             }
50         LL ans=0;
51         F(i,1,w) F(j,1,k-1)
52             (ans+=f[i][j]*f[i][k-j]+f[i][j]*f[i][k-j-1])%=P;
53         printf("%lld\n",ans);
54     }
55     return 0;
56 }
View Code

相关文章: