思路:裸的状压dp,大家可以搜索旅行商问题状压dp解法,和这个一样。
设d[ u ][ s ] 为 u 排在第一位,排队的人的集合为 s 所能得到的最优值,那么对于当前的d[ j ][ i ],我们找一个集合 i 中没有出现的同学 k,那么令 k 排在第一位得到转移方程:d[ k ][ s | (1<<k)]=max( d[ k ][ s | (1<<k)], d[ j ][ s ]+w[ k ][ k ]-w[ j ][ j ]+w[ j ][ k ] )。
#include<bits/stdc++.h>
#define ll long long
#define db double
using namespace std;
const int maxn=1<<18;
int d[18][maxn],a[20][20];
void up(int &x,int y)
{
x=max(x,y);
}
int main()
{
int n,ans=0;
cin>>n;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
cin>>a[i][j];
if(i==j)
d[i][1<<i]=a[i][j];
}
for(int i=1;i<(1<<n);i++)
for(int j=0;j<n;j++)
{
if(i==(1<<n)-1)
ans=max(ans,d[j][i]);
if(((1<<j)&i))
{
for(int k=0;k<n;k++)
if(!((1<<k)&i))
up(d[k][i|(1<<k)],d[j][i]+a[k][k]-a[j][j]+a[j][k]);
}
}
printf("%d\n",ans);
}
思路:由于树上每条边的距离都相等,我们反过来求不能组成三角形的概率,很容易想到,当三个点在树上同一条链上时,这三个点不能组成三角形,那么我树形dp,枚举每个点作为中间点所能得到非法三角形的概率即可。
#include<bits/stdc++.h>
#define ll long long
#define db double
using namespace std;
const int maxn=1e5+10;
int d[maxn],n;
ll m;
db ans;
vector<int>G[maxn];
void dfs(int u,int fa)
{
for(int i=0;i<G[u].size();i++)
{
int v=G[u][i];
if(v==fa)continue;
dfs(v,u);
ans+=1.0*d[u]*d[v]/m;
d[u]+=d[v];
}
d[u]++;
ans+=1.0*(d[u]-1)*(n-d[u])/m;
}
int main()
{
int u,v;
scanf("%d",&n);
m=1ll*n*(n-1)*(n-2)/6;
for(int i=1;i<n;i++)
{
scanf("%d%d",&u,&v);
G[u].push_back(v);
G[v].push_back(u);
}
dfs(1,0);
printf("%.10lf\n",1.0-ans);
}
G.Chino with Train to the Rabbit Town
思路:这个题和某道cf题相似,说白了,求最多有多少个不相交区间,区间异或=t,我们枚举区间右端点 r,假设我已经求了并保存了[ 1 , r],[ 2 ,r ],[ 3 , r ],,,,,[ r ,r ]的区间异或值,我现在枚举r+1,要重新更新前面 r 个区间异或值异或a [r+1] 去和 t 做对比,复杂度太大,我们可以反过来更新对比对象 t=r^a[r+1],这样就可以不用更新已经保存过的值了,那么新的 a[r+1] 就不是 a[r+1]了,是什么留给读者思考哦。
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6;
int vis[maxn];
int main()
{
int n,x,t,s=0,ans=0,mx=0;
scanf("%d%d",&n,&t);
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
t^=x;
vis[s]=i;
if(vis[t]>mx)
ans++,mx=i,vis[t]=0;
s^=x;
}
printf("%d\n",ans);
}
思路:把(x ,y ,dir )看成一个状态去bfs搜索就行,可以直接用普通队列去搜索,复杂度n*m*4。
#include<bits/stdc++.h>
#define ll long long
#define db double
using namespace std;
const int maxn=2001;
int vis[maxn][maxn][4],n,m,sx,sy,tx,ty;
int dir[4][2]={1,0,-1,0,0,1,0,-1};
char s[maxn][maxn];
struct node
{
int x,y,w;
node(int a,int b,int c)
{
x=a,y=b,w=c;
}
};
queue<node>q;
int check(int x,int y,int d)
{
if(x<1||x>n||y<1||y>m||s[x][y]=='*'||vis[x][y][d])
return 0;
return 1;
}
int bfs()
{
q.push(node(sx,sy,0));
while(!q.empty())
{
node tmp=q.front();q.pop();
for(int i=0;i<4;i++)
{
int x=tmp.x,y=tmp.y;
x+=dir[i][0];
y+=dir[i][1];
while(check(x,y,i))
{
vis[x][y][i]=1;
if(s[x][y]=='T')
return tmp.w;
q.push(node(x,y,tmp.w+1));
x+=dir[i][0];
y+=dir[i][1];
}
}
}
return -1;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%s",s[i]+1);
for(int j=1;j<=m;j++)
if(s[i][j]=='S')
sx=i,sy=j;
else if(s[i][j]=='T')
tx=i,ty=j;
}
int ans=bfs();
if(ans==-1)
puts("troil");
else printf("%d\n",ans);
}
思路:先离线建树,然后树剖,线段树保存区间或值,再然后重新按顺序查询,操作1,先用并查集判环,如果无环,线段树更新两个点的权值,操作2,直接查询两个点路径的权值或即可。
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+10;
int id[maxn],cnt,top[maxn],p[maxn],a[maxn];
int size[maxn],f[maxn],son[maxn],dep[maxn];
vector<int>G[maxn];
ll sum[maxn*4];
void dfs1(int u,int fa,int deep)
{
size[u]=1;
dep[u]=deep;
f[u]=fa;
for(int i=0;i<G[u].size();i++)
{
int v=G[u][i];
if(v==fa)continue;
dfs1(v,u,deep+1);
size[u]+=size[v];
if(!son[u]||size[v]>size[son[u]])
son[u]=v;
}
}
void dfs2(int u,int root)
{
top[u]=root;
id[u]=++cnt;
if(!son[u])return;
dfs2(son[u],root);
for(int i=0;i<G[u].size();i++)
{
int v=G[u][i];
if(v==f[u]||v==son[u])continue;
dfs2(v,v);
}
}
void update(int o,int l,int r,int k,ll v)
{
if(l==r)
{
sum[o]=v;
return;
}
int ls=o*2,rs=o*2+1,m=(l+r)/2;
if(k<=m)update(ls,l,m,k,v);
else update(rs,m+1,r,k,v);
sum[o]=sum[ls]|sum[rs];
}
ll query(int o,int l,int r,int ql,int qr)
{
if(l>=ql&&r<=qr)
return sum[o];
int ls=o*2,rs=o*2+1,m=(l+r)/2;
ll res=0;
if(ql<=m)res|=query(ls,l,m,ql,qr);
if(qr>m)res|=query(rs,m+1,r,ql,qr);
return res;
}
ll qu_path(int x,int y)
{
ll ans=0;
while(top[x]!=top[y])
{
if(dep[top[x]]<dep[top[y]])swap(x,y);
ans|=query(1,1,cnt,id[top[x]],id[x]);
x=f[top[x]];
}
if(id[x]>id[y])swap(x,y);
ans|=query(1,1,cnt,id[x],id[y]);
return ans;
}
int Count(ll x)
{
int res=0;
while(x)
{
if(x&1)res++;
x/=2;
}
return res;
}
int find(int x)
{
if(p[x]!=x)
p[x]=find(p[x]);
return p[x];
}
int op[maxn*3],x[maxn*3],y[maxn*3];
int main()
{
int n,q,u,v;
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++)
{
p[i]=i;
scanf("%d",&a[i]);
}
for(int i=1;i<=q;i++)
{
scanf("%d%d%d",&op[i],&x[i],&y[i]);
if(op[i]==1)
{
u=find(x[i]),v=find(y[i]);
if(u!=v)
{
p[u]=v;
G[x[i]].push_back(y[i]);
G[y[i]].push_back(x[i]);
}
}
}
dfs1(1,0,1);
dfs2(1,1);
for(int i=1;i<=n;i++)p[i]=i;
for(int i=1;i<=q;i++)
{
if(op[i]==1)
{
u=find(x[i]),v=find(y[i]);
if(u!=v)
{
p[u]=v;
a[x[i]]=a[y[i]]=(a[x[i]]+a[y[i]])/2;
update(1,1,n,id[x[i]],1ll<<a[x[i]]);
update(1,1,n,id[y[i]],1ll<<a[y[i]]);
}
}
else
{
u=find(x[i]);
v=find(y[i]);
if(u!=v)
{
puts("-1");
continue;
}
ll ans=qu_path(x[i],y[i]);
printf("%d\n",Count(ans));
}
}
}