2016 ACM ICPC Asia Region - Tehran
2016 ACM ICPC Asia Region - Tehran
A - Tax
题目描述:算税。
solution
模拟。
B - Key Maker
题目描述:给出\(n\)个序列,给定一个序列,问\(n\)个序列中有多少个序列满足对应位的值小于或等于给定序列的值。
solution
模拟。
C - IOI 2017 Logo
题目描述:有\(m\)件作品,\(n\)个人投票,每个人可以选三件作品,分别给\(1, 2, 3\)分,每件作品按总分排序,总分相同按得\(3\)分的数量排序,还是相同,则按得\(2\)分得数量排序。输出第一的编号(可能有多个)
solution
模拟。
D - MicroRNA Ranking
题目描述:给出\(m\)个\(n\)的全排列,求一个\(n\)的全排列\(a_i\),满足对于\(i<j\),至少在一半的全排列中,\(a_i\)排在\(a_j\)的前面,求字典序最小的\(a_i\),或无解。
solution
统计出\(cnt[i][j]\)表示\(i\)排在\(j\)前面的次数,如果\(cnt[i][j]\)小于一半,则\(j\)到\(i\)连一条边,表示\(j\)要排在\(i\)的前面。构出的图如果有环,或者不连通则无解。然后求字典序最小的拓扑序(用优先队列维护即可)。
时间复杂度:\(O(n^2(m+logn))\)
E - Möbius Strip
题目描述:普通的莫比乌斯带称为\(1\)型,拧了两次的为\(2\)型,即拧了\(n\)次为\(n\)型。现在有一个剪的操作:沿着莫比乌斯带的宽度\(1/3\)处划线,然后沿线剪下来,得到两个莫比乌斯带。现在给出两组莫比乌斯带,经过若干次剪的操作,使得两组相等。
solution
当\(n\)为偶数时,剪出来的两个莫比乌斯带都是\(n\)型;当\(n\)为奇数时,剪出来的一个是\(n\)型,一个是\(2n+2\)型。所以将奇数的都剪一遍,然后判断两组莫比乌斯带是不是相等(偶数数字相等即可,数量不必相等,奇数要完全相等)。
时间复杂度:\(O(n)\)
F - Expression
题目描述:带分式的表达式计算。
solution
大模拟。
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef vector<LL> vint;
const int cut=int(1e8);
const int maxn=510;
class gjd
{
vint a;
int len;
public:
gjd(int x=0)
{
a.clear();
a.push_back(0);
if (x==0)
{
a.push_back(0);
len=1;
return;
}
len=0;
while (x)
{
++len;
a.push_back(x%cut);
x/=cut;
}
}
gjd(string st)
{
a.clear();
a.push_back(0);
len=(st.length()-1)/8+1;
for (int i=1; i<=len; ++i)
{
int en=int(st.length())-(i-1)*8;
int start=max(0, int(st.length())-i*8);
int num=0;
for (int j=start; j<en; ++j)
num=num*10+st[j]-'0';
a.push_back(num);
}
}
inline bool checkzero()
{
return (len==1 && a[1]==0);
}
bool operator <= (const gjd &b) const
{
if (len!=b.len) return len<b.len;
for (int i=len; i; --i)
if (a[i]!=b.a[i]) return a[i]<b.a[i];
return true;
}
gjd operator + (const gjd &b) const
{
gjd ans;
ans.len=max(len, b.len);
ans.a.assign(ans.len+6, 0);
for (int i=1; i<=len; ++i) ans.a[i]=a[i];
for (int i=1; i<=b.len; ++i) ans.a[i]+=b.a[i];
for (int i=1; i<=ans.len; ++i)
{
ans.a[i+1]+=ans.a[i]/cut;
ans.a[i]%=cut;
}
if (ans.a[ans.len+1]) ++ans.len;
ans.a.resize(ans.len+1);
//ans.a.shrink_to_fit();
return ans;
}
gjd operator - (const gjd &b) const
{
gjd ans=(*this);
for (int i=1; i<=b.len; ++i) ans.a[i]-=b.a[i];
for (int i=1; i<ans.len; ++i)
{
if (ans.a[i]<0)
{
ans.a[i]+=cut;
ans.a[i+1]--;
}
}
while (ans.len>1 && ans.a[ans.len]==0) --ans.len;
ans.a.resize(ans.len+1);
//ans.a.shrink_to_fit();
return ans;
}
gjd operator * (const gjd &b) const
{
gjd ans;
ans.len=len+b.len-1;
ans.a.assign(ans.len+6, 0);
for (int i=1; i<=len; ++i)
for (int j=1; j<=b.len; ++j)
ans.a[i+j-1]+=a[i]*b.a[j];
for (int i=1; i<=ans.len; ++i)
{
ans.a[i+1]+=ans.a[i]/cut;
ans.a[i]%=cut;
}
while (ans.a[ans.len+1])
{
++ans.len;
ans.a[ans.len+1]=ans.a[ans.len]/cut;
ans.a[ans.len]%=cut;
}
while (ans.len>1 && ans.a[ans.len]==0) --ans.len;
ans.a.resize(ans.len+1);
//ans.a.shrink_to_fit();
return ans;
}
gjd operator / (const gjd &b) const
{
if (len<b.len) return gjd(0);
gjd ans;
ans.len=len-b.len+1;
ans.a.assign(ans.len+6, 0);
for (int i=ans.len; i; --i)
{
int L=0, R=cut;
while (L+1<R)
{
int mid=(L+R)>>1;
ans.a[i]=mid;
if (b*ans<=(*this)) L=mid;
else R=mid;
}
ans.a[i]=L;
}
while (ans.len>1 && ans.a[ans.len]==0) --ans.len;
ans.a.resize(ans.len+1);
//ans.a.shrink_to_fit();
return ans;
}
gjd operator % (const gjd &b) const
{
return (*this)-(*this)/b*b;
}
void print()
{
printf("%lld", a[len]);
for (int i=len-1; i; --i)
{
int x=a[i];
for (int j=cut/10; j; j/=10)
{
printf("%d", x/j);
x%=j;
}
}
}
};
gjd gcd(gjd , gjd );
class frac
{
gjd numerator, denominator;
public:
frac(int a=0, int b=1)
{
numerator=gjd(a);
denominator=gjd(b);
}
frac(gjd a, gjd b)
{
numerator=a;
denominator=b;
}
frac operator + (const frac &b) const
{
gjd tnum=numerator*b.denominator+b.numerator*denominator;
gjd tde=denominator*b.denominator;
//gjd d=gcd(tnum, tde);
//return frac(tnum/d, tde/d);
return frac(tnum, tde);
}
frac operator * (const frac &b) const
{
gjd tnum=numerator*b.numerator;
gjd tde=denominator*b.denominator;
//gjd d=gcd(tnum, tde);
//return frac(tnum/d, tde/d);
return frac(tnum, tde);
}
frac operator / (frac b) const
{
swap(b.numerator, b.denominator);
return (*this)*b;
}
void print()
{
gjd d=gcd(numerator, denominator);
(numerator/d).print();
putchar('/');
(denominator/d).print();
putchar('\n');
}
};
int n, m;
char st[maxn][maxn];
int len[maxn];
gjd gcd(gjd b, gjd c)
{
if (c.checkzero()) return b;
else return gcd(c, b%c);
}
void getst(int nid)
{
char ch;
while ((ch=getchar())!='-' && ch!='+' && ch!='*' && ch!=' ' && (ch<'0' || ch>'9'));
len[nid]=1;
st[nid][1]=ch;
while ((ch=getchar())=='-' || ch=='+' || ch=='*' || ch==' ' || (ch>='0' && ch<='9')) st[nid][++len[nid]]=ch;
}
void read()
{
m=0;
for (int i=1; i<=n; ++i)
{
getst(i);
m=max(m, len[i]);
}
for (int i=1; i<=n; ++i)
for (int j=len[i]+1; j<=m; ++j)
st[i][j]=' ';
}
frac getnum(int x, int &y, int lim)
{
string str="";
while (y<=lim && isdigit(st[x][y])) str+=st[x][y++];
return frac(gjd(str), gjd(1));
}
frac dfs(int up, int down, int le, int ri)
{
stack<frac> num;
stack<int> oper;
while (!num.empty()) num.pop();
while (!oper.empty()) oper.pop();
int nid=0;
for (int j=le; j<=ri; ++j)
{
for (int i=up; i<=down; ++i)
if (st[i][j]!=' ')
{
nid=i;
break;
}
if (nid!=0) break;
}
for (int j=le; j<=ri; )
{
if (st[nid][j]==' ') { ++j; continue; }
if (st[nid][j]>='0' && st[nid][j]<='9') num.push(getnum(nid, j, ri));
else
if (st[nid][j]=='-')
{
int k=j;
while (k<=ri && st[nid][k]=='-') ++k;
num.push(dfs(up, nid-1, j, k-1)/dfs(nid+1, down, j, k-1));
j=k;
}
else
{
int level=st[nid][j]=='*';
while (!oper.empty() && oper.top()>=level)
{
frac c=num.top();
num.pop();
frac b=num.top();
num.pop();
if (oper.top()) num.push(b*c);
else num.push(b+c);
oper.pop();
}
oper.push(level);
++j;
}
}
while (!oper.empty())
{
frac c=num.top();
num.pop();
frac b=num.top();
num.pop();
if (oper.top()) num.push(b*c);
else num.push(b+c);
oper.pop();
}
return num.top();
}
int main()
{
while (scanf("%d", &n)==1 && n)
{
read();
dfs(1, n, 1, m).print();
}
return 0;
}
G - Elections
题目描述:有\(n\)座城市,每个城市的选票为\(a_i\),能赢得每个城市的选票的概率为\(p_i\),若得到的总票数大于总票数的一半,则胜利,问胜利的概率。
solution
背包\(dp\)
时间复杂度:\(O(nm)\)(\(m\)为总票数)
H - Explosion at Cafebazaar
题目描述:有\(n\)个点,\(m\)条边的有向图,每个点初始时有\(1bit\)的信息,然后每个时刻复制自己的信息给每一个相邻的点,然后删掉自己的信息,最后接受信息。问最后有多少个点会信息爆炸。
solution
强连通缩点。如果一个点内的环有度大于\(2\)的点,则这个点信息爆炸;否则这个点能到的点(这个点要是一个环)信息爆炸。信息爆炸的点能到的点也是信息爆炸。
时间复杂度:\(O(m+n)\)
I - Linear Galaxy
题目描述:在\(x\)轴上有\(2^n+1\)个点,从中选出\(2^{n-1}+1\)个点,加边(边权为两点间的距离)构成一个简单环,求出最短边的最大值。
J - Rahyab
题目描述:给出一个有向图,一个起点,一个终点,现在有\(C\)个人要从起点走到终点,每个人确定自己走的线路,一个人的花费等于他走的边中最多人经过的边的人数的平方。求所有人总花费的最小值。
solution
显然,最终的线路只会相同,分离,不会边相交。所有可以用网络流求出从起点到终点的边不相交的路径数,将人平均分配到这几条路径上。
时间复杂度:\(O(n^2m)\)(远小于)
K - Base Stations
题目描述:给出平面上的\(n\)个点,每个点有一个\(k\)值,求一个距离最远的\(k\)值不同的点对的距离的平方。
solution
\(KD-TREE\)
时间复杂度:\(O(n\sqrt{n})\)
L - Skeletons
题目描述:不太明白。
2016 ACM ICPC Asia Region - Tehran的更多相关文章
- 2016 ACM/ICPC Asia Regional Shenyang Online 1003/HDU 5894 数学/组合数/逆元
hannnnah_j’s Biological Test Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K ...
- 2016 ACM/ICPC Asia Regional Qingdao Online 1001/HDU5878 打表二分
I Count Two Three Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others ...
- 2016 ACM/ICPC Asia Regional Shenyang Online 1009/HDU 5900 区间dp
QSC and Master Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) ...
- 2016 ACM/ICPC Asia Regional Shenyang Online 1007/HDU 5898 数位dp
odd-even number Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)T ...
- 2016 ACM/ICPC Asia Regional Dalian Online 1002/HDU 5869
Different GCD Subarray Query Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K ( ...
- 2016 ACM/ICPC Asia Regional Dalian Online 1006 /HDU 5873
Football Games Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)To ...
- HDU 5874 Friends and Enemies 【构造】 (2016 ACM/ICPC Asia Regional Dalian Online)
Friends and Enemies Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Othe ...
- HDU 5889 Barricade 【BFS+最小割 网络流】(2016 ACM/ICPC Asia Regional Qingdao Online)
Barricade Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total S ...
- HDU 5875 Function 【倍增】 (2016 ACM/ICPC Asia Regional Dalian Online)
Function Time Limit: 7000/3500 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total ...
随机推荐
- (转)Nginx图片服务器
本文转至博客http://wenxin2009.iteye.com/blog/2117079 Nginx搭建图片服务器 Nginx下载地址:http://nginx.org/en/download.h ...
- BZOJ4939 Ynoi2016掉进兔子洞(莫队+bitset)
容易发现要求三个区间各数出现次数的最小值.考虑bitset,不去重离散化后and一发就可以了.于是莫队求出每个区间的bitset.注意空间开不下,做多次即可.输出的东西错了都能调一年服了我了. #in ...
- 51nod1238 最小公倍数之和 V3 莫比乌斯函数 杜教筛
题意:求\(\sum_{i = 1}^{n}\sum_{j = 1}^{n}lcm(i, j)\). 题解:虽然网上很多题解说用mu卡不过去,,,不过试了一下貌似时间还挺充足的,..也许有时间用phi ...
- unity3D AR涂涂乐制作浅谈
unity3D AR涂涂乐制作浅谈 AR为现在是虚拟现实较为火爆的一个技术,其中有个比较炫酷的就是AR涂涂乐的玩法,这个技术可以把扫描到的图片上的纹理 粘贴到模型上实现为模型上色的功能,但是我们需要怎 ...
- NOIP2003 传染病控制 【搜索 + 卡时】
题目背景 近来,一种新的传染病肆虐全球.蓬莱国也发现了零星感染者,为防止该病在蓬莱国大范围流行,该国政府决定不惜一切代价控制传染病的蔓延.不幸的是,由于人们尚未完全认识这种传染病,难以准确判别病毒携带 ...
- 洛谷P3933 Chtholly Nota Seniorious 【二分 + 贪心 + 矩阵旋转】
威廉需要调整圣剑的状态,因此他将瑟尼欧尼斯拆分护符,组成了一个nnn行mmm列的矩阵. 每一个护符都有自己的魔力值.现在为了测试圣剑,你需要将这些护符分成 A,B两部分. 要求如下: 圣剑的所有护符, ...
- python基础----内置函数----匿名函数(lambda)
Python3版本所有的内置函数: 1. abs() 获取绝对值 >>> abs(-) >>> abs() >>> abs() >>& ...
- Oracle恢复误删除表操作语句
一.表和数据恢复 1.从回收站里查询被删除的表 select object_name,original_name,partition_name,type,ts_name,createtime,drop ...
- CCPC2018-A-Buy and Resell
Problem Description The Power Cube is used as a stash of Exotic Power. There are n cities numbered 1 ...
- UVA 1575 Factors
https://vjudge.net/problem/UVA-1575 题意: 令f(k)=n 表示 有n种方式,可以把正整数k表示成几个数的乘积的形式. 例 10=2*5=5*2,所以f(10)=2 ...