六省联考2017 Day2
2018.3.27 Test
时间：7:30~11:50
期望得分：(50+)+0+20=70
实际得分：52+5+20=77
总结
T1 看错一点题，暴力也废了很长时间。
T2 期望DP没写过不敢写，然而50分和期望没有关系，贪心什么的就行。没细看。
T3 建图死活建不出来，没想明白费用流还费了不少时间写费用流。
T1
好像很麻烦，先弃疗。
T2
代码倒很短，学学期望再写。
T3 BZOJ.4873.[六省联考2017]寿司餐厅(最小割ISAP 最大权闭合子图)
不会建图。。在最大权闭合子图和费用流徘徊，快考完发现费用流不对。。硬是把区间单独考虑了，没考虑与小区间的限制关系。
正解: 最大权闭合子图.
如果想要取一个区间，那么会把这一区间中所有的都取一遍，且一共只取一遍。
即: 选(i,j)(i<j)必选(i+1,j)和(i,j-1)；选(i,i)必选id[i]。
由源点向每个权值为正的区间连权值的边；由每个权值为负的区间向汇点连权值绝对值的边(阻碍割掉连向汇点的一边)；
用每种寿司i直接代表区间(i,i)，根据D[i,i]-id[i]的正负向汇点/源点连边；寿司向id[i]连INF的边，id[i]向汇点连m*id^2的边。
Ans = 正权和 - 最大流
//2240kb 148ms(怎么好像很慢。。)
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
const int N=1e4+5,M=7e4+5,INF=0x3f3f3f3f;
int n,m,src,des,Enum,cur[N],H[N],nxt[M],fr[M],to[M],cap[M],lev[N],num[N],que[N],pre[N];
int id[105][105],A[105],d[105][105],exist[1005];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline void AddEdge(int u,int v,int w)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
}
bool BFS()
{
for(int i=src; i<des; ++i) lev[i]=des+1;
que[0]=des, lev[des]=0; int h=0,t=1;
while(h<t)
{
int x=que[h++];
for(int i=H[x]; i; i=nxt[i])
if(lev[to[i]]==des+1 && cap[i^1])
lev[to[i]]=lev[x]+1, que[t++]=to[i];
}
return lev[src]<=des;
}
int Augment()
{
int mn=INF;
for(int i=des; i!=src; i=fr[pre[i]])
mn=std::min(mn,cap[pre[i]]);
for(int i=des; i!=src; i=fr[pre[i]])
cap[pre[i]]-=mn,cap[pre[i]^1]+=mn;
return mn;
}
long long ISAP()
{
if(!BFS()) return 0;
for(int i=src; i<=des; ++i) cur[i]=H[i],++num[lev[i]];
int x=src; long long res=0;
while(lev[src]<=des)
{
if(x==des) x=src,res+=Augment();
bool can=0;
for(int i=cur[x]; i; i=nxt[i])
if(lev[to[i]]==lev[x]-1 && cap[i])
{
can=1, cur[x]=i, pre[x=to[i]]=i;
break;
}
if(!can)
{
int mn=des;
for(int i=H[x]; i; i=nxt[i])
if(cap[i]) mn=std::min(mn,lev[to[i]]);
if(!--num[lev[x]]) break;
++num[lev[x]=mn+1];
cur[x]=H[x];
if(x!=src) x=fr[pre[x]];
}
}
return res;
}
int main()
{
n=read(),m=read(),Enum=1;
int cnt=n;
for(int i=1; i<=n; ++i) id[i][i]=i;
for(int i=1; i<n; ++i)
for(int j=i+1; j<=n; ++j) id[i][j]=++cnt;
src=0, des=cnt+n+1;
for(int i=1; i<=n; ++i)
{
A[i]=read();
if(m && !exist[A[i]]) exist[A[i]]=++cnt,AddEdge(cnt,des,/*m**/A[i]*A[i]);
if(m) AddEdge(i,exist[A[i]],INF);
}
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j) d[i][j]=read();
long long res=0;
for(int t,i=1; i<=n; ++i)
(t=d[i][i]-A[i])>0 ? res+=t,AddEdge(src,i,t) : AddEdge(i,des,-t);
for(int i=1; i<n; ++i)
for(int j=i+1; j<=n; ++j)
{
AddEdge(id[i][j],id[i+1][j],INF), AddEdge(id[i][j],id[i][j-1],INF);
d[i][j]>0 ? res+=d[i][j],AddEdge(src,id[i][j],d[i][j]) : AddEdge(id[i][j],des,-d[i][j]);
}
printf("%lld",res-ISAP());
return 0;
}
考试代码
T1
5k暴力心累（一半复制的）
刚开始想错题意了，判的点不能走。当然挺好改。
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int MAXIN=1e6;
int T,X;
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
namespace Subtask1
{
const int N=9;
int n,cnt,pos1[1<<9],pos2[1<<9],H[N],Enum,fr[N<<1],to[N<<1],nxt[N<<1],pre[N];
bool can,vis[N<<1],vis2[N<<1],visp[N],visp2[N];
inline void AddEdge(int u,int v)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum;
}
void Pre()
{
for(int t,p,i=1; i<(1<<8); ++i)
{
t=0, p=cnt+1;
for(int j=i,l=1; j; j>>=1,++l)
if(j&1){
if(++t>2) break;
if(pos1[p]) pos2[p]=l;
else pos1[p]=l;
}
if(t==2) ++cnt;
else if(t==1) ++cnt,pos2[p]=pos1[p];
else pos1[p]=pos2[p]=0;
}
// for(int i=1; i<=cnt; ++i) printf("Pre:%d -> %d\n",pos1[i],pos2[i]);
}
void DFS(int x,int f,int t)
{
if(x==t) {can=1; return;}
for(int i=H[x]; i&&!can; i=nxt[i])
if(to[i]!=f&&!vis[i]) pre[to[i]]=i,DFS(to[i],x,t);
}
void Calc_DFS(int x)
{
visp2[x]=1;
for(int i=H[x]; i; i=nxt[i])
if(!vis[i]&&!vis2[i]&&!visp[to[i]]&&!visp2[to[i]]) Calc_DFS(to[i]);
}
void Solve()
{
Pre();
while(T--)
{
n=read(), Enum=0;
for(int i=1; i<=n; ++i) H[i]=0;
for(int u,v,i=1; i<n; ++i) u=read(),v=read(),AddEdge(u,v);
int res=0;
for(int i=1; i<=cnt; ++i)
{
if(pos1[i]>n||pos2[i]>n) continue;
for(int j=1; j<=Enum; ++j) vis[j]=0;
for(int j=1; j<=n; ++j) visp[j]=0;
can=0, DFS(pos1[i],0,pos2[i]);
int t=pos2[i];
while(t!=pos1[i]) visp[t]=vis[pre[t]]=1,t=fr[pre[t]];
visp[pos1[i]]=1;
for(int ans,j=1; j<=cnt; ++j)
{
if(pos1[j]>n||pos2[j]>n) continue;
can=0, DFS(pos1[j],0,pos2[j]);
if(!can) continue;
for(int k=1; k<=Enum; ++k) vis2[k]=0;
for(int k=1; k<=n; ++k) visp2[k]=0;
t=pos2[j];
while(t!=pos1[j]) visp2[t]=vis2[pre[t]]=1,t=fr[pre[t]];
visp2[pos1[j]]=1;
ans=0;
for(int k=1; k<=n; ++k)
if(!visp[k]&&!visp2[k]) Calc_DFS(k),++ans;
if(ans>res) res=ans;
}
}
printf("%d\n",res);
}
}
}
namespace Subtask2
{
const int N=1003;
int n,p0,p1,h0,h1,H[N],Enum,fr[N<<1],to[N<<1],nxt[N<<1],pre[N];
bool can,vis[N<<1],vis2[N<<1],visp[N],visp2[N];
inline void AddEdge(int u,int v)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum;
}
void DFS(int x,int f,int t)
{
if(x==t) {can=1; return;}
for(int i=H[x]; i&&!can; i=nxt[i])
if(to[i]!=f&&!vis[i]) pre[to[i]]=i,DFS(to[i],x,t);
}
void Calc_DFS(int x)
{
visp2[x]=1;
for(int i=H[x]; i; i=nxt[i])
if(!vis[i]&&!vis2[i]&&!visp[to[i]]&&!visp2[to[i]]) Calc_DFS(to[i]);
}
void Solve()
{
while(T--)
{
n=read(),p0=read(),p1=read(),h0=read(),h1=read(), Enum=0;
for(int i=1; i<=n; ++i) H[i]=0;
for(int u,v,i=1; i<n; ++i) u=read(),v=read(),AddEdge(u,v);
for(int i=1; i<=Enum; ++i) vis[i]=0;
for(int i=1; i<=n; ++i) visp[i]=0;
can=0, DFS(p0,0,p1);
int t=p1;
while(t!=p0) visp[t]=vis[pre[t]]=1,t=fr[pre[t]];
visp[p0]=1;
can=0, DFS(h0,0,h1);
if(!can) {puts("0"); continue;}
for(int k=1; k<=Enum; ++k) vis2[k]=0;
for(int k=1; k<=n; ++k) visp2[k]=0;
t=h1;
while(t!=h0) visp2[t]=vis2[pre[t]]=1,t=fr[pre[t]];
visp2[h0]=1;
int res=0;
for(int k=1; k<=n; ++k)
if(!visp[k]&&!visp2[k]) Calc_DFS(k),++res;
printf("%d\n",res);
}
}
}
namespace Subtask3
{
const int N=1003;
int n,p0,p1,H[N],Enum,fr[N<<1],to[N<<1],nxt[N<<1],pre[N];
bool can,vis[N<<1],vis2[N<<1],visp[N],visp2[N];
inline void AddEdge(int u,int v)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum;
}
void DFS(int x,int f,int t)
{
if(x==t) {can=1; return;}
for(int i=H[x]; i&&!can; i=nxt[i])
if(to[i]!=f&&!vis[i]) pre[to[i]]=i,DFS(to[i],x,t);
}
void Calc_DFS(int x)
{
visp2[x]=1;
for(int i=H[x]; i; i=nxt[i])
if(!vis[i]&&!vis2[i]&&!visp[to[i]]&&!visp2[to[i]]) Calc_DFS(to[i]);
}
void Solve()
{
while(T--)
{
n=read(),p0=read(),p1=read(), Enum=0;
for(int i=1; i<=n; ++i) H[i]=0;
for(int u,v,i=1; i<n; ++i) u=read(),v=read(),AddEdge(u,v);
for(int i=1; i<=Enum; ++i) vis[i]=0;
for(int i=1; i<=n; ++i) visp[i]=0;
can=0, DFS(p0,0,p1);
int t=p1;
while(t!=p0) visp[t]=vis[pre[t]]=1,t=fr[pre[t]];
visp[p0]=1;
int res=0;
for(int i=1; i<=n; ++i)
for(int j=1; j<=n; ++j)
{
can=0, DFS(i,0,j);
if(!can) continue;
for(int k=1; k<=Enum; ++k) vis2[k]=0;
for(int k=1; k<=n; ++k) visp2[k]=0;
t=j;
while(t!=i) visp2[t]=vis2[pre[t]]=1,t=fr[pre[t]];
visp2[i]=1;
int ans=0;
for(int k=1; k<=n; ++k)
if(!visp[k]&&!visp2[k]) Calc_DFS(k),++ans;
if(ans>res) res=ans;
}
printf("%d\n",res);
}
}
}
int main()
{
freopen("treediagram.in","r",stdin);
freopen("treediagram.out","w",stdout);
T=read(),X=read();
if(!X) Subtask1::Solve();
else if(X==2) Subtask2::Solve();
else Subtask3::Solve();
fclose(stdin);fclose(stdout);
return 0;
}
T2
直接输出0真的可以，还有5分233→_→
#include <cstdio>
#include <cctype>
#define gc() getchar()
#define mod (100003)
const int N=1e3+5;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
int main()
{
freopen("trennen.in","r",stdin);
freopen("trennen.out","w",stdout);
putchar('0');
fclose(stdin);fclose(stdout);
return 0;
}
T3
爆搜：
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=2e4+5,M=3e6+5,INF=0x3f3f3f3f;
int n,m,Ans,Max;
int cnt,tm[35],id[105][105],A[105],l[55555],r[55555],cost[45555],sum[105],d[105][105];
bool vis[7][7];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
void DFS(int x,int s)
{
if(x>cnt)
{
memset(tm,0,sizeof tm);
memset(vis,0,sizeof vis);
int tmp=0;
for(int i=1; i<=cnt; ++i)
if((s>>i-1)&1)
for(int j=l[i]; j<=r[i]; ++j)
{
for(int k=j; k<=r[i]; ++k)
if(!vis[j][k])
tmp+=d[j][k], vis[j][k]=1;
++tm[A[j]];
}
for(int i=1; i<=Max; ++i) if(tm[i]) tmp-=m*i+tm[i]*i;
Ans=std::max(Ans,tmp);
}
else DFS(x+1,s), DFS(x+1,s|(1<<x-1));
}
int main()
{
freopen("sushi.in","r",stdin);
freopen("sushi.out","w",stdout);
n=read(),m=read();
for(int i=1; i<=n; ++i) Max=std::max(Max,A[i]=read()),sum[i]=sum[i-1]+A[i];
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j) d[i][j]=read();
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j) id[i][j]=++cnt,cost[cnt]=sum[j]-sum[i-1],l[cnt]=i,r[cnt]=j;
DFS(1,0);
printf("%d",Ans);
fclose(stdin);fclose(stdout);
return 0;
}
sb费用流。。（还是由最大流又改成费用流的）
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=2e4+5,M=3e6+5,INF=0x3f3f3f3f;
int n,m,src,des,Enum,H[N],nxt[M],fr[M],to[M],cap[M],cost[M],dis[N],pre[N];
int cnt,id[105][105],A[105],sum[105],d[105][105];
bool inq[N];
std::queue<int> q;
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline void AddEdge(int u,int v,int w,int c)
{
printf("%d->%d cap:%d c:%d\n",u,v,w,c);
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w, cost[Enum]=c;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0, cost[Enum]=-c;
}
bool SPFA()
{
memset(dis,0x3f,sizeof dis);
dis[src]=0, q.push(src);
while(!q.empty())
{
int x=q.front(); q.pop(),inq[x]=0;
for(int i=H[x]; i; i=nxt[i])
if(dis[to[i]]>dis[x]+cost[i] && cap[i])
{
dis[to[i]]=dis[x]+cost[i], pre[to[i]]=i;
if(!inq[to[i]]) q.push(to[i]),inq[to[i]]=1;
}
}
return dis[des]!=INF;
}
int MCMF()
{
int mn=INF,res=0;
for(int i=des; i!=src; i=fr[pre[i]])
mn=std::min(mn,cap[pre[i]]);
for(int i=des; i!=src; i=fr[pre[i]])
cap[pre[i]]-=mn,cap[pre[i]^1]+=mn,res+=cost[pre[i]]*mn;
for(int i=des; i!=src; i=fr[pre[i]])
printf("%d<-",i);
printf("%d flow:%d res:%d\n",src,mn,res);
return res;
}
int main()
{
freopen("sushi.in","r",stdin);
// freopen("sushi.out","w",stdout);
n=read(),m=read(),Enum=1;
for(int i=1; i<=n; ++i) A[i]=read(),sum[i]=sum[i-1]+A[i];
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j)
d[i][j]=read();
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j) id[i][j]=++cnt;
src=0, des=cnt<<1|1;
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j)
{
int t=0;
for(int k=i; k<=j; ++k)
for(int l=k; l<=j; ++l)
++t, AddEdge(id[i][j],id[k][l]+cnt,INF,0);
AddEdge(src,id[i][j],t,sum[j]-sum[i-1]);
}
for(int i=1; i<=n; ++i)
for(int j=i; j<=n; ++j) AddEdge(id[i][j]+cnt,des,1,-d[i][j]);
long long res=0,ans=1e9;
while(SPFA()) ans=std::min(ans,res+=MCMF()),printf("%I64d\n",res);
printf("%I64d",-ans);
fclose(stdin);fclose(stdout);
return 0;
}
六省联考2017 Day2的更多相关文章
- 【BZOJ4873】[六省联考2017]寿司餐厅（网络流）
[BZOJ4873][六省联考2017]寿司餐厅(网络流) 题面 BZOJ 洛谷 题解 很有意思的题目 首先看到答案的计算方法,就很明显的感觉到是一个最大权闭合子图. 然后只需要考虑怎么构图就行了. ...
- 【BZOJ4868】[六省联考2017]期末考试（贪心）
[BZOJ4868][六省联考2017]期末考试(贪心) 题面 BZOJ 洛谷 题解 显然最终的答案之和最后一个公布成绩的课程相关. 枚举最后一天的日期,那么维护一下前面有多少天可以向后移,后面总共需 ...
- 六省联考2017 Day1
目录 2018.3.18 Test T1 BZOJ.4868.[六省联考2017]期末考试 T2 T3 BZOJ.4870.[六省联考2017]组合数问题(DP 矩阵快速幂) 总结 考试代码 T1 T ...
- 洛谷 P3747 [六省联考2017]相逢是问候 解题报告
P3747 [六省联考2017]相逢是问候 题目描述 \(\text {Informatik verbindet dich und mich.}\) 信息将你我连结. \(B\) 君希望以维护一个长度 ...
- P3746 [六省联考2017]组合数问题
P3746 [六省联考2017]组合数问题 \(dp_{i,j}\)表示前\(i\)个物品,取的物品模\(k\)等于\(r\),则\(dp_{i,j}=dp_{i-1,(j-1+k)\%k}+dp_{ ...
- [luogu] P3745 [六省联考2017]期末考试 （贪心）
P3745 [六省联考2017]期末考试 题目描述 有 \(n\) 位同学,每位同学都参加了全部的 \(m\) 门课程的期末考试,都在焦急的等待成绩的公布. 第 \(i\) 位同学希望在第 \(t_i ...
- BZOJ 4872 luogu P3750 [六省联考2017]分手是祝愿
4872: [Shoi2017]分手是祝愿 Time Limit: 20 Sec Memory Limit: 512 MB[Submit][Status][Discuss] Description ...
- bzoj千题计划266：bzoj4872: [六省联考2017]分手是祝愿
http://www.lydsy.com/JudgeOnline/problem.php?id=4872 一种最优解是 从大到小灯有亮的就灭掉 最优解是唯一的,且关灯的顺序没有影响 最优解 对每个开关 ...
- bzoj千题计划265：bzoj4873: [六省联考2017]寿司餐厅
http://www.lydsy.com/JudgeOnline/problem.php?id=4873 选a必选b,a依赖于b 最大权闭合子图模型 构图: 1.源点 向 正美味度区间 连 流量为 美 ...
随机推荐
- python命名空间与闭包函数详解
python命名空间与闭包函数详解 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 本篇博客主要介绍的知识点如下: 1>.三元运算 2>.命名空间 3>.globa ...
- bzoj千题计划281：bzoj4558: [JLoi2016]方
http://www.lydsy.com/JudgeOnline/problem.php?id=4558 容斥原理 全部的正方形-至少有一个点被删掉的+至少有两个点被删掉的-至少有3个点被删掉的+至少 ...
- hdu 5956 The Elder
http://acm.hdu.edu.cn/showproblem.php?pid=5956 转移方程:dp[i]=(dis[i]-dis[j])*(dis[i]-dis[j])+P+dp[j] 斜率 ...
- plsql免安装客户端的配置
不安装oracle,在安装了plsql之后,需要连接数据库,连接数据库需要在tns中tnsnames.ora中配置 首先需要两个文件: network instantclient-basic-win3 ...
- 洛谷 P1609 最小回文数 题解
这题其实并不难,重点在你对回文数的了解,根本就不需要高精度. 打个比方: 对于一个形如 ABCDEFGH 的整数 有且仅有一个比它大的最小回文数 有且仅有一个比它小的最大回文数 而整数 ABCDDCB ...
- jinja模板语言
模板 要了解jinja2,那么需要先理解模板的概念.模板在Python的web开发中广泛使用,它能够有效的将业务逻辑和页面逻辑分开,使代码可读性增强.并且更加容易理解和维护. 模板简单来说就是一个其中 ...
- insert into与insert ignore以及replace into的区别
insert ignore表示,如果表中已经存在相同的记录,则忽略当前新数据: INSERT INTO有无数据都插入,如果主键则不插入; REPLACE INTO 如果是主键插入则会替换以前的数据; ...
- Html设置html与body元素高度问题
为什么要设置HTML和body的高度? 在一些比较正规的网站经常见到会设置html与body的高度.是为了后面的div高度正确的显示. div的100%是从其上一级div的宽高继承来的,有一点很关键, ...
- flask基础之Response响应对象（九）
前言 Response对象负责对客户端的响应,每一个请求都会有一个Response对象,那么它在一个请求的声明周期内是怎么发挥作用的呢? Response对象 响应发生的位置 先回顾一下http请求的 ...
- find中的-print0和xargs中-0的奥妙【转】
find cygnus/firmware_cygnus/target/linux/brcm5830/files/arch/arm/mach-iproc/pm_iproc/ -name "*. ...