题目链接:https://jzoj.net/senior/#main/show/6099

考虑直接统计某个点到其它所有点的距离和

我们先把整个团当成一个点建图,处理出任意两个团之间的距离\(dis(i,j)\),注意这里的\(dis\)表示的是两个不位于团的相交部分的两点之间距离,即路程一定是某个点\(a->相交部分->b->……\)这样的

接下来我们枚举点\(x\),统计所有点到\(x\)的最短路,我们考虑从\(x\)到所有团的最短路,并按照该最短路从小到大的顺序考虑它们对答案的贡献

对于\(x\)和任意一点\(y\),其最短路一定是由团之间的距离产生的答案,我们假设\(x\)和\(y\)分别属于联通块\(i,j\),那么\(dis(i,j)\)想要对答案产生贡献的前提条件就是在集合\(j\)中所有比其小的\(dis\)仍未未出现

即它对答案的贡献就是那些属于当前集合\(j\)中的除了已经计算贡献的点以外的点数

我们记\(cnt[i][sta]\)表示团\(i\)中的点一定存在,\(sta\)中为0的点一定不存在的点的数量

接下来暴力统计答案即可

#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
#include<queue>
#include<bitset>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x)&(-x)
#define rep(i,a,b) for (int i=a;i<=b;i++)
#define per(i,a,b) for (int i=a;i>=b;i--)
#define maxd 1000000007
typedef long long ll;
const int N=100000;
const double pi=acos(-1.0);
int n,m,cost[20],in[100100],id[20];
ll dis[20][20],cnt[20][1<<20],all,sum[1<<20],d[20];
bitset<100100> block[20]; 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*10+(ch-'0');ch=getchar();}
return x*f;
} bool cmp(int x,int y)
{
return d[x]<d[y];
} void init()
{
n=read();m=read();all=(1<<m)-1;
rep(i,1,m)
{
cost[i]=read();int len=read();
rep(j,1,len)
{
int x=read();
in[x]|=(1<<(i-1));
block[i][x]=1;
}
}
rep(i,1,n) sum[in[i]]++;
rep(i,1,m)
{
rep(j,0,all) if (j&(1<<(i-1))) cnt[i][j]=sum[j];
rep(j,1,m)
{
rep(k,0,all)
{
if (k&(1<<(j-1))) cnt[i][k]+=cnt[i][k^(1<<(j-1))];
}
}
}
} void floyd()
{
rep(i,1,m) rep(j,1,m) dis[i][j]=dis[j][i]=maxd;
rep(i,1,m) dis[i][i]=cost[i];
rep(i,1,m)
{
rep(j,1,m)
{
if (i==j) continue;
if ((block[i]&block[j]).any()) {dis[i][j]=cost[i]+cost[j];dis[j][i]=dis[i][j];}
}
}
rep(k,1,m)
{
rep(i,1,m)
{
rep(j,1,m)
{
if ((k==i) || (k==j) || (i==j)) continue;
ll nowdis=dis[i][k]+dis[k][j];
if ((k!=i) && (k!=j)) nowdis-=cost[k];
dis[i][j]=min(nowdis,dis[i][j]);
}
}
}
} void work()
{
ll ans=0;
rep(i,1,n)
{
rep(j,1,m) {id[j]=j;d[j]=maxd;}
rep(j,1,m)
{
if (in[i]&(1<<(j-1)))
{
rep(k,1,m) d[k]=min(d[k],dis[j][k]);
}
}
sort(id+1,id+1+m,cmp);
int sta=0;
rep(j,1,m)
{
ans+=d[id[j]]*cnt[id[j]][all^sta];
sta|=(1<<(id[j]-1));
}
ans-=d[id[1]];
}
printf("%lld",ans/2);
} void file()
{
freopen("dist.in","r",stdin);
freopen("dist.out","w",stdout);
} int main()
{
file();
init();
floyd();
work();
return 0;
}

jzoj6099. 【GDOI2019模拟2019.4.1】Dist的更多相关文章

  1. jzoj6101. 【GDOI2019模拟2019.4.2】Path

    题目链接:https://jzoj.net/senior/#main/show/6101 记\(f_i\)为从\(i\)号点走到\(n\)号点所花天数的期望 那么根据\(m\)条边等可能的出现一条和一 ...

  2. [JZOJ6075]【GDOI2019模拟2019.3.20】桥【DP】【线段树】

    Description N,M<=100000,S,T<=1e9 Solution 首先可以感受一下,我们把街道看成一行,那么只有给出的2n个点的纵坐标是有用的,于是我们可以将坐标离散化至 ...

  3. [jzoj 6101] [GDOI2019模拟2019.4.2] Path 解题报告 (期望)

    题目链接: https://jzoj.net/senior/#main/show/6101 题目: 题解: 设$f_i$表示从节点$i$到节点$n$的期望时间,$f_n=0$ 最优策略就是如果从$i, ...

  4. [jzoj 6093] [GDOI2019模拟2019.3.30] 星辰大海 解题报告 (半平面交)

    题目链接: https://jzoj.net/senior/#contest/show/2686/2 题目: 题解: 说实话这题调试差不多花了我十小时,不过总算借着这道题大概了解了计算几何的基础知识 ...

  5. [jzoj 6080] [GDOI2019模拟2019.3.23] IOer 解题报告 (数学构造)

    题目链接: https://jzoj.net/senior/#main/show/6080 题目: 题意: 给定$n,m,u,v$ 设$t_i=ui+v$ 求$\sum_{k_1+k_2+...+k_ ...

  6. [jzoj 6092] [GDOI2019模拟2019.3.30] 附耳而至 解题报告 (平面图转对偶图+最小割)

    题目链接: https://jzoj.net/senior/#main/show/6092 题目: 知识点--平面图转对偶图 在求最小割的时候,我们可以把平面图转为对偶图,用最短路来求最小割,这样会比 ...

  7. [jzoj 6086] [GDOI2019模拟2019.3.26] 动态半平面交 解题报告 (set+线段树)

    题目链接: https://jzoj.net/senior/#main/show/6086 题目: 题解: 一群数字的最小公倍数就是对它们质因数集合中的每个质因数的指数取$max$然后相乘 这样的子树 ...

  8. [jzoj 4528] [GDOI2019模拟2019.3.26] 要换换名字 (最大权闭合子图)

    题目链接: https://jzoj.net/senior/#contest/show/2683/0 题目: 题解: 不妨枚举一个点,让两颗树都以这个点为根,求联通块要么点数为$0$,要么包括根(即联 ...

  9. [jzoj 6087] [GDOI2019模拟2019.3.26] 获取名额 解题报告 (泰勒展开+RMQ+精度)

    题目链接: https://jzoj.net/senior/#main/show/6087 题目: 题解: 只需要统计$\prod_{i=l}^r (1-\frac{a_i}{x})$ =$exp(\ ...

随机推荐

  1. arcgis api 3.x for js入门开发系列九热力图效果(附源码下载)

    前言 关于本篇功能实现用到的 api 涉及类看不懂的,请参照 esri 官网的 arcgis api 3.x for js:esri 官网 api,里面详细的介绍 arcgis api 3.x 各个类 ...

  2. WebLogic使用总结

    一. WebLogic安装 WebLogic使用总结(一)——WebLogic安装 一.下载WebLogic 到Oracle官网http://www.oracle.com/ 下载WebLogic(根据 ...

  3. jupyter notebook安装、登录

    pip install jupyter 提示pip需要升级(本人装的是anaconda) 输入:python -m pip install --upgrade pip 安装完成. 运行jupyter ...

  4. android找不到aar包

    转载请标明出处,维权必究:https://www.cnblogs.com/tangZH/p/9939663.html  在做项目的时候引入aar包,编译的时候却提示错误(这个错误大概说的是...... ...

  5. vue父子组件之间传值

    vue父子组件进行传值 vue中的父子组件,什么是父组件什么是子组件呢?就跟html标签一样,谁包裹着谁谁就是父组件,被包裹的元素就是子组件. 父组件向子组件传值 下面用的script引入的方式,那种 ...

  6. redis快照持久化和aof日志持久化

    持久化就是即使断电/重启需要存储的数据不会丢失,即将数据存储在设备中,一般存在硬盘内 redis的持久化有2种方式 :1-rdb快照  2-aof日志,可以通过配置redis.conf文件进行配置 r ...

  7. 【笔记】基于Python的数字图像处理

    [博客导航] [Python相关] 前言 基于Python的数字图像处理,离不开相关处理的第三方库函数.搜索网络资源,列出如下资源链接. Python图像处理库到底用哪家 python计算机视觉编程— ...

  8. Python操作db2

    官方文档:https://www.ibm.com/support/knowledgecenter/en/SSEPGG_9.5.0/com.ibm.db2.luw.apdv.python.doc/doc ...

  9. Django 【认证系统】auth

    本篇内容 介绍Django框架提供的auth 认证系统 方法: 方法名 备注 create_user 创建用户 authenticate 登录验证 login 记录登录状态 logout 退出用户登录 ...

  10. TCP/IP及内核参数优化调优

    Linux下TCP/IP及内核参数优化有多种方式,参数配置得当可以大大提高系统的性能,也可以根据特定场景进行专门的优化,如TIME_WAIT过高,DDOS攻击等等.如下配置是写在sysctl.conf ...