HDU 6513 Reverse It(组合数学+消去重复度)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6513

题目大意

给定一个01矩阵,现在可以选择至多两个子矩阵

把值进行转换(0-1,1-0),

问能产生总共多少种不同的花纹数.

题目分析 

首先考虑两个子矩阵和一个子矩阵,

不难发现一个子矩阵就是:C(n+1,2)*C(m+1,2),

因为n个长度其格点有n+1个.

而两个子矩阵就是在前面的一个子矩阵中选择两个,

主要是重复度的消去,

就是说对于同一种花纹即可以通过一个子矩阵计算,也可以通过两个子矩阵计算出来.

主要是几种情况:

1,两个子矩阵在同行或同列,这样的情况组合数是:C(n+1,2)*C(m+1,4)+C(m+1,2)*C(n+1,4),

重复度是:三个,不难发现我们可以通过三种方式组合这样的花纹,所以消去2倍即可.

2,L型,这样的花纹可以由三种情况组成,而对于每种组合有四种旋转,

所以组合数是C(n+1,3)*C(m+1,3),重复度是:8.

3,对角型,这样的情况分析和上面是一样的,组合数是C(n+1,3)*C(m+1,3),

重复度是:4.

4,单个矩阵型,对于单个小矩阵,即C(n+1,2)*C(m+1,2),

不难枚举发现我们可以有n+m-2种方式得到这样的情况,为方便计算,

我们直接 令答案总数为两个子矩阵情况+1,

所以这样的形状重复度是(n+m-3).

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=1e2+100;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){
    if(y==0) return x;
    return gcd(y,x%y);
}
/*
题目大意:
给定一个01矩阵,现在可以选择至多两个子矩阵
把值进行转换(0-1,1-0),
问能产生总共多少种不同的花纹数.

题目分析:
首先考虑两个子矩阵和一个子矩阵,
不难发现一个子矩阵就是:C(n+1,2)*C(m+1,2),
因为n个长度其格点有n+1个.
而两个子矩阵就是在前面的一个子矩阵中选择两个,
主要是重复度的消去,
就是说对于同一种花纹即可以通过一个子矩阵计算,也可以通过两个子矩阵计算出来.
主要是几种情况:
1,两个子矩阵在同行或同列,这样的情况组合数是:C(n+1,2)*C(m+1,4)+C(m+1,2)*C(n+1,4),
重复度是:三个,不难发现我们可以通过三种方式组合这样的花纹,所以消去2倍即可.
2,L型,这样的花纹可以由三种情况组成,而对于每种组合有四种旋转,
所以组合数是C(n+1,3)*C(m+1,3),重复度是:8.
3,对角型,这样的情况分析和上面是一样的,组合数是C(n+1,3)*C(m+1,3),
重复度是:4.
4,单个矩阵型,对于单个小矩阵,即C(n+1,2)*C(m+1,2),
不难枚举发现我们可以有n+m-2种方式得到这样的情况,为方便计算,
我们直接 令答案总数为两个子矩阵情况+1,
所以这样的形状重复度是(n+m-3).
*/
ll n,m;
char s[maxn][maxn];
ll C(ll x,int y){
    if(x<1LL*y) return 0;
    ll ans=1,tmp=1;
    for(int i=1;i<=y;i++){
        ans=ans*(x-i+1);
        tmp*=i;
    }
    return ans/tmp;
}
int main(){
    while(scanf("%lld%lld",&n,&m)!=EOF){
        rep(i,0,n) scanf("%s",s[i]);
        ll ret=C(n+1,2)*C(m+1,2);
        ret=C(ret,2)+1;
        ret-=4LL*C(n+1,3)*C(m+1,3);///对角
        ret-=8LL*C(n+1,3)*C(m+1,3);///斜率
        ret-=(2LL*C(n+1,2)*C(m+1,4)+2LL*C(n+1,4)*C(m+1,2));///分割
        ret-=(n+m-3)*(C(n+1,2)*C(m+1,2));
        printf("%lld\n",ret);
    }
    return 0;
}