2019-3-22 杂题选讲
T1 CF1137C 图论好题
一个国家有 、 n n n 个城市,通过 m m m 条单向道路相连。有趣的是,在这个国家,每周有 d d d 天,并且每个城市恰好有一个博物馆。
已知每个博物馆一周的营业情况(开门或关门)和 m m m 条单向道路,由于道路的设计,每条道路都需要恰好一个晚上的时间通过。你需要设计一条旅游路线,使得从首都: 1 1 1 号城市开始,并且当天是本周的第一天。每天白天,如果当前城市的博物馆开着门,旅行者可以进入博物馆参观展览,否则什么也做不了,这一天的晚上,旅行者要么结束行程,要么通过一条道路前往下一个城市。当然,旅行者可以多次经过一个城市。
要求让旅行者能够参观的不同博物馆数量尽量多(同一个城市的博物馆参观多次仅算一次),请你求出这个最大值。
正解
按日期拆点后
T
a
r
j
a
n
Tarjan
Tarjan 缩点,在
D
A
G
DAG
DAG 上
D
P
DP
DP 找一条点权最大的路。
唯一要注意的就是在给强连通分量赋新点权时注意一下不重计的细节,一定是该强连通分量里面不同的博物馆的个数。
直接
D
F
S
DFS
DFS 下去搞不定 , 我偷了个小懒,直接加上
i
n
l
i
n
e
inline
inline 就过了。
#include<bits/stdc++.h>
using namespace std;
const int N=100005;
const int D=55;
int n,m,d,in1[N],in2[N];
int point=0,num[N][D],G[N][D];
int tot=0,tot2=0,first[N*D],first2[N*D],nxt[N*D],nxt2[N*D],val[N*D];
int pt=0,cnt=0,top=0,dfn[N*D],vis[N*D],st[N*D],w[N*D],low[N*D],bl[N*D];
int du[N*D],line[N*D],l=1,r=0,dp[N*D],ans=0,tag[N*D];
char s[N][D];
struct Edge
{
int u,v;
}edge[N*D],edge2[N*D];
inline int Read()
{
int x=0,f=1;
char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
return f*x;
}
inline void Add(int u,int v)
{
tot++;
nxt[tot]=first[u];
first[u]=tot;
edge[tot]=(Edge){u,v};
return;
}
inline void Add2(int u,int v)
{
tot2++;
nxt2[tot2]=first2[u];
first2[u]=tot2;
edge2[tot2]=(Edge){u,v};
du[v]++;
return;
}
inline void Init()
{
for (int i=1;i<=n;i++)
for (int j=1;j<=d;j++) num[i][j]=++point;
for (int i=1;i<=n;i++)
for (int j=1;j<=d;j++) G[i][j]=(s[i][j]=='1')?1:0,val[num[i][j]]=G[i][j];
for (int i=1;i<=m;i++)
for (int j=1;j<=d;j++) Add(num[in1[i]][j],num[in2[i]][(j%d)+1]);
return;
}
inline void Tarjan(int u)
{
dfn[u]=low[u]=++cnt;
st[++top]=u;
vis[u]=1;
for (int j=first[u];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (!dfn[v])
{
Tarjan(v);
low[u]=min(low[v],low[u]);
}
else if (vis[v]) low[u]=min(low[u],dfn[v]);
}
if (low[u]==dfn[u])
{
pt++;
while (st[top+1]!=u)
{
vis[st[top]]=0;
bl[st[top]]=pt;
top--;
}
}
return;
}
inline void TP_DP()
{
for (int i=1;i<=pt;i++) if (!du[i]) line[++r]=i,dp[i]=w[i];
while (l<=r)
{
int tmp=line[l++];
for (int j=first2[tmp];j!=-1;j=nxt2[j])
{
int v=edge2[j].v;
dp[v]=max(dp[v],dp[tmp]+w[v]);
du[v]--;
if (!du[v]) line[++r]=v;
}
}
ans=dp[bl[1]];
return;
}
int main()
{
// freopen("testdata.in","r",stdin);
memset(first,-1,sizeof(first));
memset(first2,-1,sizeof(first2));
n=Read(),m=Read(),d=Read();
for (int i=1;i<=m;i++) in1[i]=Read(),in2[i]=Read();
for (int i=1;i<=n;i++) scanf("%s",s[i]+1);
Init();
for (int i=1;i<=point;i++) if (!dfn[i]) Tarjan(i);
for (int i=1;i<=point;i++)
for (int j=first[i];j!=-1;j=nxt[j])
{
int v=edge[j].v;
int from=bl[i],to=bl[v];
if (from!=to) Add2(to,from);
}
for (int i=1;i<=n;i++)
{
for (int j=1;j<=d;j++)
if ((!tag[bl[num[i][j]]])&&(G[i][j])) tag[bl[num[i][j]]]=1,w[bl[num[i][j]]]++;
for (int j=1;j<=d;j++) tag[bl[num[i][j]]]=0;
}
TP_DP();
printf("%d\n",dp[bl[1]]);
return 0;
}
T2 看守 脑筋急转弯
高维曼哈顿 给出n个d维的点,求任意一对点之间曼哈顿距离最大值,高
维的曼哈顿距离定义为每一维坐标的差的绝对值。
n
≤
1
,
000
,
000
n\leq1,000,000
n≤1,000,000
d
≤
4
d\leq4
d≤4
正解
对于每个绝对值符号,枚举所有情况,一共情况数
2
d
2^d
2d
对于每个点把它原本
d
d
d 维的坐标化简成一个数,将一组情况的最大两个数加起来即可。
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
const int D=5;
int n,d,tot;
long long x[N],y[N],ans=0,mx1,mx2,a[N][D];
long long Read()
{
long long x=0,f=1;
char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
return f*x;
}
int main()
{
scanf("%d%d",&n,&d);
for (int i=1;i<=n;i++)
for (int j=1;j<=d;j++) a[i][j]=Read();
tot=(1<<d);
for (int i=0;i<tot;i++)
{
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
for (int j=1;j<=d;j++)
{
int tag=(i>>(j-1))&1;
if (!tag) tag=-1;
for (int k=1;k<=n;k++) x[k]+=a[k][j]*tag,y[k]-=a[k][j]*tag;
}
mx1=-1e18;
mx2=-1e18;
for (int j=1;j<=n;j++) mx1=max(mx1,x[j]),mx2=max(mx2,y[j]);
ans=max(ans,mx1+mx2);
}
printf("%lld\n",ans);
return 0;
}
T3 ???
T4 Power Tree原题 好题
给定一棵树,根是1,每个节点有一个价格,现在你的任务是购买一些节点,使得每个(非根的)叶子节点的数值你可以随意操控(数值不是价格)。对于每个你购买的节点,你可以将这个节点的子树的数值统一变化一个整数
x
x
x。
问:最少要花多少钱可以完成任务;哪些点是某个最优方案可能购买到的。
思考过程
然而给的“差分”提示没有用上,是我想复杂了
玩一下可以发现一个很简单但是对这题没有用 的性质:每个最优方案的买点数等于叶子节点数。
对于第一问的最小花费显然是可以DP的,DP方程也很简单。
设
d
p
[
i
]
[
0
/
1
]
dp[i][0/1]
dp[i][0/1] 表示控制以
i
i
i 为根的子树最小花费是多少。其中 0 表示
i
i
i 的任何一个祖先都没有购买 ,1表示购买了某一个i的祖先。
0和1的状态非常关键,如果是0,代表所有儿子都要“平起平坐”,自己搞好自己的子树,如果是1,则意味某个祖先可以通过祖先的力量调整整颗子树从而使得某一个儿子不需要自己控制自己,可以通过该祖先的力量和其他兄弟的力量 差 使自己的相对数值变化。(可能没说清楚)
所以写出方程:
d
p
[
u
]
[
0
]
=
m
i
n
(
∑
d
p
[
s
o
n
]
[
0
]
,
s
[
u
]
+
d
p
[
s
o
n
]
[
0
]
+
d
p
[
s
o
n
]
[
1
]
+
c
[
u
]
)
dp[u][0]=min ({\sum_ {}dp[son][0], s[u]+dp[son][0]+dp[son][1]+c[u]})
dp[u][0]=min(∑dp[son][0],s[u]+dp[son][0]+dp[son][1]+c[u])
d
p
[
u
]
[
1
]
=
m
i
n
(
s
[
u
]
+
d
p
[
s
o
n
]
[
0
]
+
d
p
[
s
o
n
]
[
1
]
)
dp[u][1]=min ({ s[u]+dp[son][0]+dp[son][1]})
dp[u][1]=min(s[u]+dp[son][0]+dp[son][1])
答案是
d
p
[
1
]
[
0
]
dp[1][0]
dp[1][0]
对于第二问,哪些点可能被选入最优方案。刚才 DP 是从下往上合并的,而决策是从上至下的,于是从上到下再搜一遍,看看哪些状态可能被选中,这里有可能存在多种状态有着同样的最优花费,都要考虑,细节较多。
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+5;
int n,tot=0,first[N],nxt[N<<1],lf[N],d[N],pre[N],vis[N][2],cnt=0,tag[N];
long long c[N],s[N],dp[N][2];
queue <int> q1,q2;
struct Edge
{
int u,v;
}edge[N<<1];
int Read()
{
int x=0,f=1;
char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
return f*x;
}
long long Readl()
{
long long x=0,f=1;
char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
return f*x;
}
void Add(int u,int v)
{
tot++;
nxt[tot]=first[u];
first[u]=tot;
edge[tot]=(Edge){u,v};
d[u]++;
return;
}
void DFS(int u,int fa)
{
if (lf[u])
{
dp[u][0]=c[u];
dp[u][1]=0;
return;
}
for (int j=first[u];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (v==fa) continue;
pre[v]=u;
DFS(v,u);
s[u]+=dp[v][0];
}
dp[u][0]=s[u];
for (int j=first[u];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (v==fa) continue;
dp[u][1]=min(dp[u][1],s[u]-dp[v][0]+dp[v][1]);
dp[u][0]=min(dp[u][0],s[u]-dp[v][0]+dp[v][1]+c[u]);
}
return;
}
void BFS()
{
q1.push(1),q2.push(0);
while (q1.size())
{
int tmp1=q1.front(),tmp2=q2.front();
q1.pop(),q2.pop();
if (lf[tmp1])
{
if (!tmp2) tag[tmp1]=1;
continue;
}
if (dp[tmp1][tmp2]==s[tmp1]&&tmp2==0)
{
for (int j=first[tmp1];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (v==pre[tmp1]) continue;
if (!vis[v][0]) q1.push(v),q2.push(0),vis[v][0]=1;
}
}
int only=-1;
for (int j=first[tmp1];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (v==pre[tmp1]) continue;
if (tmp2==0)
{
if (dp[tmp1][tmp2]==s[tmp1]-dp[v][0]+dp[v][1]+c[tmp1])
{
if (!vis[v][1]) q1.push(v),q2.push(1),vis[v][1]=1;
tag[tmp1]=1;
if (only==1) only=0;
if (only==-1) only=1;
}
if (dp[tmp1][tmp2]!=s[tmp1]-dp[v][0]+dp[v][1]+c[tmp1]) if (!vis[v][0]) q1.push(v),q2.push(0),vis[v][0]=1;
}
if (tmp2==1)
{
if (dp[tmp1][tmp2]==s[tmp1]-dp[v][0]+dp[v][1])
{
if (!vis[v][1]) q1.push(v),q2.push(1),vis[v][1]=1;
if (only==1) only=0;
if (only==-1) only=1;
}
if (dp[tmp1][tmp2]!=s[tmp1]-dp[v][0]+dp[v][1]) if (!vis[v][0]) q1.push(v),q2.push(0),vis[v][0]=1;
}
}
if (only==0)
{
for (int j=first[tmp1];j!=-1;j=nxt[j])
{
int v=edge[j].v;
if (v==pre[tmp1]) continue;
if (!vis[v][0]) q1.push(v),q2.push(0),vis[v][0]=1;
}
}
}
return;
}
int main()
{
// freopen("testdata.in","r",stdin);
memset(first,-1,sizeof(first));
n=Read();
for (int i=1;i<=n;i++) c[i]=Readl(),dp[i][1]=dp[i][0]=1e18;
for (int i=1;i<n;i++)
{
int u,v;
u=Read(),v=Read();
Add(u,v),Add(v,u);
}
for (int i=2;i<=n;i++) if (d[i]==1) lf[i]=1;
DFS(1,0);
BFS();
for (int i=1;i<=n;i++) if (tag[i]) cnt++;
printf("%lld %d\n",dp[1][0],cnt);
for (int i=1;i<=n;i++) if (tag[i]) printf("%d ",i);
return 0;
}