hdu 6129

这题的规律找起来找到昏迷
假设变化后的数字为 A B C
原先的数字为abc
hdu 6129_第1张图片
括号里的数字为异或的次数
在计算括号的次数的时候就可以发现这个相加规律和杨辉三角类似
然后发现 横的一排其实就是
hdu 6129_第2张图片

观察一下规律
发现第一个数被应用到了每一个组合数 从小到大排过去
第二个数被应用到了(n-1)个数,从本身的位置往右排
也就是说位置很好确定,只要判断组合数为奇==异或的有效改变


补题的时候我傻乎乎的直接算组合数了,实在是太弱智了,m超大的好吗
比赛的时候学姐用lucas过了,去学了一下这个定理
emmmm没搞懂,但知道应用方法了,在网上扒拉了一个lucas模板直接用了(厚颜无耻)
http://www.cnblogs.com/vongang/archive/2012/12/02/2798138.html


放个代码存一下lucas模板

#include
using namespace std;
const int maxn=200005;
long long ans[maxn],a[maxn];
typedef long long LL;


LL exp_mod(LL a, LL b, LL p)
{
    LL res = 1;
    while(b != 0)
    {
        if(b&1) res = (res * a) % p;
        a = (a*a) % p;
        b >>= 1;
    }
    return res;
}

LL Comb(LL a, LL b, LL p)
{
    if(a < b)   return 0;
    if(a == b)  return 1;
    if(b > a - b)   b = a - b;

    LL ans = 1, ca = 1, cb = 1;
    for(LL i = 0; i < b; ++i)
    {
        ca = (ca * (a - i))%p;
        cb = (cb * (b - i))%p;
    }
    ans = (ca*exp_mod(cb, p - 2, p)) % p;
    return ans;
}

LL Lucas(int n, int m, int p)
{
    LL ans = 1;

    while(n&&m&&ans)
    {
        ans = (ans*Comb(n%p, m%p, p)) % p;
        n /= p;
        m /= p;
    }
    return ans;
}

int main()
{
    int cases,i,n,wz,j;
    long long m,now;

    scanf("%d",&cases);

    while (cases--)
    {

        scanf("%d%lld",&n,&m);
        ans[0]=0;
        for(i=1; i<=n; i++)
        {
            scanf("%lld",&a[i]);
            ans[i]=0;
        }


        for (i=0; i//  printf("now=%lld\n",now);
            if (Lucas(m-1+i,m-1,2))
            {
                wz=1;
                for (j=i+1; j<=n; j++)
                {
                    ans[j]=ans[j]^a[wz];
                    wz++;
                }
            }
        }
        printf("%lld",ans[1]);
        for (i=2; i<=n; i++)
            printf(" %lld",ans[i]);
        printf("\n");

    }
    return 0;
}

你可能感兴趣的:(比赛补题,lucas)