消失之物

Posted 2022-05-31 znsbc-13

几种方法

容斥，时间二分，dp+数据结构强艹，主席树

不要拘泥自己思路，思路放开其实不难

主要讲容斥和时间二分

容斥（其实不是容斥，只是类似于容斥）

我们算出来总值，然后减去当前贡献即可

注意一下枚举顺序

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define A 2100
ll f[A],x[A],w[A],ans[A];
const ll mod=10;
ll n,m;
int main()
    scanf("%lld%lld",&n,&m);
    f[0]=1;
    for(ll i=1;i<=n;i++)
        scanf("%lld",&w[i]);
    for(ll i=1;i<=n;i++)
        for(ll j=m;j>=w[i];j--)
            (f[j]+=f[j-w[i]])%=10;
    for(ll i=1;i<=n;i++)
    
        memset(ans,0,sizeof(ans));
        for(ll j=0;j<=m;j++)
            ans[j]=f[j];
        for(ll j=w[i];j<=m;j++)
            ans[j]=(ans[j]-ans[j-w[i]]+10)%10;
        for(ll j=1;j<=m;j++)
            printf("%lld",ans[j]);
        cout<<endl;
    

 

时间二分

skyh代码%%%

比较神的方法

#include<iostream>
#include<cstdio>
#define int short
using namespace std;
const int N=2010;
int n,m,w[N],dp[15][N];
inline void add(int &a,int b)
    a+=b;
    if(a>=10) a-=10;

void solve(int dep,int l,int r)
    if(l==r)
        for(int i=1;i<=m;++i) printf("%hd",dp[dep-1][i]);
        puts("");
        return ;
    
    int mid=l+r>>1;
    for(int i=0;i<=m;++i) dp[dep][i]=dp[dep-1][i];
    for(int i=mid+1;i<=r;++i) for(int j=m;j>=w[i];--j) add(dp[dep][j],dp[dep][j-w[i]]);
    solve(dep+1,l,mid);
    for(int i=0;i<=m;++i) dp[dep][i]=dp[dep-1][i];
    for(int i=l;i<=mid;++i) for(int j=m;j>=w[i];--j) add(dp[dep][j],dp[dep][j-w[i]]);
    solve(dep+1,mid+1,r);

signed main()
    scanf("%hd%hd",&n,&m);
    for(int i=1;i<=n;++i) scanf("%hd",&w[i]);
    dp[0][0]=1; solve(1,1,n);
    return 0;