题面

传送门(bzoj)

传送门(CF)

\(llx\)身边妹子成群,这天他需要从\(n\)个妹子中挑出\(3\)个出去浪,但是妹子之间会有冲突,表现为\(i,j\)之间连有一条边\((i,j)\),定义一种选择方案的权值为\(Ai+Bj+Ck,i<j<k\),求所有选择方案的权值之和

题解

容斥,至少\(0\)条边相连的方案\(-\)至少\(1\)条边相连的方案\(+\)至少\(2\)条边相连的方案\(-\)至少\(3\)条边相连的方案

至少\(3\)条边相连的方案最难数,是个三元环计数,和这题方法一样

虽然我感觉自己去数啥都数不出来

//minamoto
#include<bits/stdc++.h>
#define R register
#define ull unsigned long long
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
const int N=2e5+5;
struct eg{int v,nx;}e[N<<1];int head[N],tot;
inline void add(R int u,R int v){e[++tot]={v,head[u]},head[u]=tot;}
struct EG{int u,v;}E[N];
vector<int>to[N],Pre[N],sum[N];int deg[N],sz[N],ss[N],vis[N];
int n,m,tim;ull res,A,B,C,ps[N];
inline ull bin(R int n){return 1ll*n*(n-1)>>1;}
inline ull calc(R int l,R int r){return 1ll*(l+r)*(r-l+1)>>1;}
int main(){
// freopen("testdata.in","r",stdin);
n=read(),m=read(),A=read(),B=read(),C=read();
for(R int i=1,u,v;i<=m;++i){
u=read(),v=read(),u>v?swap(u,v),0:0;
++sz[u],++ss[v],ps[v]+=u,++deg[u],++deg[v],to[u].push_back(v),Pre[v].push_back(u);
E[i].u=u,E[i].v=v;
}
fp(i,0,n-1)sort(to[i].begin(),to[i].end()),sort(Pre[i].begin(),Pre[i].end());
fp(i,1,m)deg[E[i].u]>deg[E[i].v]||(deg[E[i].u]==deg[E[i].v]&&E[i].u>E[i].v)?add(E[i].u,E[i].v):add(E[i].v,E[i].u);
fp(i,0,n-3)res+=A*i*bin(n-1-i);
fp(i,1,n-2)res+=B*i*i*(n-1-i);
fp(i,2,n-1)res+=C*i*bin(i);
fp(i,0,n-1){
sum[i].resize(sz[i]);sz[i]?sum[i][0]=to[i][0]:0;
fp(j,1,sz[i]-1)sum[i][j]=sum[i][j-1]+to[i][j];
}
for(R int i=1,u,v;i<=m;++i){
u=E[i].u,v=E[i].v;
res-=(A*u+B*v)*(n-1-v)+C*calc(v+1,n-1),
res-=(B*u+C*v)*u+A*calc(0,u-1),
res-=(A*u+C*v)*(v-u-1)+B*calc(u+1,v-1);
}
fp(u,0,n-1)fp(i,0,sz[u]-1){
R int v=to[u][i];
res+=(A*u+B*v)*(sz[u]-i-1)+C*(sum[u][sz[u]-1]-sum[u][i]);
if(sz[v])res+=(A*u+B*v)*sz[v]+C*sum[v][sz[v]-1];
}
for(R ull i=0,sum=0;i<n;++i,sum=0)
fp(j,0,ss[i]-1){
R int u=Pre[i][j];
res+=(B*u+C*i)*j+A*sum,
sum+=Pre[i][j];
}
fp(u,0,n-1){
++tim;
go(u)vis[v]=tim;
for(R int k=head[u];k;k=e[k].nx)go(e[k].v)if(vis[v]==tim){
R int x=u,y=e[k].v,z=v;
if(x>y)swap(x,y);if(x>z)swap(x,z);if(y>z)swap(y,z);
res-=A*x+B*y+C*z;
}
}
printf("%llu\n",res);
//printf("%I64u\n",res); 要交到CF上的话得写这个格式
return 0;
}

BZOJ.5407.girls/CF985G. Team Players(三元环计数+容斥)的更多相关文章

  1. [CF985G]Team Players

    题意:给出一个图,求$\sum\limits_{\substack{i\lt j\lt k\\\nexists(i,j),(j,k),(i,k)}}Ai+Bj+Ck$ 挺好的一道题==,就是稍微毒了点 ...

  2. 【BZOJ5332】[SDOI2018]旧试题(数论,三元环计数)

    [BZOJ5332][SDOI2018]旧试题(数论,三元环计数) 题面 BZOJ 洛谷 题解 如果只有一个\(\sum\),那么我们可以枚举每个答案的出现次数. 首先约数个数这个东西很不爽,就搞一搞 ...

  3. Codechef SUMCUBE Sum of Cubes 组合、三元环计数

    传送门 好久没有做过图论题了-- 考虑\(k\)次方的组合意义,实际上,要求的所有方案中导出子图边数的\(k\)次方,等价于有顺序地选出其中\(k\)条边,计算它们在哪一些图中出现过,将所有方案计算出 ...

  4. loj#6076「2017 山东一轮集训 Day6」三元组 莫比乌斯反演 + 三元环计数

    题目大意: 给定\(a, b, c\),求\(\sum \limits_{i = 1}^a \sum \limits_{j = 1}^b \sum \limits_{k = 1}^c [(i, j) ...

  5. [hdu 6184 Counting Stars(三元环计数)

    hdu 6184 Counting Stars(三元环计数) 题意: 给一张n个点m条边的无向图,问有多少个\(A-structure\) 其中\(A-structure\)满足\(V=(A,B,C, ...

  6. LOJ2565 SDOI2018 旧试题 莫比乌斯反演、三元环计数

    传送门 这道题的思路似乎可以给很多同时枚举三个量的反演题目提供一个很好的启发-- 首先有结论:\(d(ijk) = \sum\limits_{x|i}\sum\limits_{y|j}\sum\lim ...

  7. BZOJ.5407.girls(容斥 三元环)

    题目链接 CF 原题 \(Description\) 有n个点,其中有m条边连接两个点.每一个没有连边的三元组\((i,j,k)(i<j<k)\)对答案的贡献为\(A*i+B*j+C*k\ ...

  8. HDU 6184 Counting Stars 经典三元环计数

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6184 题意: n个点m条边的无向图,问有多少个A-structure 其中A-structure满足V ...

  9. HDU6184【Counting Stars】(三元环计数)

    题面 传送门 给出一张无向图,求 \(4\) 个点构成两个有公共边的三元环的方案数. 题解 orz余奶奶,orz zzk 首先,如果我们知道经过每条边的三元环个数\(cnt_i\),那么答案就是\(\ ...

随机推荐

  1. Simple Cubemap Reflection

    [Simple Cubemap Reflection] Cubemap加在MainTex上,所以Property需要按如下定义: 注意_Cubemap的类型是CUBE. 使用Cubemap,需要计算反 ...

  2. python 文件的读取&更新

    [python 文件的读取&更新] 任务抽象: 读取一个文件, 更新内容后, 重新写入文件. 实际应用: 磁盘上的一个配置文件, 读入内存后为一个dict, 对dict更新后重新写入磁盘. d ...

  3. 【UVALive2965】Jurassic Remains

    题意 题意给定n个大写字母组成的字符串,选择尽量多的串,使得每个大写字母都能出现偶数次.n<=24,每个字符串中每个字母最多出现一次. 分析 这是训练指南上的一道中途相遇法的简单题,但是好像也能 ...

  4. 117. Populating Next Right Pointers in Each Node II 计算右边的附属节点

    [抄题]: Given a binary tree struct TreeLinkNode { TreeLinkNode *left; TreeLinkNode *right; TreeLinkNod ...

  5. java校验银行卡号

    public class CheckBankCard { /* 校验过程: 1.从卡号最后一位数字开始,逆向将奇数位(1.3.5等等)相加. 2.从卡号最后一位数字开始,逆向将偶数位数字,先乘以2(如 ...

  6. log4j-over-slf4j工作原理详解

    log4j-over-slf4j工作原理详解 摘自:https://blog.csdn.net/john1337/article/details/76152906 置顶 2017年07月26日 17: ...

  7. CENTOS7 YUM安装BOOST1.53(静态版本)

    按照之前的博文更新163的源之后,执行: yum install boost-static.i686 yum install boost-devel.i686 yum install boost-do ...

  8. NSPredicate过滤数组数据

    NSPredicate编写软件时,经常需要获取一个对象集合,然后删除不满足条件的对象,保留符合条件的对象,从而提供一些有意义的对象.Cocoa提供了一个名为NSPredicate的类,他用于指定过滤器 ...

  9. 关于Qt官方下载页的最新变动

    时间过得很快,现在Qt已经迎来了5.10版本,但是当我们去下载页下载对应安装包的时候,已经找不到之前的offline安装包了.你能够看到的只有在线安装包,并且我自己有做过测试,国内的网络基本上没有机会 ...

  10. Gym - 100989H (贪心)

    After the data structures exam, students lined up in the cafeteria to have a drink and chat about ho ...