hdu 6129(找规律)

来源:互联网 发布:电脑版地图软件 编辑:程序博客网 时间:2024/05/20 11:24




题目大意:(转载自:http://blog.csdn.net/mengxiang000000/article/details/77200451)

设定b【i】=a【1】^a【2】^a【3】^..................a【i】;

每进行一次,我们可以从a数组得到一个b数组。问进行m次的结果。


思路:


我们随手写下四项的前两次结果,不难看出,我们设定ans【i】【j】表示进行到第i次,第j个位子的答案的话,ans【i】【j】有推导式:

ans【i】【j】=ans【i-1】【j】^ans【i】【j-1】;



那么很显然,对于每一项,他的系数就是杨辉三角的值,那么如果当前位子系数为奇数的话,结果就会有贡献。

同样很显然,我们第i行,第j列的答案,其系数为C(i+j-2,j-1)【此时只考虑a的系数】;


那么我们只需要考虑第一项(a)对所有位子的结果的影响即可(因为b就相当于向后挪了一下递推即可)。

那么根据上述公式,考虑第一项(a)对所有位子的结果的影响然后递推一下就行了。



PS:当初找的时候,只是直接将结果求出来,并没有保留中间的过程,然后gg。要利用到判断C(n,m)的性质),当 n & m == m时,为奇。



#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef pair<int,int> P;#define fi first#define se second#define INF 0x3f3f3f3f#define clr(x,y) memset(x,y,sizeof x)#define PI acos(-1.0)const int maxn = 2e5 + 10;const int Mod = 1e9 + 7;const int N = 2;int a[maxn];int b[maxn];int n,m;void solve(){    clr(b,0);    for(int i = 1; i <= n;i ++)    {        if(( (m + i - 2) & (i - 1)) == (i - 1) )        {            for(int j = i; j <= n; j ++)            {                b[j] ^= a[j - i + 1];            }        }    }    for(int i = 1; i <= n; i ++)        printf("%d%c",b[i],i < n ? ' ':'\n');}int main(){    int Tcase;scanf("%d",&Tcase);    while( Tcase --)    {        scanf("%d%d",&n,&m);for(int i = 1; i <= n;i ++)scanf("%d",&a[i]);        solve();    }    return 0;}



原创粉丝点击