Codeforces Round #411 (Div. 2)

本文提供了CodeForces 805比赛的A-F题详细解答过程及代码实现,涵盖贪心算法、字符串处理、图论、树结构等核心算法知识点。

拿游子小号和大家打了一场codeforces,感觉自己被肝爆了qwqqq,wxl求轻虐...qwq



A

貌似是让求区间[l,r]中出现次数最多的约数...那就是2呗...当然l==r时是l

//codeforces 805A
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
int n,l,r;
int main()
{
 int i,j;
 scanf("%d%d",&l,&r);
 if(l==r) printf("%d\n",l);
 else printf("2\n");
 return 0;
}




B

构造一个ab序列,不能出现长度为3的回文

那就aabbaabb无限循环...

//codeforces 805B
#include<cstdio>
#include<cstring>//codeforces 805D
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mod=1e9+7;
const int mxn=1000005;
ll n,num,cnt,res,ans;
ll pw[mxn];
char s[mxn];
int main()
{
 int i,j;
 scanf("%s",s+1);
 n=strlen(s+1);
 pw[0]=1;
 fo(i,1,n) pw[i]=(pw[i-1]*2)%mod;
 fo(i,1,n)
 {
  if(s[i]=='a') cnt++;
  else if(s[i]=='b' && (s[i+1]=='a' || i==n))
  {
   res++;
   ans=(ans+res*(pw[cnt]-1)%mod)%mod;
   res=0;
  }
  else res++;
 }
 printf("%lld\n",ans);
 return 0;
}




C

贪心的来吧...1->n->2->n-1->3->n-2...

大力观察得出这样是最优解→_→



//codeforces 805C
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
int n;
int main()
{
 int i,j;
 scanf("%d",&n);
 if(n&1) printf("%d\n",n/2);
 else printf("%d\n",n/2-1);
 return 0;
}


D

给一个ab序列,如果出现ab就变成bba,问变几次就不需要再变了...

大力观察么...

//codeforces 805D
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mod=1e9+7;
const int mxn=1000005;
ll n,num,cnt,res,ans;
ll pw[mxn];
char s[mxn];
int main()
{
 int i,j;
 scanf("%s",s+1);
 n=strlen(s+1);
 pw[0]=1;
 fo(i,1,n) pw[i]=(pw[i-1]*2)%mod;
 fo(i,1,n)
 {
  if(s[i]=='a') cnt++;
  else if(s[i]=='b' && (s[i+1]=='a' || i==n))
  {
   res++;
   ans=(ans+res*(pw[cnt]-1)%mod)%mod;
   res=0;
  }
  else res++;
 }
 printf("%lld\n",ans);
 return 0;
}


E

这个题...观察到联通块的性质,所以暴力dfs+贪心染色就好啦

//codeforces 805E. Ice cream coloring
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=300005;
bool bo[mxn],wxl[mxn];
int n,m,cnt,tot,ans=1;
vector <int> e[mxn];
int head[mxn],color[mxn];
struct node {int to,next;} f[mxn<<1];

inline void add(int u,int v)
{
 f[++cnt].to=v,f[cnt].next=head[u],head[u]=cnt;
}

inline void dfs(int u,int fa)
{
 int tmp=1,v;
 for(int i=0;i<e[u].size();i++)
 {
  v=e[u][i];
  bo[color[v]]=1;
 }
 for(int i=0;i<e[u].size();i++)
 {
  v=e[u][i];
  while(bo[tmp]) tmp++;
  if(!color[v]) color[v]=tmp++;
 }
 for(int i=0;i<e[u].size();i++)
 {
  v=e[u][i];
  bo[color[v]]=0;
 }
 for(int i=head[u];i;i=f[i].next)
 {
  v=f[i].to;
  if(v==fa) continue;
  dfs(v,u);
 }
}

int main()
{
 int i,j,x,u,v;
 scanf("%d%d",&n,&m);
 fo(i,1,n)
 {
  scanf("%d",&x);
  while(x--)
  {
   scanf("%d",&v);
   e[i].push_back(v);
  }
 }
 fo(i,2,n)
 {
  scanf("%d%d",&u,&v);
  add(u,v),add(v,u);
 }
 dfs(1,0);
 fo(i,1,m) ans=max(ans,color[i]);
 printf("%d\n",ans);
 fo(i,1,m) printf("%d ",max(color[i],1));
 printf("\n");
 return 0;
}




F

调了一上午嗯嗯嗯

这题超讨厌...树的题完全不会...真的好难啊qwqqq

做出来还是不会证明复杂度qwqqq

总的来说:map大法好,vector把你保


//codeforces 805F
#include<map>
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=100005;
double x,ans;
vector <int> e[mxn],res[mxn];
map < pair<int,int>,double > mp;
bool vis[mxn],dia[mxn];
int n,m,q,fff,cnt,tot,mx;
int head[mxn],p[mxn][3],len[mxn],last[mxn];
int father[mxn],pre[mxn],dist[mxn],size[mxn];
struct node {int to,next,from;} f[mxn<<1];
inline void add(int u,int v)
{
 f[++cnt].to=v,f[cnt].next=head[u],head[u]=cnt;
}
inline int find(int x)
{
 if(x==father[x]) return father[x];
 return father[x]=find(father[x]);
}
inline void dfs1(int u,int lenth)
{
 vis[u]=1,size[fff]++;
 if(lenth>=mx) mx=lenth,p[fff][1]=u;
 for(int i=head[u];i;i=f[i].next)
 {
  int v=f[i].to;
  if(!vis[v]) dfs1(v,lenth+1);
 }
}
inline void dfs2(int u,int lenth)
{
 vis[u]=1;
 if(lenth>=len[fff]) p[fff][2]=u;
 len[fff]=max(len[fff],lenth);
 for(int i=head[u];i;i=f[i].next)
 {
  int v=f[i].to;
  if(!vis[v]) dfs2(v,lenth+1);
 }
}
inline void dfs3(int u)  //标记最长链
{
 vis[u]=1;
 if(u==p[fff][2])
 {
  dia[u]=1,dist[u]=len[fff];
  int num=0;
  while(pre[u])
    num++,u=pre[u],dia[u]=1,dist[u]=max(num,len[fff]-num);
  return;
 }
 for(int i=head[u];i;i=f[i].next)
 {
  int v=f[i].to;
  if(v==pre[u]) continue;
  pre[v]=u;
  dfs3(v);
 }
}
inline void dfs4(int u,int lenth)
{
 vis[u]=1,dist[u]=lenth;
 e[fff].push_back(lenth);
 for(int i=head[u];i;i=f[i].next)
 {
  int v=f[i].to;
  if(dia[v] || vis[v]) continue;
  dfs4(v,lenth+1);
 }
}
inline double solve(int t1,int t2)
{
 int i,j,l,r;
 ll sum=0,tmp=max(len[t1],len[t2]);
 for(l=0;l<(int)e[t1].size();l++)
 {
  r=upper_bound(e[t2].begin(),e[t2].end(),tmp-e[t1][l]-1)-e[t2].begin();
  sum+=r*tmp+res[t2][size[t2]-1]+(size[t2]-r)*(e[t1][l]+1);
  if(r) sum-=res[t2][r-1];
 }
 return (double)sum/((double)size[t1]*size[t2]);
}
int main()
{
 int i,j,u,v,t1,t2,z;
 scanf("%d%d%d",&n,&m,&q);
 fo(i,1,n) father[i]=i;
 fo(i,1,m)
 {
  scanf("%d%d",&u,&v);
  add(u,v),add(v,u);
  int t1=find(u),t2=find(v);
  if(t1!=t2) father[t1]=t2;
 }
 fo(i,1,n) if(!vis[i])
   fff=find(i),mx=0,dfs1(i,0);
 memset(vis,0,sizeof vis);
 fo(i,1,n) if(!vis[i]) fff=find(i),dfs2(p[fff][1],0);
 memset(vis,0,sizeof vis);
 fo(i,1,n) if(!vis[i]) fff=find(i),dfs3(p[fff][1]);
 memset(vis,0,sizeof vis);
 fo(i,1,n) if(!vis[i] && dia[i]) fff=find(i),dfs4(i,dist[i]);
 fo(fff,1,n) if(father[fff]==fff)
 {
  sort(e[fff].begin(),e[fff].end());
  res[fff].push_back(e[fff][0]);
  for(j=1;j<e[fff].size();j++)
    res[fff].push_back(res[fff][j-1]+e[fff][j]);
 }
 while(q--)
 {
  scanf("%d%d",&u,&v);
  int t1=find(u),t2=find(v);
  if(size[t2]<size[t1]) swap(t1,t2);
  if(t1==t2)
  {
   printf("-1\n");
   continue;
  }
  if(mp.count(make_pair(t1,t2)))
  {
   printf("%.10lf\n",mp[make_pair(t1,t2)]);
   continue;
  }
  ans=mp[make_pair(t1,t2)]=solve(t1,t2);
  printf("%.10lf\n",ans);
 }
 return 0;
}
/*
5 2 1
2 4
4 3
2 5
*/


内容概要:本文围绕“考虑电动汽车聚合可调节能力的含波动性电源电氢耦合系统多目标优化运行”展开研究,提出了一种基于Matlab代码实现的多目标优化模型。该模型深度融合电-氢耦合系统与高比例波动性可再生能源(如风电、光伏),充分挖掘电动汽车(EV)集群作为移动储能单元的灵活调节潜力,通过聚合调控提升系统对新能源的消纳能力与运行经济性。研究系统构建了电动汽车可调度能力、电解水制氢与储氢动态过程、多能源协同互补的优化调度框架,并结合智能优化算法实现经济性、低碳性与运行稳定性等多重目标的协同优化。文中配套提供了完整的Matlab仿真代码、相关数据及可能的论文支撑材料,极大地方便了模型的复现、验证与后续深化研究。; 适合人群:具备电力系统、综合能源系统、优化理论或新能源技术等相关领域基础知识的研究生、科研人员,以及从事新型电力系统规划、清洁能源消纳与智慧能源管理的工程技术人员。; 使用场景及目标:①开展高渗透率可再生能源接入下的综合能源系统多目标优化调度研究;②探究电动汽车集群在电网削峰填谷、平抑新能源出力波动及提供辅助服务方面的应用价值与潜力;③学习并掌握电氢耦合系统的建模方法、多目标优化求解技术及其在Matlab/Simulink环境下的仿真实现流程。; 阅读建议:此资源不仅提供可运行的代码,更蕴含了前沿的科研思路与创新方法,建议读者结合所提供的代码、数据与可能的论文文档,系统性地学习从问题建模、算法设计到仿真分析的完整科研过程,并重点关注其中关于需求侧资源聚合、多能互补协同与绿色低碳运行的核心理念。
内容概要:本文档名为《经济学期刊论文复现:数字化转型能促进企业的高质量发展吗》,表面上聚焦于经济学领域中数字化转型对企业高质量发展影响的研究,实则是一份涵盖多学科交叉的科研仿真代码资源合集。资源以Matlab、Simulink、Python为主要工具,系统整合了电力系统仿真、微电网优化调度、路径规划、信号处理、图像处理、机器学习预测模型等方向的可复现算法与仿真模型。尽管标题指向经济学实证分析,但内容重心在于提供顶级期刊论文的复现代码,如企业全要素生产率(TFP)测算方法(OL、FE、LP、OP、GMM)、风光储氢系统优化、需求响应与综合能源系统调度等,并融合智能优化算法与深度学习技术进行数据建模与预测分析,体现出极强的工程化与科研实用性。; 适合人群:具备一定编程基础,熟练掌握Matlab/Simulink/Python等仿真工具,从事工程仿真、经济实证研究或交叉学科科研工作的研究生、高校教师及科研人员。; 使用场景及目标:① 复现经济学顶刊论文中的计量经济模型,深入探究数字化转型对企业全要素生产率的影响机制;② 借助提供的代码资源开展电力系统故障仿真、微电网优化、多能系统调度等科研项目的算法验证与仿真分析;③ 应用机器学习与深度学习模型完成负荷预测、风电光伏出力预测、电池健康状态评估等典型实证任务; 阅读建议:此资源虽冠以经济学论文之名,实质为多领域高价值仿真代码集成,建议读者依据自身研究方向筛选适配内容,优先关注“顶刊复现”“论文复现”类项目,结合配套数据与代码进行实证推演,并通过公众号“荔枝科研社”获取完整资料与持续技术支持。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值