也是CF985G。。。

容斥+三元环计数

CF数据太弱啦

vis没赋初值-1竟然过了QAQ

所以又调了我半个小时才搞掉QAQ

数数真难QAQ

记得要写#include<vector>!!!

Dev给加的奇奇怪怪的编译选项会给你自动填补的!!!

QAQ长个记性QAQ

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<vector>
#define ull unsigned long long
#define ll long long
#define inf 20021225
#define N 400010
#define it vector<int>::iterator
using namespace std; int sz[N],deg[N],n,m;
vector<ull> sum[N];
vector<int> ed[N];
struct edge{int to,lt;}e[N];
int in[N],cnt; int vis[N];
void add(int x,int y){e[++cnt].to=y;e[cnt].lt=in[x];in[x]=cnt;}
void go(int x){for(int i=in[x];i;i=e[i].lt) vis[e[i].to]=x;}
ull calc(int x){if(x<) return ; return (ll)x*(x-)/;}
ull up(int l,int r){if(l>r) return ; return (ll)(r+l)*(r-l+)/;}
int main()
{
scanf("%d%d",&n,&m); ull A,B,C,res=;
scanf("%llu%llu%llu",&A,&B,&C); int x,y;
//scanf("%I64d%I64d%I64d",&A,&B,&C); int x,y;
for(int i=;i<=m;i++)
{
scanf("%d%d",&x,&y);
ed[x].push_back(y);
ed[y].push_back(x);
}
//ull edg=(ull)n*(n-1)>>1;
//at least 0
for(int i=;i<n;i++)
res+=A*i*calc(n-i-),
res+=B*i*i*(n-i-),
res+=C*i*calc(i),
sort(ed[i].begin(),ed[i].end()),
sum[i].push_back();
//printf("%llu\n",res);
for(x=;x<n;x++) for(it i=ed[x].begin();i!=ed[x].end();i++)
sz[x]=(*i)>x?sz[x]:sz[x]+,++deg[x],sum[x].push_back(sum[x][deg[x]-]+(*i));
//at least 1
for(x=;x<n;x++) for(it i=ed[x].begin();i!=ed[x].end();i++)
if((*i)<x)
{
y=*i;
//x as B
res-=B*x*(n-x-);
res-=A*y*(n-x-);
res-=C*up(x+,n-);
//x as C
res-=A*up(,y-);
res-=B*y*y;
res-=B*up(y+,x-);
res-=A*y*(x-y->?x-y-:);
res-=C*x*(x-);
}
//printf("%llu\n",res);
//at least 2
//int dig;
//printf("%llu\n",res);
for(x=;x<n;x++) for(int i=;i<=deg[x];i++)
{
//printf("%d %llu\n",x,sum[x][i]);
//x has 2 edges out
y=ed[x][i-];
if(y>x)
{
//x is A
res+=A*x*(deg[x]-i); res+=B*y*(deg[x]-i);
res+=C*(sum[x][deg[x]]-sum[x][i]);
//printf("%d %d %llu\n",x,y,res);
}
else// if(y<x)
{
//x is B
res+=B*x*(deg[x]-sz[x]); res+=A*y*(deg[x]-sz[x]);
res+=C*(sum[x][deg[x]]-sum[x][sz[x]]);
//x is C
res+=C*x*(i-); res+=B*y*(i-);
res+=A*sum[x][i-];
}
}
//printf("%llu\n",res);
//at least 3
/**for(int i=1;i<=n;i++)
{
if(deg[i]<=top)
{
for(int j=0;j<deg[i];j++)
{
for(int k=j+1;k<deg[i];k++)
{
x=i; y=ed[x][j]; z=ed[x][k];
if(x>y) swap(x,y);
if(x>z) swap(x,z);
res-=A*x+B*y+C*z;
}
}
}
else big.push_back(i);
}*/
for(x=;x<n;x++) for(int i=;i<deg[x];i++)
if(deg[ed[x][i]]<deg[x] || (deg[ed[x][i]]==deg[x]&&ed[x][i]<x))
add(x,ed[x][i]);
memset(vis,-,sizeof(vis));
for(x=;x<n;x++)
{
go(x);
for(int i=in[x];i;i=e[i].lt)
{
y=e[i].to;
for(int j=in[y];j;j=e[j].lt)
{
if(vis[e[j].to]==x)
{
int xx=x,yy=y,zz=e[j].to;
if(xx>yy) swap(xx,yy);
if(xx>zz) swap(xx,zz);
if(yy>zz) swap(yy,zz);
res-=A*xx+B*yy+C*zz;
}
}
}
}
printf("%llu\n",res);
//printf("%I64u\n",res);
return ;
}
/**
4 1
2 3 4
1 0
*/

由于这个题调的我太过痛苦 写作业都写不进去

bug列举一下

1.我不会冒泡排序排三个数

2.(0,n-1)看成(1,n)重调系数非常蛋疼

3.初值就是忘了赋

4.系数和值少一个

5.手贱

最新文章

  1. Python: 利用Python进行数据分析 学习记录
  2. laravel中TokenMismatchException异常处理
  3. 旷世奇坑!!!spring 不能自动注入
  4. [数据库]漫谈ElasticSearch关于ES性能调优几件必须知道的事(转)
  5. sql 查询最近30分钟或者自定义时间数据
  6. QT pro文件解析
  7. svo:一个半直接单目视觉里程计(长期更新)
  8. Ubuntu 14.4 使用中遇到的问题汇总
  9. ListView添加节点
  10. HDU 2092 整数解
  11. 读书笔记之《Java编程思想》
  12. js对象的合并
  13. JVM学习一:JVM之类加载器概况
  14. Java访问级别修饰符
  15. 使用 JavaScript 将网站后台的数据变化实时更新到前端
  16. SmtpClient SSL 发送邮件异常排查
  17. tkinter widget
  18. Wannafly挑战赛9 E - 组一组
  19. USACO 4.3 Street Race
  20. 机器学习实战:用nodejs实现人脸识别

热门文章

  1. swan.after
  2. hdu 4641K-string SAM的O(n^2)算法 以及 SAM+并查集优化
  3. Python_011(生成器)
  4. 解决webpack打包vue项目后,部署完成后,刷新页面页面404
  5. Java机试题目
  6. grammar_action
  7. Myeclipse优化配置
  8. Vagrant 手册之 Multi-machine 多机器
  9. memcpy复制字符串的注意事项/memcpy不能用来拷贝类类型
  10. hive环境