USACO 2021 January Contest, Platinum Problem 1. Sum of Distances 铂金组T1

2024-02-04 14:32

本文主要是介绍USACO 2021 January Contest, Platinum Problem 1. Sum of Distances 铂金组T1,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

USACO 2021 January Contest, Platinum Problem 1. Sum of Distances 铂金组T1

题意概述
给出K张图,设K元组,(1,1,1,1,……,1)出发,(a1,a2,a3,……,ak)与(b1,b2,b3,……,bk)有边当且仅当a1与b1,a2与b2,……,ak与bk均有边连接
求出发可到达的距离和

solution
对于多元组a1-ak,认为是第i张图从1走到ai,假如某几张图比其他图更早到达,则可以通过在某条边上反复横跳来等待其他图到达,所以多元组内奇偶性相同,则可以分别求出奇偶最短路,再用奇偶方案相加再减去同时的方案,现在考虑统计答案,我们枚举最短路的值域,由于最短路的长度不会超过节点数,再枚举每个值有哪几张图,乘上它的出现次数,假如图第一次出现就增加计数,否则就除以上一次出现时乘的数,假如所有图都出现过后统计答案,答案乘以当前最短路值,因为对于长度为i的最短路来说,方案数为其他图小于等于i的方案数

#include<bits/stdc++.h>
#define ll long long
#define N 100005
#define mo 1000000007
using namespace std;
int K,g[N][3],inv[N],add[N];
basic_string<int> w[3][N+5];
int getans(basic_string<int> *w){int cnt=0,res=0,t=1;memset(add,0,sizeof add);for(int i=0;i<N;++i){for(int j:w[i]){if(add[j])t=1ll*t*inv[add[j]]%mo;else++cnt;if(cnt==K){res=(res+1ll*t*i)%mo;}t=1ll*t*(++add[j])%mo;}}return res;
}
int main(){inv[1]=1;for (int i=2;i<=N;++i) {inv[i]=1ll*(mo-mo/i)*inv[mo%i]%mo;}scanf("%d",&K);for(int o=1;o<=K;++o){int n,m;basic_string<int> odd,even;vector<basic_string<int> > t;queue< pair<int,bool> > que;scanf("%d%d",&n,&m);odd.assign(n+1,N),even.assign(n+1,N),t.resize(n+1);for(int i=1,u,v;i<=m;++i){scanf("%d%d",&u,&v);t[u]+=v,t[v]+=u;}que.push({1,0});even[1]=0;while(!que.empty()){pair<int,bool> i=que.front();que.pop();if(i.second){for(int j:t[i.first]){if(even[j]==N)even[j]=odd[i.first]+1,que.push({j,0});}}elsefor(int j:t[i.first]){if(odd[j]==N)odd[j]=even[i.first]+1,que.push({j,1});}}for(int i=1;i<=n;++i){w[0][odd[i]]+=o,w[1][even[i]]+=o,w[2][max(even[i],odd[i])]+=o;}}printf("%d",((getans(w[0])+getans(w[1])-getans(w[2]))%mo+mo)%mo);
}

这篇关于USACO 2021 January Contest, Platinum Problem 1. Sum of Distances 铂金组T1的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/677818

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

usaco 1.3 Calf Flac(暴搜)

思路是暴搜。 需要注意的地方是输入的方法,以及输出时的换行。 代码: /*ID: who jayLANG: C++TASK: calfflac*/#include<stdio.h>#include<string.h>#include<math.h>int main(){freopen("calfflac.in","r",stdin);freopen("calfflac.ou

usaco 1.3 Barn Repair(贪心)

思路:用上M块木板时有 M-1 个间隙。目标是让总间隙最大。将相邻两个有牛的牛棚之间间隔的牛棚数排序,选取最大的M-1个作为间隙,其余地方用木板盖住。 做法: 1.若,板(M) 的数目大于或等于 牛棚中有牛的数目(C),则 目测 给每个牛牛发一个板就为最小的需求~ 2.否则,先对 牛牛们的门牌号排序,然后 用一个数组 blank[ ] 记录两门牌号之间的距离,然后 用数组 an

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

usaco 1.2 Palindromic Squares(进制转化)

考察进制转化 注意一些细节就可以了 直接上代码: /*ID: who jayLANG: C++TASK: palsquare*/#include<stdio.h>int x[20],xlen,y[20],ylen,B;void change(int n){int m;m=n;xlen=0;while(m){x[++xlen]=m%B;m/=B;}m=n*n;ylen=0;whi

usaco 1.2 Name That Number(数字字母转化)

巧妙的利用code[b[0]-'A'] 将字符ABC...Z转换为数字 需要注意的是重新开一个数组 c [ ] 存储字符串 应人为的在末尾附上 ‘ \ 0 ’ 详见代码: /*ID: who jayLANG: C++TASK: namenum*/#include<stdio.h>#include<string.h>int main(){FILE *fin = fopen (

usaco 1.2 Milking Cows(类hash表)

第一种思路被卡了时间 到第二种思路的时候就觉得第一种思路太坑爹了 代码又长又臭还超时!! 第一种思路:我不知道为什么最后一组数据会被卡 超时超了0.2s左右 大概想法是 快排加一个遍历 先将开始时间按升序排好 然后开始遍历比较 1 若 下一个开始beg[i] 小于 tem_end 则说明本组数据与上组数据是在连续的一个区间 取max( ed[i],tem_end ) 2 反之 这个

usaco 1.2 Transformations(模拟)

我的做法就是一个一个情况枚举出来 注意计算公式: ( 变换后的矩阵记为C) 顺时针旋转90°:C[i] [j]=A[n-j-1] [i] (旋转180°和270° 可以多转几个九十度来推) 对称:C[i] [n-j-1]=A[i] [j] 代码有点长 。。。 /*ID: who jayLANG: C++TASK: transform*/#include<

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b