reimu hakurei

题目名意义

博丽灵梦
本弱表示从来自机都选的是博丽灵梦QAQ
（当然像红魔乡附加关那种就只能选雾雨魔理沙了

题目描述

构造

首先很显然
$f[i]=\sum_{j|i}b[j]*d(\frac{i}{j})$
d是除数个数函数。
先把d处理出来。
一个显然的结论：一定能构造出价值为满的
b都是正整数，改变i的选或不选状态只对后面产生影响。
从前往后做，如果f[i]和r[i]相同，就选上i。
于是构造成功。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int maxn=200000+10;
int b[maxn],d[maxn];
ll f[maxn],r[maxn];
bool bz[maxn];
int i,j,k,l,t,n,m,sum;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
int main(){
    n=read();
    fo(i,1,n)
        fo(j,1,n/i)
            d[i*j]++;
    fo(i,1,n) b[i]=read();
    fo(i,1,n) r[i]=read();
    fo(i,1,n){
        t=read();
        sum+=t;
    }
    printf("%d\n",sum);
    fo(i,1,n)
        if (f[i]==r[i]){
            bz[i]=1;
            fo(j,1,n/i)
                f[i*j]+=(ll)b[i]*d[j];
        }
    fo(i,1,n) printf("%d%c",bz[i],i==n?'\n':' ');
}