BZOJ 4806 - 4809 象棋四题

Posted You Siki

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 4806 - 4809 象棋四题相关的知识,希望对你有一定的参考价值。

4806: 炮

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 103  Solved: 72
[Submit][Status][Discuss]

Description

众所周知,双炮叠叠将是中国象棋中很厉害的一招必杀技。炮吃子时必须隔一个棋子跳吃,即俗称"炮打隔子"。 
炮跟炮显然不能在一起打起来,于是rly一天借来了许多许多的炮在棋盘上摆了起来……他想知道,在N×M的矩形
方格中摆若干炮(可以不摆)使其互不吃到的情况下方案数有几种。
棋子都是相同的。

 

Input

一行,两个正整数N和M。
N<=100,M<=100

 

Output

一行,输出方案数mod 999983。

 

Sample Input

1 3

Sample Output

7

HINT

 

Source

[Submit][Status][Discuss]

 

 1 #include<cstdio>
 2 #define N 105
 3 #define P 999983
 4 int n,m;long long f[N][N][N],ans;
 5 inline long long C(int a,int b){
 6     return b==1?a:(1LL*a*(a-1)/2)%P;
 7 }
 8 main(){
 9     scanf("%d%d",&n,&m),f[0][m][0]=1;
10     for(int i=0;i<n;++i)
11         for(int j=0;j<=m;++j)
12             for(int k=0;k<=m;++k)
13                 if(f[i][j][k]){
14                     (f[i+1][j][k]+=f[i][j][k])%=P;
15                     if(j>0)(f[i+1][j-1][k+1]+=f[i][j][k]*C(j,1))%=P;
16                     if(j>1)(f[i+1][j-2][k+2]+=f[i][j][k]*C(j,2))%=P;
17                     if(j>0)(f[i+1][j-1][k]+=f[i][j][k]*C(j,1)%P*C(k,1))%=P;
18                     if(k>0)(f[i+1][j][k-1]+=f[i][j][k]*C(k,1))%=P;
19                     if(k>1)(f[i+1][j][k-2]+=f[i][j][k]*C(k,2))%=P;
20                 }
21     for(int i=0;i<=m;++i)
22         for(int j=0;j<=m;++j)
23             ans=(ans+f[n][i][j])%P;
24     printf("%lld\n",ans);
25 }

 

4807: 車

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 126  Solved: 45
[Submit][Status][Discuss]

Description

众所周知,車是中国象棋中最厉害的一子之一,它能吃到同一行或同一列中的其他棋子。車跟車显然不能在一起打
起来,于是rly一天又借来了许多许多的車在棋盘上摆了起来……他想知道,在N×M的矩形方格中摆最多个数的車
使其互不吃到的情况下方案数有几种。但是,由于上次摆炮摆得实在太累,他为了偷懒,打算增加一个条件:对于
任何一个車A,如果有其他一个車B在它的上面(車B行号小于車A),那么車A必须在車B的右边(車A列号大于車B)
棋子都是相同的。

 

Input

一行,两个正整数N和M。
N<=1000000,M<=1000000

 

Output

一行,输出方案数的末尾50位(不足则直接输出)。
 

 

Sample Input

2 2

Sample Output

1

HINT

 

Source

[Submit][Status][Discuss]

 

 1 #include<cstdio>
 2 #define N 1000005
 3 int n,m,t,pri[N],que[N],cnt[N],ans[64],tot=1,mor;
 4 void calc(int lim,int flg){
 5     for(int i=2;i<=lim;++i)
 6         for(int num=i;num!=1;)
 7             cnt[pri[num]]+=flg,num/=pri[num];
 8 }
 9 main(){
10     scanf("%d%d",&n,&m);
11     if(n<m)n^=m^=n^=m;
12     for(int i=2;i<=n;++i){
13         if(!pri[i])que[++t]=pri[i]=i;
14         for(int j=1;j<=t&&i*que[j]<=n;++j){
15             pri[que[j]*i]=que[j];
16             if(!(i%que[j]))break;
17         }
18     }
19     ans[1]=1;
20     calc(n,+1);
21     calc(m,-1);
22     calc(n-m,-1);
23     for(int i=1;i<=t;++i)
24         for(int j=1;j<=cnt[que[i]];++j){tot+=7;
25             for(int k=1;k<=tot;++k)ans[k]*=que[i];
26             for(int k=1;k<=tot;++k)ans[k+1]+=ans[k]/10,ans[k]%=10;
27             while(tot&&!ans[tot])--tot;if(tot>50)tot=50;
28         }
29     while(tot&&!ans[tot])--tot;
30     for(int i=tot;i;--i)putchar(0+ans[i]);
31 }

 

4808: 马

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 106  Solved: 43
[Submit][Status][Discuss]

Description

众所周知,马后炮是中国象棋中很厉害的一招必杀技。"马走日字"。本来,如果在要去的方向有别的棋子挡住(俗
称"蹩马腿"),则不允许走过去。为了简化问题,我们不考虑这一点。马跟马显然不能在一起打起来,于是rly在
一天再次借来了许多许多的马在棋盘上摆了起来……但这次,他实在没兴趣算方案数了,所以他只想知道在N×M的
矩形方格中摆马使其互不吃到的情况下的最多个数。但是,有一个很不幸的消息,rly由于玩得太Happy,质量本来
就不好的棋盘被rly弄坏了,不过幸好只是破了其中的一些格子(即不能再放子了),问题还是可以继续解决的。

 

Input

一行,两个正整数N和M。
接下来N行,每行M个数,要么为0,表示没坏,要么为1,表示坏了。
N<=200,M<=200

 

Output

一行,输出最多的个数。

 

Sample Input

2 3
0 1 0
0 1 0

Sample Output

2

HINT

 

Source

[Submit][Status][Discuss]
  1 #include <cstdio>
  2  
  3 template <class T>
  4 inline T min(const T &a, const T &b)
  5 {
  6     return a < b ? a : b;
  7 }
  8  
  9 inline int nextInt(void)
 10 {
 11     int r = 0;
 12     int f = 0;
 13     int c = getchar();
 14      
 15     for (; c < 48; c = getchar())
 16         if (c == -)f = 1;
 17      
 18     for (; c > 47; c = getchar())
 19         r = r * 10 + c - 48;
 20      
 21     return f ? -r : r;
 22 }
 23  
 24 const int mxn = 40005;
 25 const int mxm = 2000005;
 26 const int inf = 1000000007;
 27  
 28 int s;
 29 int t;
 30 int tt;
 31 int hd[mxn];
 32 int nt[mxm];
 33 int to[mxm];
 34 int fl[mxm];
 35  
 36 inline void add(int x, int y)
 37 {
 38     nt[tt] = hd[x], to[tt] = y, fl[tt] = 1, hd[x] = tt++;
 39     nt[tt] = hd[y], to[tt] = x, fl[tt] = 0, hd[y] = tt++;
 40 }
 41  
 42 int dep[mxn];
 43 int que[mxn];
 44 int cur[mxn];
 45  
 46 inline bool bfs(void)
 47 {
 48     for (int i = s; i <= t; ++i)
 49         dep[i] = 0;
 50      
 51     int l = 0, r = 1;
 52     dep[s] = 1;
 53     que[0] = s;
 54      
 55     while (l != r)
 56     {
 57         int u = que[l++], v;
 58          
 59         for (int i = hd[u]; ~i; i = nt[i])
 60             if (v = to[i], !dep[v] && fl[i])
 61                 dep[que[r++] = v] = dep[u] + 1;
 62     }
 63      
 64     return dep[t];
 65 }
 66  
 67 int dfs(int u, int f)
 68 {
 69     if (!f || u == t)
 70         return f;
 71      
 72     int used = 0, flow, v;
 73      
 74     for (int &i = cur[u]; ~i; i = nt[i])
 75         if (v = to[i], fl[i] && dep[v] == dep[u] + 1)
 76         {
 77             flow = dfs(v, min(f - used, fl[i]));
 78              
 79             used += flow;
 80             fl[i] -= flow;
 81             fl[i^1] += flow;
 82              
 83             if (used == f)
 84                 return used;
 85         }
 86      
 87     if (!used)
 88         dep[u] = 0;
 89      
 90     return used;
 91 }
 92  
 93 inline int flow(void)
 94 {
 95     int ret = 0, tmp;
 96      
 97     while (bfs())
 98     {
 99         for (int i = s; i <= t; ++i)
100             cur[i] = hd[i];
101          
102         while (tmp = dfs(s, inf))
103             ret += tmp;
104     }
105      
106     return ret;
107 }
108  
109 int n, m, id[205][205];
110  
111 const int mov[8][2] =
112 {
113     {+2, +1},
114     {+2, -1},
115     {-2, +1},
116     {-2, -1},
117     {+1, +2},
118     {+1, -2},
119     {-1, +2},
120     {-1, -2}
121 };
122  
123 signed main(void)
124 {
125     n = nextInt();
126     m = nextInt();
127      
128     for (int i = 1; i <= n; ++i)
129         for (int j = 1; j <= m; ++j)
130             if (!nextInt())id[i][j] = ++t;
131      
132     ++t;
133      
134     for (int i = s; i <= t; ++i)
135         hd[i] = -1;
136      
137     for (int i = 1; i <= n; ++i)
138         for (int j = 1; j <= m; ++j)
139             if (id[i][j])
140             {
141                 if ((i ^ j) & 1)
142                     add(s, id[i][j]);
143                 else
144                     add(id[i][j], t);
145             }
146      
147     for (int i = 1; i <= n; ++i)
148         for (int j = 1; j <= m; ++j)
149             if (((i ^ j) & 1) && id[i][j])
150                 for (int k = 0; k < 8; ++k)
151                 {
152                     int x = i + mov[k][0];
153                     int y = j + mov[k][1];
154                      
155                     if (x < 1 || x > n)continue;
156                     if (y < 1 || y > m)continue;
157                      
158                     if (!id[x][y])continue;
159                      
160                     add(id[i][j], id[x][y]);
161                 }
162      
163     printf("%d\n", t - 1 - flow());
164 }

 

4809: 皇后

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 123  Solved: 62
[Submit][Status][Discuss]

Description

众所不知,rly现在不会玩国际象棋。但是,作为一个OIer,rly当然做过八皇后问题。这里再啰嗦几句,皇后可以攻击到同行同列同对角线,在n*n的方格中摆n个皇后使其互不攻击到,求不同的解的数量,这就是经典的n皇后问题。现在问题推广到n皇后问题,这个问题对于你而言实在是小菜一叠。但因为上一次rly把棋盘弄破了,又拿不出新的,所以rly打算难一点点,问题就是破棋盘上的n皇后问题。他想知道……(你们懂的)。
棋子都是相同的。

 

Input

一行,一个正整数N。
接下来N行,每行N个数,要么为0,表示没坏,要么1,表示坏了。
N<=16

 

Output

一行,输出不同的解的数量。

 

Sample Input

4
1 0 1 1
1 1 1 0
0 1 1 1
1 1 0 1

Sample Output

1

HINT

 

Source

[Submit][Status][Discuss]

 

 1 #include <cstdio>
 2 
 3 inline int nextInt(void)
 4 {
 5     int r = 0;
 6     int c = getchar();
 7     
 8     for (; c < 48; c = getchar());
 9     for (; c > 47; c = getchar())
10         r = r * 10 + c - 48;
11     
12     return r;
13 }
14 
15 const int mxn = 20;
16 
17 typedef unsigned long long lnt;
18 
19 int n, ans, map[mxn][mxn];
20 
21 void search(int x, lnt a, lnt b, lnt c)
22 {
23     if (x == n)++ans;
24     else
25     {
26         lnt d = a | b | c;
27         
28         for (register int y = 0; y < n; ++y)
29             if (!map[x][y] && !((d >> y) & 1))
30                 search(x + 1,
31                     (a | (1ULL << y)),
32                     (b | (1ULL << y)) >> 1,
33                     (c | (1ULL << y)) << 1);
34     }
35 }
36 
37 signed main(void)
38 {
39     n = nextInt();
40     
41     for (int i = 0; i < n; ++i)
42         for (int j = 0; j < n; ++j)
43             map[i][j] = nextInt();
44     
45     search(0, 0, 0, 0);
46     
47     printf("%d\n", ans);
48 }

 

 

@Author: YouSiki

 

以上是关于BZOJ 4806 - 4809 象棋四题的主要内容,如果未能解决你的问题,请参考以下文章

[bzoj 4809]皇后 (dfs)

bzoj4809 -- 组合数

bzoj4806 炮

bzoj 4806: 炮dp

BZOJ 4806 : 炮

JZOJ 1667 ( bzoj 1801 ) [ AHOI 2009 ] 中国象棋 —— DP