BC div2补题以及 复习模除 逆元__BestCoder Round #78 (div.2)
第一题没话说 智商欠费 加老柴辅导终于过了
需要在意的是数据范围为2的63次方-1 三个数相加肯定爆了
四边形的定义 任意边小于其余三边之和
换句话说就是 最长边小于其余三边之和
这样的话问题转化为 最长边依次减其余三边的结果是否小于等于0
还有一点是题目出现0边 即最小边不为0 想得太多反而把0也算为合法。。。。
问题只需要 sort一下 判断a[0]==0||a[3]-a[2]-a[1]-a[0]>=0 输出NO //存在0边且最大边大于其他边之和
第二题 好多种姿势 题目链接http://bestcoder.hdu.edu.cn/contests/contest_chineseproblem.php?cid=683&pid=1002
官方题解
我们令dp[i][j]表示在前i个数中,选出若干个数使得它们的gcd为j的方案数,于是只需要枚举第i+1个数是否被选中来转移就可以了
令第i+1个数为v,当考虑dp[i][j]的时候,我们令$dp[i+1][j] += dp[i]j,dp[i+1][gcd(j,v)] += dp[i]j
复杂度O(N*MaxV) MaxV 为出现过的数的最大值
其实有O(MaxV *log(MaxV))的做法,我们考虑记f[i]表示从这些数中选择若干个数,使得他们的gcd是i的倍数的方案数。假如有K个数是i的倍数,则 f[i]=2^K-1,再用g[i]表示从这些数中选择若干个数,使得他们的gcd是i的方案数,则g[i]=f[i] - g[j] (对于所有j是i的倍数)。
由调和级数可以得到复杂度为O(MaxV *log(MaxV))
DP之二维数组转移
我们把dp[i][j]作为考虑了第i个数GCD为j的方案数
直接gcd会超时 所以我们打个表GCD
那么dp[i][j]+=dp[i-1][j]; dp[i][GCD[j][v[i]]]+=dp[i-1][j]; 然后就可以转移辣;
#include<cstdio>
#include<map>
//#include<bits/stdc++.h>
#include<vector>
#include<stack>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
#include<queue>
#include<cstdlib>
#include<climits>
#define PI acos(-1.0)
#define INF 0x3fffffff
using namespace std;
typedef long long ll;
typedef __int64 int64;
const ll mood=1e9+;
const int64 Mod=;
const double eps=1e-;
const int N=;
const int MAXN=;
typedef int rl;
inline void r(rl&num){
num=;rl f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<='')num=num*+ch-'',ch=getchar();
num*=f;
}
int gcd(int a,int b)
{
return b==?a:gcd(b,a%b);
}
int v[N];
int GCD[N][N];
int64 dp[N][N];
int main()
{
for(int i=;i<;i++)
{
for(int j=;j<=i;j++)
{
GCD[i][j]=GCD[j][i]=gcd(i,j);
}
}
int ci;
r(ci);
while(ci--)
{
int n;
r(n);
int mx=-;
for(int i=;i<=n;i++)
{
r(v[i]);
mx=max(mx,v[i]);
dp[i][v[i]]=;
}
for(int i=;i<=n;i++)
{
for(int j=;j<=mx;j++)
{
dp[i][j]+=dp[i-][j];
dp[i][j]%=Mod;
dp[i][GCD[j][v[i]]]+=dp[i-][j];
dp[i][GCD[j][v[i]]]%=Mod;
}
}
int64 ans=; for(int i=;i<=mx;i++)
{ ans+=(dp[n][i]*i)%Mod;
ans%=Mod;
}
memset(dp,,sizeof(dp));
memset(v,,sizeof(v));
printf("%I64d\n",ans);
}
return ;
}
二维
仔细想了一下 觉得可以优化为滚动数组 试了好久不对 最后瞎蒙
每个数都多考虑了一次 所以/2需要乘逆元 正好1e8+7是素数
Mod为素数,那么还可以根据费马小定理得到逆元为 2的(Mod-2)次方%Mod
即除2等于乘2的(Mod-2)次方%Mod
所以加了一个快速幂 但是优化为滚动数组后 时间增加了一丢丢 但空间大幅度减少
16757862 | 2016-04-03 12:45:34 | Accepted | 5656 | 2511MS | 5504K | 1925 B | G++ | zxMrlc |
16755798 | 2016-04-03 00:36:10 | Accepted | 5656 | 2449MS | 13404K | 1722 B | G++ | zxMrlc |
但是姿势老感觉有问题 等wtw学长指点后我再改改 还有官方的第二个姿势还没有学会。。。衰
#include<cstdio>
#include<map>
//#include<bits/stdc++.h>
#include<vector>
#include<stack>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
#include<queue>
#include<cstdlib>
#include<climits>
#define PI acos(-1.0)
#define INF 0x3fffffff
using namespace std;
typedef long long ll;
typedef __int64 int64;
const ll mood=1e9+;
const int64 Mod=;
const double eps=1e-;
const int N=;
const int MAXN=;
typedef int rl;
inline void r(rl&num){
num=;rl f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<='')num=num*+ch-'',ch=getchar();
num*=f;
}
int gcd(int a,int b)
{
return b==?a:gcd(b,a%b);
}
int v[N];
int GCD[N][N];
int64 dp[N];
int main()
{
int64 xx=Mod-;
int64 an=,t=;
while(xx>)
{
if(xx&) an*=t;
xx/=;
an%=Mod;
t*=t;
t%=Mod;
}
for(int i=;i<;i++)
{
for(int j=;j<=i;j++)
{
GCD[i][j]=GCD[j][i]=gcd(i,j);
}
} int ci;
r(ci);
while(ci--)
{
int n;
r(n);
int mx=-;
for(int i=;i<=n;i++)
{
r(v[i]);
mx=max(mx,v[i]);
} for(int i=;i<=n;i++)
{
dp[v[i]]++;
for(int j=;j<=mx;j++)
{
// dp[i][j]+=dp[i-1][j];
dp[j]%=Mod;
dp[GCD[j][v[i]]]+=dp[j];
dp[GCD[j][v[i]]]%=Mod;
}
}
int64 ans=;
//for(int i=1;i<=mx;i++) cout<<dp[i]<<endl;
for(int i=;i<=mx;i++)
{ ans+=(dp[i]*i)%Mod;
ans%=Mod;
}
memset(dp,,sizeof(dp));
memset(v,,sizeof(v));
printf("%I64d\n",ans*an%Mod);
}
return ;
}
滚动数组
我们每次加入的数据会导致翻倍 所以刚才改为加完/2;
因为添加的v[i]导致的影响就是 当前位置dp[v[i]]多1 即方案数多了选自己的 所以在循环结尾-1就ok了 。。。根本不需要模除 但时间特么变大了
还是有点模糊的 不太清楚到底怎么回事。
16758163 | 2016-04-03 13:17:49 | Accepted | 5656 | 2636MS | 5504K | 1730 B | G++ | zxMrlc |
#include<cstdio>
#include<map>
//#include<bits/stdc++.h>
#include<vector>
#include<stack>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
#include<queue>
#include<cstdlib>
#include<climits>
#define PI acos(-1.0)
#define INF 0x3fffffff
using namespace std;
typedef long long ll;
typedef __int64 int64;
const ll mood=1e9+;
const int64 Mod=;
const double eps=1e-;
const int N=;
const int MAXN=;
typedef int rl;
inline void r(rl&num){
num=;rl f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<='')num=num*+ch-'',ch=getchar();
num*=f;
}
int gcd(int a,int b)
{
return b==?a:gcd(b,a%b);
}
int v[N];
int GCD[N][N];
int64 dp[N];
int main()
{
for(int i=;i<;i++)
{
for(int j=;j<=i;j++)
{
GCD[i][j]=GCD[j][i]=gcd(i,j);
}
} int ci;
r(ci);
while(ci--)
{
int n;
r(n);
int mx=-;
for(int i=;i<=n;i++)
{
r(v[i]);
mx=max(mx,v[i]);
} for(int i=;i<=n;i++)
{
dp[v[i]]++;
for(int j=;j<=mx;j++)
{
// dp[i][j]+=dp[i-1][j];
dp[j]%=Mod;
dp[GCD[j][v[i]]]+=dp[j];
dp[GCD[j][v[i]]]%=Mod;
}
dp[v[i]]--;
}
int64 ans=;
for(int i=;i<=mx;i++)
{ ans+=(dp[i]*i)%Mod;
ans%=Mod;
}
memset(dp,,sizeof(dp));
memset(v,,sizeof(v));
printf("%I64d\n",ans%Mod);
}
return ;
}
滚动第二次优化
BC div2补题以及 复习模除 逆元__BestCoder Round #78 (div.2)的更多相关文章
- codeforces round 422 div2 补题 CF 822 A-F
A I'm bored with life 水题 #include<bits/stdc++.h> using namespace std; typedef long long int LL ...
- codeforces round 421 div2 补题 CF 820 A-E
A Mister B and Book Reading O(n)暴力即可 #include<bits/stdc++.h> using namespace std; typedef lon ...
- Codeforces round 419 div2 补题 CF 816 A-E
A Karen and Morning 水题 注意进位即可 #include<bits/stdc++.h> using namespace std; typedef long long i ...
- codeforces 447 A-E div2 补题
A DZY Loves Hash 水题 #include<iostream> #include<cstdio> #include<cstdlib> #include ...
- codeforces round 418 div2 补题 CF 814 A-E
A An abandoned sentiment from past 水题 #include<bits/stdc++.h> using namespace std; int a[300], ...
- codeforces round 417 div2 补题 CF 812 A-E
A Sagheer and Crossroads 水题略过(然而被Hack了 以后要更加谨慎) #include<bits/stdc++.h> using namespace std; i ...
- codeforces round 416 div2 补题 CF 811 A B C D E
A. Vladik and Courtesy 水题略过 #include<cstdio> #include<cstdlib> #include<cmath> usi ...
- codeforces round 420 div2 补题 CF 821 A-E
A Okabe and Future Gadget Laboratory 暴力 #include<bits/stdc++.h> using namespace std; typedef l ...
- 2018 HDU多校第三场赛后补题
2018 HDU多校第三场赛后补题 从易到难来写吧,其中题意有些直接摘了Claris的,数据范围是就不标了. 如果需要可以去hdu题库里找.题号是6319 - 6331. L. Visual Cube ...
随机推荐
- monkey基本命令及脚本编写
Monkey 是Android自带的黑盒测试工具,一般通过随机触发界面事件,来确定应用是否会发生异常,多用于android应用的稳定性.压力测试 基本命令: adb shell monkey [op ...
- 使用git命令将本地项目上传到GIthub
1.(先进入项目文件夹)通过命令 git init 把这个目录变成git可以管理的仓库 2.把文件添加到版本库中,使用命令 git add .添加到暂存区里面去,不要忘记后面的小数点“.”,意为添加文 ...
- python 并发编程之IO 模型
首先说一下 IO 发生时涉及的对象和步骤.以read 为例,会经历两个阶段: 1)等待数据准备 2)将数据从内核拷贝到进程中 二,阻塞Io(blocking IO) 在 Linux中 默认情况下所有 ...
- CCF 201604-1 折点计数 (水题,暴力)
问题描述 给定n个整数表示一个商店连续n天的销售量.如果某天之前销售量在增长,而后一天销售量减少,则称这一天为折点,反过来如果之前销售量减少而后一天销售量增长,也称这一天为折点.其他的天都不是折点.如 ...
- Apache CXF简介
Apache CXF是一个开源的,全功能的,容易使用的Web服务框架.CXF是由Celtix和XFire合并,在Apache软件基金会共同完成的.CXF的名字来源于"Celtix" ...
- Codeforces 176B【计数DP】
题意: 给你两个串s1,s2和一个K, 有一种操作是在一个串切开然后交换位置, 问s1有多少种方法经过K次这样的操作变成s2: 思路: (从来没接触过计数DP...还是太菜...参考了[大牛blog] ...
- 萌新学习图的强连通(Tarjan算法)笔记
--主要摘自北京大学暑期课<ACM/ICPC竞赛训练> 在有向图G中,如果任意两个不同顶点相互可达,则称该有向图是强连通的: 有向图G的极大强连通子图称为G的强连通分支: Tarjan算法 ...
- 51nod1181【素数筛】
思路: 直接就是筛出素数,然后我很撒比的从那个地方往后for找一个位置也是质数的输出: #include <bits/stdc++.h> using namespace std; type ...
- 3D游戏中人物换装解决方案
换装基本上是每个网游都必须有的一个功能,每种网游的做法都各有不同,有些是换掉整个模型,有些则是通过可以换掉模型的一个部分完成.前者属于整体换,相对简单些:后者则是通过部分替换实现,目前用的比较多,本文 ...
- c# UserControl 如何在父窗口中显示Scroll
1. UserControl 不应该是Dock = Fill, 而应该使用 Anchor 2. 父窗口(或者父控件,或者容器)的AutoScroll = True 3.UserControl 加入容器 ...