link
C-Half and Half
题意：三种pizza，可以花(A)价钱买一个A-pizza，花(B)价钱买一个B-pizza，花(C*2)价钱买A-pizza和B-pizza各一个，需要(x)个A-pizza和(y)个B-pizza，最小化代价
枚举买多少个AB-pizza即可
#include<bits/stdc++.h>
#define ll long long
#define dbg1(x) cerr<<#x<<"="<<(x)<<" "
#define dbg2(x) cerr<<#x<<"="<<(x)<<"
"
#define dbg3(x) cerr<<#x<<"
"
using namespace std;
#define reg register
inline 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<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
ll ans,A,B,C,X,Y;
ll val(ll i){return i*2ll*C+max(X-i,0ll)*A+max(Y-i,0ll)*B;}
int main()
{
A=read(),B=read(),C=read(),
X=read(),Y=read();
reg int i;
ans=val(0);
for(i=1;i<=X||i<=Y;++i)
ans=min(ans,val(i));
return 0*printf("%lld
",ans);
}
D-Static Sushi
题意：环形吧台，从指定点开始，移动需要耗体力，吃寿司可涨体力，求最大收益
分别枚举顺时针走，逆时针走，先顺后逆，先逆后顺
#include<bits/stdc++.h>
#define ll long long
#define dbg1(x) cerr<<#x<<"="<<(x)<<" "
#define dbg2(x) cerr<<#x<<"="<<(x)<<"
"
#define dbg3(x) cerr<<#x<<"
"
using namespace std;
#define reg register
inline ll read()
{
ll 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<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
const int MN=1e5+5;
ll N,C,X[MN],W[MN];
ll p1,p2,s1,s2,p1m[MN],p2m[MN],s1m[MN],s2m[MN];
int main()
{
N=read(),C=read();
reg int i;
for(i=1;i<=N;++i) X[i]=read(),W[i]=read();
X[N+1]=C;
for(i=1;i<=N;++i)
{
p1+=W[i]+X[i-1]-X[i];
p2=p1-X[i];
p1m[i]=max(p1m[i-1],p1);
p2m[i]=max(p2m[i-1],p2);
}
for(i=N;i;--i)
{
s1+=W[i]-X[i+1]+X[i];
s2=s1-C+X[i];
s1m[i]=max(s1,s1m[i+1]);
s2m[i]=max(s2,s2m[i+1]);
}
ll ans=0;
for(i=1;i<=N;++i)
ans=max(p2m[i]+s1m[i+1],ans),
ans=max(p1m[i],ans);
for(i=1;i<=N;++i)
ans=max(s2m[i]+p1m[i-1],ans),
ans=max(s1m[i],ans);
return 0*printf("%lld
",ans);
}
E-Everything on it
题意：有(n)个元素，一种合法的方案为任意多个集合，每个集合互不相同，(n)个元素的出现次数不少于(2)次，求方案数，对(M)（是质数）取模，(nle 3000)
容斥。
[ans=sum_{i=0}^n(-1)^iinom{n}{i}W(i)
]
(W(i))表示表示至少有(i)个元素出现小于(2)次的方案数
[W(i)=sum_{j=0}^if(i,j)2^{(N-i)j}2^{2^{N-i}}
]
(f(i,j))表示球标号，盒子不标号，球可以不放，盒子非空的方案数，递推可得
#include<bits/stdc++.h>
#define ll long long
#define dbg1(x) cerr<<#x<<"="<<(x)<<" "
#define dbg2(x) cerr<<#x<<"="<<(x)<<"
"
#define dbg3(x) cerr<<#x<<"
"
using namespace std;
#define reg register
inline 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<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
int N,M;
const int MN=3005;
int fac[MN],inv[MN],pw[9000005];
int Mul(int x,int y){return (1ll*x*y)%M;}
int Add(int x,int y){return (x+y)%M;}
int C(int x,int y){return Mul(fac[x],Mul(inv[y],inv[x-y]));}
int f[MN][MN],g[MN],ans;
int main()
{
N=read();M=read();
reg int i,j;
for(pw[0]=i=1;i<=N*N;++i) pw[i]=Mul(2,pw[i-1]);
for(fac[0]=i=1;i<=N;++i)fac[i]=Mul(fac[i-1],i);
for(inv[0]=inv[1]=1,i=2;i<=N;++i)inv[i]=Mul(inv[M%i],(M-M/i));
for(i=2;i<=N;++i)inv[i]=Mul(inv[i-1],inv[i]);
f[0][0]=1;
for(i=1;i<=N;++i)for(f[i][0]=1,j=1;j<=i;++j)
f[i][j]=Add(f[i-1][j],Add(f[i-1][j-1],Mul(f[i-1][j],j)));
reg int tmp=2;
for(i=N;~i;--i,tmp=Mul(tmp,tmp))for(j=0;j<=i;++j)
g[i]=Add(g[i],Mul(f[i][j],Mul(tmp,pw[(N-i)*j])));
for(i=0;i<=N;++i)
ans=Add(ans,Mul((i&1?M-1:1),Mul(C(N,i),g[i])));
return 0*printf("%d
",ans);
}
F-Sweet Alchemy
题意：有棵树，每个点都有代价，要求从上面取点，单个点可取多次，孩子取的次数大于等于父亲的次数且小于等于父亲次数(+D)，代价和不超过(X)，求最多取多少次。
很容易转化题面为：
每个点的代价为子树代价和，价值为子树大小，除(1)点可取inf次外，其它都只能选最多(D)个
除了(n)，以及价值很小（(le 50)），其它都很大
考虑贪心做法：每次取性价比高的，取到不能取为止
但是这么做是有问题的
假设有一种做法，性价比高的还剩下50个，而性价比低已经选了50个，那么可以将他们互换，数量为对方的价值大小，这样贡献不变，体积减小
因为如果性价比高的数有(>50)个没选，那么性价比低的数选的数量一定(<50)，所以将物品分为两个部分，第一个部分每个物品为50个（或者少于50个），(f[i],ile N^3)，表示达到这样价值需要的最少体积，第二个部分为剩下部分，贪心选择性价比高的就可以了
最后枚举第一个部分贡献的价值，求第二个部分的最大价值即可
发现答案一定再枚举的方案中
#include<bits/stdc++.h>
#define ll long long
#define dbg1(x) cerr<<#x<<"="<<(x)<<" "
#define dbg2(x) cerr<<#x<<"="<<(x)<<"
"
#define dbg3(x) cerr<<#x<<"
"
using namespace std;
#define reg register
inline ll read()
{
ll 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<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
const int MN=51;
ll N,X,D,v[MN],p[MN],z[MN],id[MN],nm[MN];
vector<int>s[MN];
void dfs(int x)
{
reg int i;p[x]=1;
for(i=s[x].size()-1;~i;--i)
dfs(s[x][i]),v[x]+=v[s[x][i]],p[x]+=p[s[x][i]];
}
ll f[MN*MN*MN];
bool cmp(const int&x,const int&y){return p[x]*v[y]>p[y]*v[x];}
signed main()
{
N=read(),X=read(),D=read();
reg ll i,j,k;ll lim=0;
for(v[1]=read(),i=2;i<=N;++i)
v[i]=read(),s[read()].push_back(i),z[i]=min(D,N);
dfs(1);z[1]=N;
for(i=1;i<=N;++i)lim+=z[i]*p[i];
for(i=1;i<=lim+1;++i)f[i]=1e18;
for(i=1;i<=N;++i)for(j=lim;j;--j)
{
ll r=min(z[i],j/p[i]);
for(k=1;k<=r;++k)f[j]=min(f[j],f[j-p[i]*k]+v[i]*k);
f[j]=min(f[j],f[j+1]);
}
for(z[1]=1e18,i=2;i<=N;++i)z[i]=D-z[i];
for(i=1;i<=N;++i)id[i]=i;
std::sort(id+1,id+N+1,cmp);
ll ans=0;
for(j=X,i=1;i<=N;++i)
{
nm[i]=min(z[id[i]],j/v[id[i]]);
ans+=p[id[i]]*nm[i];
j-=v[id[i]]*nm[i];
}
ll tmp=ans;
for(j=X-j,k=i,i=1;i<=lim;++i)
{
if(f[i]>X)break;
ll most=X-f[i];
while(j>most)
{
while(!nm[k])--k;
ll n=(j-most+v[id[k]]-1)/v[id[k]];
n=min((ll)n,nm[k]);nm[k]-=n;
j-=v[id[k]]*n;tmp-=p[id[k]]*n;
}
ans=max(ans,tmp+i);
}
printf("%lld
",ans);
return 0;
}
Blog来自PaperCloud，未经允许，请勿转载，TKS！