思路:设d[ i ][ j ][ p ]为在 i 天时,一共工作了 j 天,且连续工作了 p 天直到第 i 天时花费的最小体力。那么转移方程很简单:
update:已经被hacked
如果第 i 天我不工作,那么对于所有的合法的 p,d[ i ][ j ][ 0 ]=max( d[ i-1 ][ j ][ p ] ),如果我第 i 天工作,那么对于所有合法的 j 和 p,d[ i ][ j ][ p ]=d[ i-1 ][ j-1 ][ p-1 ]+p,然后对于所有的d[ i ][ j ][ p ]<=k,取一个最大的 j 当答案即可。
#include<bits/stdc++.h>
using namespace std;
int d[405][405][35],n,k,ans;
char s[160000];
int main()
{
cin>>n>>k>>s+1;
memset(d,-1,sizeof(d));
d[0][0][0]=0;
for(int i=1;i<=n;i++)
{
if(s[i]=='1')
{
for(int j=1;j<=n;j++)
for(int p=1;p*(p+1)/2<=n;p++)
if(d[i-1][j-1][p-1]!=-1)
{
d[i][j][p]=d[i-1][j-1][p-1]+p;
if(d[i][j][p]<=k)
ans=max(ans,j);
}
}
for(int j=0;j<=n;j++)
for(int p=0;p*(p+1)/2<=n;p++)
if(d[i-1][j][p]!=-1)
{
if(d[i][j][0]==-1)d[i][j][0]=d[i-1][j][p];
else
d[i][j][0]=min(d[i-1][j][p],d[i][j][0]);
if(d[i][j][0]<=k)
ans=max(ans,j);
}
}
cout<<ans;
}
update:B题代码被hack了,感谢sugarbliss的指出,hack数据:5 1000 11111,检查发现确实有细节写错了(吐槽牛客数据比较水),还是一样的思路,不过这个地方要改:枚举的p时候是 p*(p+1)/2<=n,改成p*(p+1)/2<=k即可,不过这样的话,dp数组第三维也要开400,会超空间,所以要用小技巧,cur范围为0和1,每次到了一次新的一天,cur取反即可,那么!cur就表示上一天的所有状态,由!cur推cur即可。
#include<bits/stdc++.h>
using namespace std;
int d[2][405][405],n,k,ans,cur;
char s[160000];
int main()
{
cin>>n>>k>>s+1;
memset(d,-1,sizeof(d));
d[0][0][0]=0;
for(int i=1;i<=n;i++)
{
cur^=1;
memset(d[cur],-1,sizeof(d[cur]));
if(s[i]=='1')
{
for(int j=1;j<=n;j++)
for(int p=1;p*(p+1)/2<=k;p++)
if(d[!cur][j-1][p-1]!=-1)
{
d[cur][j][p]=d[!cur][j-1][p-1]+p;
if(d[cur][j][p]<=k)
ans=max(ans,j);
}
}
for(int j=0;j<=n;j++)
for(int p=0;p*(p+1)/2<=k;p++)
if(d[!cur][j][p]!=-1)
{
if(d[cur][j][0]==-1)d[cur][j][0]=d[!cur][j][p];
else
d[cur][j][0]=min(d[!cur][j][p],d[cur][j][0]);
if(d[cur][j][0]<=k)
ans=max(ans,j);
}
}
cout<<ans;
}
思路:比赛时看这题毫无头绪,赛后看别人代码也是看了好久才会,设d[ i ]表示用 i 天完成k个任务得到的最大满意度,那么我可以从k+1枚举到w求解,但是有个很大的问题,初始化d[ k ]的值后,怎么求d[ k+1 ]?我只能通过某一天完成 k-1 个任务才能进行状态转移,但是d[ i ]表示用 i 天完成 k 个任务的最大值,看起来这样设方程是无解的,但是有一个技巧,我令所有的 vi (vi下标从0开始)都减去 v0 ,那么d[ k+1 ]=d[ k ]+v1的含义就可以表示 k-1 天完成 k-1 个任务时(因为去掉了第一个任务),再完成第二个任务所能得到的最大满意度,这样就可以解决状态转移方程了,变成水题了。方程是:d[ i ]=max( d[ i-j ]+v[ j ] ,d[ i ])。
#include<bits/stdc++.h>
using namespace std;
const int maxn=4005,inf=1e8;
int d[maxn],a[maxn],n,k,w;
int main()
{
cin>>n>>k>>w;
for(int i=0;i<n;i++)
cin>>a[i];
for(int i=1;i<=w;i++)d[i]=-inf;
d[k]=k*a[0];
for(int i=1;i<n;i++)a[i]-=a[0];
for(int i=k+1;i<=w;i++)
for(int j=1;j<min(n,i-j+1);j++)
d[i]=max(d[i],d[i-j]+a[j]);
cout<<d[w];
}
update:在北大的一位同学指点下成功补了这题。其实这个题的核心算法是:暴力,每次修改一个区间,根据线段树区间最大值最小值快速找到所有不同的数然后暴力修改他们出现的次数,你或许会有问题:万一区间是 1 到 n,而且区间的每个数都不相同,这样一修改就是遍历了整颗线段树,一次操作复杂度飙升4*n?,但是你要知道,你想要构造这种区间每个数不同的情况,也需要 n次操作才行(这是重点,想一想为什么),那么复杂度均摊下来,还是q*log(n)。(太神仙了,学到了学到了)。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e5+10;
int vis[maxn],tag[maxn*4],ans;
int mx[maxn*4],mn[maxn*4];
void pushdown(int o,int ls,int rs,int l,int r)
{
if(tag[o])
{
mx[o]=mn[o]=tag[o];
if(l!=r)
{
tag[ls]=tag[rs]=tag[o];
mx[ls]=mx[rs]=tag[o];
mn[ls]=mn[rs]=tag[o];
}
tag[o]=0;
}
}
void up(int o,int l,int r,int ql,int qr,int v)
{
int ls=o*2,rs=o*2+1,m=(l+r)/2;
pushdown(o,ls,rs,l,r);
if(l>=ql&&r<=qr&&mx[o]==mn[o])
{
vis[mx[o]]-=(r-l+1);
vis[v]+=(r-l+1);
tag[o]=v;
pushdown(o,ls,rs,l,r);
return;
}
if(ql<=m)up(ls,l,m,ql,qr,v);
if(qr>m)up(rs,m+1,r,ql,qr,v);
mx[o]=max(mx[ls],mx[rs]);
mn[o]=min(mn[ls],mn[rs]);
}
int main()
{
int n,c,q,x,y;
ll a,b;
scanf("%d%d%d",&n,&c,&q);
tag[1]=1;
vis[1]=n;
while(q--)
{
scanf("%d%d%lld%lld",&x,&y,&a,&b);
int p=vis[x];
int l=(a+(p+b)%n*(p+b)%n)%n+1;
int r=(a+p*b%n*p%n*b%n)%n+1;
if(l>r)swap(l,r);
up(1,1,n,l,r,y);
}
for(int i=1;i<=c;i++)
ans=max(ans,vis[i]);
printf("%d\n",ans);
}
求任意凸多边形的面积公式:(逆时针顺序遍历点)
[(x1*y2+x2*y3+x3*y4+x4*y5+x5*y1)-(x1*y5+x2*y1+x3*y2+x4*y3+x5*y4)]/2
这题先求出xi*yi+1,以及yi*xi+1的前缀和数组sum1,sum2,然后怎么快速求 s t 切割的多边形面积呢?其实很简单,随便画一下就知道知道了,我就不写了。
#include<bits/stdc++.h>
#define db double
using namespace std;
const int maxn=1e5+10;
db x[maxn],y[maxn],sum1[maxn],sum2[maxn],ans;
int main()
{
int n,q,s,t;
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++)
scanf("%lf%lf",&x[i],&y[i]);
y[n+1]=y[1],x[n+1]=x[1];
for(int i=1;i<=n;i++)
{
sum1[i]=sum1[i-1]+x[i]*y[i+1];
sum2[i]=sum2[i-1]+y[i]*x[i+1];
}
db S=(sum1[n]-sum2[n])/2;
while(q--)
{
scanf("%d%d",&s,&t);
if(s>t)swap(s,t);
db t1=sum1[t-1]-sum1[s-1]+x[t]*y[s];
db t2=sum2[t-1]-sum2[s-1]+y[t]*x[s];
db tmp=(t1-t2)/2;
ans=max(ans,min(tmp,S-tmp));
}
printf("%.2lf\n",ans);
}