[hdu6314]Matrix

最新推荐文章于 2021-11-17 16:57:55 发布

StaroForgin

最新推荐文章于 2021-11-17 16:57:55 发布

阅读量139

点赞数 1

分类专栏： ------数学------ # 组合

本文链接：https://blog.csdn.net/Tan_tan_tann/article/details/108433751

版权

------数学------ 同时被 2 个专栏收录

76 篇文章 0 订阅

订阅专栏

组合

42 篇文章 0 订阅

订阅专栏

Matrix

题解

~~本来打着很快的，但是卡常硬是把这道题给卡了一个下午，于是就叫它“贞卡常”了~~

先令 $f_{x,y}$ 表示将一个 $x\cdot y$ 大小的矩阵，每个格子都涂成黑白，没有任何行与列全为黑的方案数。

根据容斥，枚举其中全为黑的方格数，可以得到 $f_{x,y}= \sum_{i=0}^{x}\sum_{j=0}^{y}(-1)^{i+j}C_{x}^{i}C_{y}^{j}2^{(n-i)(m-j)}$ 。

于是答案便为 $ans=\sum_{i=a}^{n}\sum_{j=b}^{m}C_{n}^{i}C_{m}^{j}f_{n-i,m-j}$

将 $f$ 的表达式带进去，可以得到

$ans=\sum_{i=a}^{n}\sum_{j=b}^{m}C_{n}^{i}C_{m}^{j}\sum_{p=0}^{n-i}\sum_{q=0}^{m-j}(-1)^{p+q}C_{n-i}^{p}C_{m-j}^{q}2^{(n-i-p)(m-j-q)}$

$= \sum_{i=a}^{n}\sum_{j=b}^{m}\sum_{p=0}^{n-i}\sum_{q=0}^{m-j}((-1)^pC_{n}^{i}C_{n-i}^{p})((-1)^{q}C_{m}^{j}C_{m-j}^{q})2^{(n-i-p)(m-j-q)}$

设 $t1=i+p,t2=j+q$ ，有

$ans= \sum_{t1=a}^{n}\sum_{t2=b}^{m}2^{(n-t1)(m-t2)}\sum_{i=a}^{t1}\sum_{j=b}^{t2}((-1)^{t1-i}C_{n}^{i}C_{n-i}^{t1-i})((-1)^{t2-q}C_{m}^{j}C_{m-j}^{t2-j})$

$=\sum_{t1=a}^{n}\sum_{t2=b}^{m}2^{(n-t1)(m-t2)}(\sum_{i=a}^{t1}(-1)^{t1-i}C_{n}^{i}C_{n-i}^{t1-i})(\sum_{j=b}^{t2}(-1)^{t2-q}C_{m}^{j}C_{m-j}^{t2-j})$

记 $fa_{t}=\sum_{i=a}^{t}(-1)^{t-i}C_{n}^{i}C_{n-i}^{t-i}$ ， $fb_{t}= \sum_{i=b}^{t}(-1)^{t-i}C_{m}^{i}C_{m-i}^{t-i}$ ，这两者可以分别用 $O\left(n^2 \right )$ 和 $O\left(m^2 \right )$ 处理出来，

带入原式得 $ans=\sum_{i=a}^{n}\sum_{j=b}^{m}fa_{i}fb_{j}$ 。时间复杂度 $O\left(n^2+m^2+nm \right )$ 。

然后你就惊奇的发现你T了。

卡常卡半天还是过不了。

其实 $fa$ 与 $fb$ 的式子可以进一步简化得 $fa_{x}=-C_{n}^{x}\sum_{i=a}^{n} \frac{C_{i-1}^{j-1}fa_{i}}{C_{n}^i}$ ， $fb$ 同理

我们将那个 $C_{n}^{x}$ 提出来，得到 $fa_x=\sum_{i=a}^{n}C_{i-1}^{j-1}fa_{i}$ ，最后再将 $C_{n}^{x}$ 乘进去，就极大的减小了常数。

得到 $ans=\sum_{i=a}^{n}\sum_{j=b}^{m}C_{n}^{i}C_{m}^{j}fa_{i}fb_{j}$

但你发现你还是得卡半天的常才过得去。

源码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
using namespace std;
#define MAXN 3005
typedef long long LL;
typedef pair<int,int> pii;
const double eps=1e-9;
const int mo=998244353;
const int N=3000;
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
    _T f=1;x=0;char s=getchar();
    while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
    while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
    x*=f;
} 
int n,m,a,b,pow2[MAXN*MAXN];
int fa[MAXN],fb[MAXN],ans,C[MAXN][MAXN];
int add(const int x,const int y){return x+y>=mo?x+y-mo:x+y;}
void init(){
    pow2[0]=1;for(int i=1;i<N*N;i++)pow2[i]=add(pow2[i-1],pow2[i-1]);
    for(int i=0;i<=N;i++)C[i][0]=C[i][i]=1;
    for(int i=1;i<=N;i++)
        for(int j=1;j<i;j++)
        C[i][j]=add(C[i-1][j-1],C[i-1][j]);
}
signed main(){
	//freopen("test.in","r",stdin);
    init();
    while(~scanf("%d %d %d %d",&n,&m,&a,&b)){
        fa[a]=fb[b]=1;ans=0;
        for(int i=a+1;i<=n;i++)
            for(int j=a;j<i;j++)
                fa[i]=(fa[i]-1ll*C[i-1][j-1]*fa[j])%mo;
        for(int i=b+1;i<=m;i++)
            for(int j=b;j<i;j++)
                fb[i]=(fb[i]-1ll*C[i-1][j-1]*fb[j])%mo;
        for(int i=a;i<=n;i++)fa[i]=1ll*fa[i]*C[n][i]%mo;
        for(int i=b;i<=m;i++)fb[i]=1ll*fb[i]*C[m][i]%mo;
        for(int i=a;i<=n;i++) 
            for(int j=b;j<=m;j++)
                ans=(ans+1ll*pow2[(n-i)*(m-j)]*fa[i]%mo*fb[j])%mo;
        if(ans<0)ans+=mo;printf("%d\n",ans);
        for(int i=a;i<=n;i++)fa[i]=0;
        for(int i=b;i<=m;i++)fb[i]=0;
    }
    return 0;
}

谢谢！！！

StaroForgin

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
[hdu6314]Matrix

Matrix题解本来打着很快的，但是卡常硬是把这道题给卡了一个下午，于是就叫它“贞卡常”了先令表示将一个大小的矩阵，每个格子都涂成黑白，没有任何行与列全为黑的方案数。根据容斥，枚举其中全为黑的方格数，可以得到。于是答案便为将的表达式带进去，可以得到设，有记，，这两者可以分别用和处理出来，带入原式得...
复制链接

扫一扫