Codeforces 241D : Numbers

【大意】
给定一个1到n的排列a和一个素数p,构造一个子序列,使得元素异或和为0且首尾相接形成的大整数是p的倍数。n,p 50000。
【解答】
做法是把原数列中小于32的数拿出来构造这个序列,无视掉序列中大于等于32的数

小于等于 32的数中,异或和为0的组合有 (312)4
当p不为2或者5时,首尾相接拼起来模p为0可以看做是一个随机函数,答案出现的概率是很高的
当p为2或者5时,手玩发现是存在解的。

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

int n,p,tp,tot,w[60],ys[N],a[N],b[N],ans[N];
bool F[35][35][N][2];
inline int rd() {int r;scanf("%d",&r);return r;}

inline int qp(int a,int b) {
    int ret = 1;
    while (b) {
        if (b&1) ret = 1LL * ret * a % p;
        b >>= 1, a = 1LL * a * a % p;
    }
    return ret;
}

void o(int x,int y,int z,int c) {
    if (!x) return ;
    if (F[x-1][y][z][c]) {
        o(x-1,y,z,c);
        return ;
    }

    ans[++tot] = ys[x];
    for (int i=0;i<p;i++)
        if ((i*w[ b[x] ]+b[x])%p == z) {
            if (F[x-1][y^b[x]][i][0]) {
                o(x-1,y^b[x],i,0);
                return ;
            }
            if (F[x-1][y^b[x]][i][1]) {
                o(x-1,y^b[x],i,1);
                return ;
            }
        }
}

int calc(int x) {
    int ret = 10;
    while (x>=10) x /= 10, ret*=10;
    return ret;
}

int main() {
    n = rd(), p = rd();
    for (int i=1;i<=n;i++) a[i] = rd();

    for (int i=1;i<=n;i++) if (a[i]<32) b[++tp] = a[i], ys[tp] = i;
    for (int i=0;i<=32;i++) w[i] = calc(i);

    F[0][0][0][0] = 1;
    for (int _=0;_<tp;_++) {
        int x = b[_+1];
        for (int i=0;i<=32;i++) {
            for (int j=0;j<p;j++) {
                if (F[_][i][j][1]) {
                    F[_+1][i^x][(j*w[x]+x)%p][1] = 1;
                    F[_+1][i][j][1] = 1;
                }
                if (F[_][i][j][0]) {
                    F[_+1][i^x][(j*w[x]+x)%p][1] = 1;
                    F[_+1][i][j][0] = 1;
                }
            }
        }
    }

    for (int _=1;_<=tp;_++) if (F[_][0][0][1]) {
        puts("Yes");
        o(_,0,0,1);
        printf("%d\n",tot);
        for (int _=tot;_>=1;_--) printf("%d%c",ans[_],_==1?'\n':' ');
        return 0;
    }
    puts("No");
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值