码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ 2302: [HAOI2011]Problem c( dp )

时间:2016-01-10 22:30:14      阅读:330      评论:0      收藏:0      [点我收藏+]

标签:

技术分享

dp(i, j)表示从i~N中为j个人选定的方案数, 状态转移就考虑选多少人为i编号, 然后从i+1的方案数算过来就可以了. 时间复杂度O(TN^2) 

---------------------------------------------------------------------

#include<cstdio>
#include<algorithm>
#include<cstring>
 
using namespace std;
 
typedef long long ll;
#define c(i) (cnt[i] - cnt[i + 1])
 
const int maxn = 309;
const int maxT = 11;
 
int dp[maxn][maxn], N, cnt[maxn];
int V[maxn][maxn], Comb[maxT][maxn][maxn];
int n[maxT], m[maxT], P[maxT], q[maxT][maxn], T, _T;
 
inline void upd(int &x, int t) {
if((x += t) >= P[_T])
x -= P[_T];
}
 
void Init() {
for(_T = 0; _T < T; _T++) {
Comb[_T][0][0] = 1;
for(int i = 1; i <= n[_T]; i++) {
Comb[_T][i][0] = 1;
for(int j = 1; j <= i; j++) {
Comb[_T][i][j] = Comb[_T][i - 1][j - 1] + Comb[_T][i - 1][j];
if(Comb[_T][i][j] >= P[_T])
Comb[_T][i][j] -= P[_T];
}
}
}
}
 
void Read() {
scanf("%d%d%d", n + _T, m + _T, P + _T);
for(int i = 0, p; i < m[_T]; i++)
scanf("%d%d", &p, q[_T] + i);
}
 
inline int C(int n, int m) {
if(n < m)
return 0;
return Comb[_T][n][m];
}
 
void Work() {
N = n[_T];
memset(cnt, 0, sizeof cnt);
for(int i = 0; i < m[_T]; i++)
cnt[q[_T][i]]++;
for(int i = N; i; i--)
if((cnt[i] += cnt[i + 1]) > N - i + 1) {
puts("NO");
return;
}
memset(dp, 0, sizeof dp);
if(cnt[N]) {
dp[N][0] = 0;
dp[N][1] = 1;
} else {
dp[N][0] = 1;
dp[N][1] = N - cnt[1];
}
for(int i = N; --i; ) {
for(int j = N - i + 1; j >= cnt[i]; j--)
for(int k = c(i); k <= j; k++)
upd(dp[i][j], ll(dp[i + 1][j - k]) * C(N - j + k - cnt[1] + cnt[i + 1], k - c(i)) % P[_T]);
}
printf("YES %d\n", dp[1][N]);
}
 
int main() {
scanf("%d", &T);
for(_T = 0; _T < T; _T++) Read();
Init();
for(_T = 0; _T < T; _T++) Work();
return 0;
}

--------------------------------------------------------------------- 

 

2302: [HAOI2011]Problem c

Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 407  Solved: 221
[Submit][Status][Discuss]

 

Description

n个人安排座位,先给每个人一个1~n的编号,设第i个人的编号为ai(不同人的编号可以相同),接着从第一个人开始,大家依次入座,第i个人来了以后尝试坐到ai,如果ai被占据了,就尝试ai+1ai+1也被占据了的话就尝试ai+2,……,如果一直尝试到第n个都不行,该安排方案就不合法。然而有m个人的编号已经确定(他们或许贿赂了你的上司...),你只能安排剩下的人的编号,求有多少种合法的安排方案。由于答案可能很大,只需输出其除以M后的余数即可。

 

Input

第一行一个整数T,表示数据组数

对于每组数据,第一行有三个整数,分别表示nmM

m不为0,则接下来一行有m对整数,p1q1p2q2 ,…, pmqm,其中第i对整数piqi表示第pi个人的编号必须为qi

 

Output

对于每组数据输出一行,若是有解则输出YES,后跟一个整数表示方案数mod M,注意,YES和数之间只有一个空格,否则输出NO

 

Sample Input

2

4 3 10

1 2 2 1 3 1

10 3 8882

7 9 2 9 5 10


Sample Output



YES 4

NO



HINT



100%的数据满足:1≤T≤10,1≤n≤300,0≤m≤n,2≤M≤109,1≤pi、qi≤n   且保证pi互不相同。



 

Source

BZOJ 2302: [HAOI2011]Problem c( dp )

标签:

原文地址:http://www.cnblogs.com/JSZX11556/p/5119141.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!