HDU6129 Just do it[组合数学]

Just do it

  HDU - 6129

 

题意:

给出T组数据，每组含有一个n和m，n代表有n个数，m代表进行m次操作。接下来给出这n个数的值a[i]。

操作是b[i]=a[1]^a[2]^...^a[i]，求出m次操作后b[1]到b[n]的值，每次操作后，用b[i]覆盖a[i]。

题解:

我们以 n=4  a={4,6,5,7}来举例子。(①②...表示下标为 i 的数的个数)

(1) b[1]=4     ① 1

     b[2]=4^6   ① 1   ② 1

     b[3]=4^6^5    ① 1  ② 1  ③ 1

     b[4]=4^6^5^7  ① 1  ② 1  ③ 1  ④1

(2) b[1]=4     ① 1

     b[2]=4^4^6   ① 2   ② 1

     b[3]=4^4^4^6^6^5    ① 3  ② 2  ③ 1

     b[4]=4^4^4^4^6^6^6^5^5^7  ① 4  ② 3  ③ 2  ④1

(3) b[1]=4     ① 1

     b[2]=4^4^4^6   ① 3   ② 1

     b[3]=4^4^4^4^4^4^6^6^6^5    ① 6  ② 3  ③ 1

     b[4]=4^4^4^4^4^4^4^4^4^4^6^6^6^6^6^6^5^5^5^7  ① 10  ② 6  ③ 3  ④1

首先，因为异或的性质，相同的两个数异或会抵消变为0，所以代表着如果这个数在该数的贡献为偶数，那么就不需要进行计算。

其次，我们根据这前三项，可以看出每个数的贡献，是根据杨辉三角的斜边。

从斜边看，第一层到第三层，都是符合我当前的每个数的贡献的。并且我们可以知道 i+m-2就是这斜边的开始第一个行，且从1开始。

普及一个 判断组合数C(n,m)为奇偶的办法 如果n&m==m 代表C(n,m)为奇数，否则为偶数。

那么根据这个，我们可以直接根据组合数的奇偶，来进行异或运算。当当前组合数是奇数时，才去异或。

异或的时候 就等于 b[j]^=a[j-i+1] 因为是对它自身和往前的数有关系。实在弄不懂，仔细观察上面给的三组数据。

#pragma comment(linker, "/STACK:102400000,102400000")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<string>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const int N=2e5+5;
const int M=21;
int a[N],b[N];
int c[M][M];
void initial()
{
    for (int i=0 ; i<M ; ++i)
    {
        c[i][0]=c[i][i]=1;
        for (int j=1 ; j<i ; ++j)
            c[i][j]=c[i-1][j]+c[i-1][j-1];
    }
    for (int i=0 ; i<M ; ++i)
    {
        for (int j=0 ; j<=i ; ++j)
            printf("%6d ",c[i][j]);
        cout<<endl;
    }
}
int main()
{
    //initial();
	int T;
	scanf("%d",&T);
	while (T--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        for (int i=1 ; i<=n ; ++i)
        {
            scanf("%d",&a[i]);
            b[i]=0;
        }
        for (int i=1 ; i<=n ; ++i)
        {
            int up=i+m-2;
            int down=i-1;
            if((up&down)==down)
                for (int j=i ; j<=n ; ++j)
                    b[j]^=a[j-i+1];
        }
        for (int i=1 ; i<=n ; ++i)
        {
            if (i!=1)
                printf(" ");
            printf("%d",b[i]);
        }
        printf("\n");
    }
	return 0;
}