暑假集训日记——7.7（月赛+scc习题）

最新推荐文章于 2020-10-15 02:08:56 发布
最新推荐文章于 2020-10-15 02:08:56 发布
阅读量408 收藏
点赞数
分类专栏： 线性筛 暑假集训 scc 思维题型/数论相关
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43870114/article/details/94961545
版权

A - 迷宫城堡
题解：
两两可达，就是判断一下这个图是否是一个强联通图，模板题。

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;
#define mp make_pair

const int N=1e5+10;
const long long INF=1e18;
const double eps=0.0000001;
const ll mod=1e9+7;
int n,m;
vector<int>edge[10005];
int dfn[N], low[N];///dfn[]表示深搜的步数，low[u]表示u或u的子树能够追溯到的最早的栈中节点的次序号
int sccno[N];///缩点数组，表示某个点对应的缩点值  有的写成belong
int step; ///dfs的步数
int scc_cnt;///强连通分量个数
//vector<int> scc[10005];///得出来的缩点，scc[i]里面存i这个缩点具体缩了哪些点
stack<int> S; ///判断后向边的栈

void dfs(int o)
{
    dfn[o]=low[o]=++step;
    S.push(o);
    for(int i=0;i<edge[o].size();i++)
    {
        int v=edge[o][i];
        if(!dfn[v])
        {
            dfs(v);
            low[o] = min(low[o], low[v]);
        }
        if(!sccno[v])
        {
            low[o] = min(low[o], dfn[v]);
        }
    }
    if(dfn[o]==low[o])
    {
        scc_cnt++;
        //scc[scc_cnt].clear();
        while(1)
        {
            int x = S.top();
            S.pop();
            //if (sccno[x] != scc_cnt) scc[scc_cnt].push_back(x);
            sccno[x] = scc_cnt;
            if (x == o) break;
        }
    }
}

void tarjan(int n)
{
    memset(sccno, 0, sizeof(sccno));
    memset(dfn, 0, sizeof(dfn));
    step = scc_cnt = 0;
    for (int i = 1; i <=n; i++)
        if (!dfn[i]) dfs(i);
}


int main()
{
   ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
   int x,y;
   while(cin>>n>>m)
   {
       if(n==0&&m==0)
       break;
       for(int i=1;i<=n;i++)
        edge[i].clear();
       for(int i=0;i<m;i++)
       {
           cin>>x>>y;
           edge[x].push_back(y);
       }
       tarjan(n);
       //cout<<scc_cnt<<endl;
       //for(int i=1;i<=n;i++)
       // cout<<low[i]<<" "<<dfn[i]<<endl;
       if(scc_cnt==1)
        cout<<"Yes"<<endl;
       else
        cout<<"No"<<endl;

   }

}

B - Summer Holiday
题解：（缩点）
对于强联通的联通分支来说，只要联系强联通的一个人就可以联系到所有人，所以把强联通缩点，代价取联通中最小的。

由以上操作后，图中就没有强联通关系了，而要使得联系人数目最小，就只联系那些 不得不联系的人（入度为0），因为如果自己不联系的话，其他人也无法联系，而入度不为一的话，一定可以有（入度为0）人能联系到他，因此统计联系入度为零人的代价即为所求。


#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;
#define mp make_pair

const int N=2005;
const long long INF=1e18;
const double eps=0.0000001;
const ll mod=1e9+7;
int n,m;
vector<int>edge[2005];
int pri[N];//联系每个的代价
int x1[N],y1[N],in[N];//in数组存入度
int minpri[N];
int dfn[N], low[N];///dfn[]表示深搜的步数，low[u]表示u或u的子树能够追溯到的最早的栈中节点的次序号
int sccno[N];///缩点数组，表示某个点对应的缩点值  有的写成belong
int step; ///dfs的步数
int scc_cnt;///强连通分量个数
//vector<int> scc[10005];///得出来的缩点，scc[i]里面存i这个缩点具体缩了哪些点
stack<int> S; ///判断后向边的栈

void dfs(int o)
{
    dfn[o]=low[o]=++step;
    S.push(o);
    for(int i=0;i<edge[o].size();i++)
    {
        int v=edge[o][i];
        if(!dfn[v])
        {
            dfs(v);
            low[o] = min(low[o], low[v]);
        }
        if(!sccno[v])
        {
            low[o] = min(low[o], dfn[v]);
        }
    }
    if(dfn[o]==low[o])
    {
        scc_cnt++;
        //scc[scc_cnt].clear();
        while(1)
        {
            int x = S.top();
            S.pop();
            //if (sccno[x] != scc_cnt) scc[scc_cnt].push_back(x);
            minpri[scc_cnt]=min(minpri[scc_cnt],pri[x]);
            sccno[x] = scc_cnt;
            if (x == o) break;
        }
    }
}

void tarjan(int n)
{
    memset(sccno, 0, sizeof(sccno));
    memset(dfn, 0, sizeof(dfn));
    memset(in, 0, sizeof(in));
    step = scc_cnt = 0;
    for (int i = 1; i <=n; i++)
        if (!dfn[i]) dfs(i);
}


int main()
{
   ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
   int x,y;
   while(cin>>n>>m)
   {
       for(int i=1;i<=n;i++)
       {
            cin>>pri[i];
            edge[i].clear();
            minpri[i]=N;
       }
       for(int i=1;i<=m;i++)
       {
           cin>>x1[i]>>y1[i];
           edge[x1[i]].push_back(y1[i]);
       }
       tarjan(n);
       //cout<<scc_cnt<<endl;
       //for(int i=1;i<=n;i++)
       // cout<<low[i]<<" "<<dfn[i]<<endl;
       for(int i=1;i<=m;i++)
       {
           if(sccno[y1[i]]!=sccno[x1[i]])
           in[sccno[y1[i]]]++;
       }
       int sum=0,cnt=0;
        for(int i=1;i<=scc_cnt;i++)
       {
            if(in[i]==0)
            sum+=minpri[i],cnt++;
       }
       cout<<cnt<<" "<<sum<<endl;


   }

}

月赛补题

问题 D: 猫抓老鼠
题解：欧拉函数
欧拉函数是求小于n的互质数…我的天哪！！！！我是会的，居然没反应过来，把欧拉函数和线性筛搞混了，我的天哪！！！当场去世，太菜了…

有n只老鼠，对于二分之一象限来说，横坐标的取值范围为[w/(n+1)+1，w/n]
易证（然后我证了半个点…都是泪

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;
#define mp make_pair

const int N=5e6+10;
const long long INF=1e18;
const double eps=0.0000001;
const ll mod=1e9+7;
ll n,m,k;
int f[N];
int prime[N];
int pre[N];
ll cnt;

void init(){
    f[1]=1;
    prime[0]=prime[1]=1;
    for(int i=2;i<=N-10;i++){
        if(prime[i]==0){
            pre[++cnt]=i;
            f[i]=i-1;
        }
        for(int j=1;j<=cnt&&pre[j]*i<=N-10;j++){
            prime[i*pre[j]]=1;
            f[i*pre[j]]=f[i]*f[pre[j]];
            if(i%pre[j]==0){
                f[i*pre[j]]=f[i]*(f[pre[j]]+1);
                break;
            }
        }
    }
}

int main()
{
   ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
   init();
   cin>>n>>m;
   ll sum=0;
   int l=n/(m+1),r=n/m;
   if(n==m)
    cout<<8<<endl;
   else
    {
        for(int i=l+1;i<=r;i++)
        {
            sum+=f[i]*8;
        }
        cout<<sum<<endl;
    }
}

问题 F: 夏令营
题解：
组合数+错排
第一轮的情况数为：从n中选出n/2个元素作为主场，剩下的作为客场，客场与主场组合，有n/2的全排列种。化简也就是
for(int i=n/2+1;i<=n;i++) ans=ians%mod
然后第二轮的情况就是第一轮的错排
公式：F[N]=(N-1)(F[N-1]+F[N-2])

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;
#define mp make_pair

const int N=5e6+10;
const long long INF=1e18;
const double eps=0.0000001;
const ll mod=1e9+7;
ll n,m,k;
ll f[N];
ll cnt;

void init(int x)
{
    f[1]=0,f[2]=1;
    for(int i=3;i<=x;i++)
    {
        f[i]=(i-1)*(f[i-1]+f[i-2])%mod;
    }
}

int main()
{
   ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
   cin>>n;
   init(n/2);
   ll ans=1;
   for(int i=n/2+1;i<=n;i++)
    ans=i*ans%mod;
   ans=ans*f[n/2]%mod;
   cout<<ans<<endl;


}

问题 C: 一脸严肃的声明: 作业作假的行为是可耻的行为!
题解：
模拟就可以了，一开始没考虑如果数组中没有中位数值的话怎么办，还有记录中位数值的下标要初始化。

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;
#define mp make_pair

const int N=1e6+10;
const long long INF=1e18;
const double eps=0.0000001;
const ll mod=1e9+7;
ll n,m,k;
ll a[N];
ll sum[N];
ll cnt[N];

ll mul(ll a,ll b){return (a%mod*b%mod)%mod;}
ll pre(ll a,ll b){return (a%mod/b%mod)%mod;}

int main()
{
   ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
   cin>>n>>m>>k;
    for(int i=1;i<=n*m;i++)
    {
        cin>>a[i];
        if(a[i]>k)
        sum[(i-1)/m+1]++;
        if(a[i]<k)
        cnt[(i-1)/m+1]++;
    }
    sort(sum+1,sum+n+1);
    sort(cnt+1,cnt+n+1);
    sort(a+1,a+(n*m)+1);
    ll ans,flag=0;
    int pos=(n*m)%2?n*m/2+1:n*m/2;
        if(k==a[pos])
        {
            cout<<0<<" "<<0<<endl;
            return 0;
        }
        else if(k<a[pos])
        {
            ans=0;
            for(int i=1;i<=n*m;i++)
            {
                if(a[i]<=k)
                {
                    ans=i;
                }
                else
                    break;
            }
            flag=1;
        }
        else
        {
            ans=n*m+1;
            for(int i=n*m;i>=1;i--)
            {
                if(a[i]>=k)
                {
                    ans=i;
                }
                else
                    break;
            }
        }
        int ppp=abs(ans-pos),cnts=0;
        if(flag==1)
            for(int i=n;i>=1;i--)
            {
                ppp-=sum[i];
                cnts++;
                if(ppp<=0)
                    break;
            }
        else
        {
             for(int i=n;i>=1;i--)
            {
                ppp-=cnt[i];
                cnts++;
                if(ppp<=0)
                    break;
            }
        }
        cout<<cnts<<" "<<abs(ans-pos)<<endl;

}
THE END GAME
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
SCC 算法设计 c++
06-02
SCC 算法设计 c++ 源码 SCC 算法设计 c++ 源码
Scc Puzzle(水题)
笃行
01-24 467
版权声明：在校学生一枚，文章均为个人心得，有不足之处请不吝赐教！ https://mp.csdn.net/mdeditor#
Scc Puzzle（AtCoder-2333）
Alex_McAvoy的博客
05-30 446
Problem Description Snuke loves puzzles. Today, he is working on a puzzle usingS- andc-shaped pieces. In this puzzle, you can combine twoc-shaped pieces into oneS-shaped piece, as shown in the ...
WUST 1946 繁忙的都市（最小生成树+克鲁斯卡尔算法）
HJY
07-15 357
城市C是一个非常繁忙的大都市，城市中的道路十分的拥挤，于是市长决定对其中的道路进行改造。城市C的道路是这样分布的：城市中有n个交叉路口，有些交叉路口之间有道路相连，两个交叉路口之间最多有一条道路相连接。这些道路是双向的，且把所有的交叉路口直接或间接的连接起来了。每条道路都有一个分值，分值越小表示这个道路越繁忙，越需要进行改造。但是市政府的资金有限，市长希望进行改造的道路越少越好，于是他提出下面的要求： 改造的那些道路能够把所有的交叉路口直接或间接的连通起来。 在满足要求1的情况下，改造的道路尽量少。 在满足
例题1.19 计算器谜题 UVa11549
Skyline
03-15 649
1.题目描述：点击打开链接 2.解题思路：本题是一道普通的模拟题，根据题意易知计算器显示的数将会出现循环，因为显示的n位数一共就有限种情况，而平方的次数是无限多的。所以不妨一个个的模拟，每次都看新得到的数是否以前出现过，如果出现过就跳出循环。但如何判断是否出现过呢？第一种方法是利用STL中的set，同时写一个Next函数求出下一个k值。下面一共给出三个不同的代码，来分别体会一下效率上的巨大差异。
div+scc demo初学
03-19
NULL 博文链接：https://yj10864.iteye.com/blog/753495
TriloTracker：用于PSG + SCC，PSG + FM和PSG（SN7）+ FM的MSX音乐跟踪器
02-13
PSG（AY3）+ Konami SCC的跟踪器。 TTFM（ttrackerFM.asm）： PSG（AY3）+ FM（OPLL）的跟踪器。 TTSMS（ttrackerSMS.asm）： PSG（SN7）+ FM（OPLL）的跟踪器。 此配置适用于Sega Master System声音芯片。 可以...
SCC_SerialVC++_Vc_strengthim4_
10-04
VC++基于串口的监视程序，通过串口发消息的形式，来监视串口的各相操作，在发送方式，可选择循环发送、触发式发送，还可设置两者混合发送。
scc：脱机备忘单和HTML，CSS和JS的快速参考命令行工具
02-03
scc 很棒的备忘单和HTML，CSS和JS快速参考命令行工具。 正在安装 下载项目并运行安装脚本。 ./install.sh 用法 scc [-h] [ -h tml [HTML] | -c ss [CSS] | -js [JS] ] | [-rand {html,css,js}] 例子 scc -js array....
scc：Suil代码编译器-将.scc代码转换为c ++
02-19
scc suil代码编译器项目。 提供用于将suil代码（C ++扩展名）转换为C ++代码的工具
Kosarju模拟：SCC + 例题-班长竞选（缩点）
qq_44654498的博客
04-14 262
Vjudge:kosarju算法讲解+例题
SCC BCC经典题目小结
qq_43555854的博客
10-15 282
hdu 3836 题意： 最少加几条边使得整张图scc。 思路： 先求出强连通分量对scc重新建图，预处理出新图每个点的入度和出度，统计出新图中出度为0和入度为0的点的个数分别为out0和in0,答案就是max(out0,in0)。 poj 2186 题意： 有m对奶牛之间的喜欢关系，喜欢关系有传递性，现在问有多少牛被所有的牛喜欢？（或有哪些牛） 思路： 先求出强连通分量对scc重新建图，预处理出新图每个点的出度，设出度为0的点数为cnt，如果cnt>1,就无解，否则cnt就是1，有解，遍历每个点属于
hdu 1827 Summer Holiday（传递的最小费用）
WEYuLi的专栏
08-09 846
Summer Holiday Time Limit: 10000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1190    Accepted Submission(s): 529 Problem Description To see a World
等比数列三角形
weixin_43870114的博客
05-15 1553
题解： （3|12）3能被12整除 #include<bits/stdc++.h> using namespace std; #define LL long long #define endl '\n' const int mod=1e9+7; LL n; double e=(1.0+sqrt(5.0))/2; LL b[5000]; int main() { whi...
普通递归关系
weixin_43870114的博客
05-29 544
注意：0 0的特判 纠结好久…哎 #include<algorithm> #include<cstdio> #include<cmath> #include<iostream> #include<cstring> #include<set> #include<map> #include<iterator&g...
求解前n个数的约数和
weixin_43870114的博客
04-03 461
#include<cstdio> int n,ans; int main() { scanf("%d",&n); for(int i=1,j;i<=n;i=j+1)//节约时间 { j=n/(n/i); ans+=(n/i)*(j-i+1); } printf("%d",ans); retur...
欧拉函数
weixin_43870114的博客
04-25 406
φ函数的值 　 通式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)……(1-1/pn) 当中p1, p2……pn为x的全部质因数。x是不为0的整数。φ(1)=1（唯一和1互质的数(小于等于1)就是1本身）。 (注意：每种质因数仅仅一个。比方12=223那么φ（12）=12*（1-1/2）*(1-1/3)=4 若n是质数p的k次幂。φ(n)=p^k - p^(k...
暑假集训日记——7.9/7.10（scc习题+codeforce）
weixin_43870114的博客
07-10 365
F. Simple Cycles Edges 题解： 要求求简单环，所以可以用tarjan求出点双连通分量，然后判断这个点双连通分量中的点的个数是否等于边的条数（有且仅有边的条数等于点的个数，才是个简单环），如果是的话把连通分量中的所有边存到一个数组中，最后找出所有边之后对边的编号进行排序输出就好了。 #include<bits/stdc++.h> using namespace st...
SCC2000 驱动
最新发布
05-06
SCC2000 是一种高性能的步进电机驱动器，它是由日本三菱电机公司生产的。SCC2000 驱动器采用闭环控制技术，能够实现高速、高精度、高效率的步进电机控制。其特点主要有以下几个方面： 1. 高速控制：SCC2000 驱动器...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值