2016多校第四题 GCD(hdu 5726)(*)
2016-07-21 09:57
239 查看
题意:给一个数组a,大小为n,接下来有m个询问,每次询问给出l、r,定义f[l,r]=gcd(al,al+1,...,ar),问f[l,r]的值 和 有多少对(l',r')使得f[l',r']=f[l,r]。n<=10万,m<=10万,1<=l<=r<=n,1<=l'<=r'。
题解:官方题解很6啊,对于每个点维护以这个点为右端点的所有区间,将这些区间按照GCD一段段合并在一块,然后递推维护下一个点的信息,将结果保存在map中,时间复杂度就是O(nlog(n))。
#include <set>
#include <map>
#include <queue>
#include <deque>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cassert>
#include <iostream>
#include <algorithm>
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF 0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 100010
#define MOD 1000000007
#define zero(x) (((x)>0?(x):-(x))<eps)
vector <pair <int, int> > gcd[maxn];
map <int, long long> mp;
int n;
int num[maxn];
int main()
{
int t,C = 1;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(int i = 1; i <= n; i++)
scanf("%d",&num[i]);
mp.clear();
for(int k = 1; k <= n; k++)
{
int last = -1;
for(int i = 0; i < gcd[k-1].size(); i++)
{
int u = gcd[k-1][i].first;
int v = gcd[k-1][i].second;
int g = __gcd(u,num[k]);
if(g == last)
continue;
last = g;
gcd[k].push_back(make_pair(g,v));
}
if(num[k] != last)
gcd[k].push_back(make_pair(num[k],k));
for(int i = 0; i < gcd[k].size(); i++)
mp[gcd[k][i].first] += (i + 1 == gcd[k].size()?k+1:gcd[k][i+1].second) - gcd[k][i].second;
}
int q;
scanf("%d",&q);
printf("Case #%d:\n",C++);
while(q--)
{
int s,t,k;
scanf("%d%d",&s,&t);
for(k = 0; k < gcd[t].size(); k++)
if(gcd[t][k].second > s)
break;
printf("%d %I64d\n",gcd[t][k-1].first,mp[gcd[t][k-1].first]);
}
for (int i = 0; i <= n; ++i)
gcd[i].clear();
}
return 0;
}
题解:官方题解很6啊,对于每个点维护以这个点为右端点的所有区间,将这些区间按照GCD一段段合并在一块,然后递推维护下一个点的信息,将结果保存在map中,时间复杂度就是O(nlog(n))。
#include <set>
#include <map>
#include <queue>
#include <deque>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cassert>
#include <iostream>
#include <algorithm>
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF 0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 100010
#define MOD 1000000007
#define zero(x) (((x)>0?(x):-(x))<eps)
vector <pair <int, int> > gcd[maxn];
map <int, long long> mp;
int n;
int num[maxn];
int main()
{
int t,C = 1;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(int i = 1; i <= n; i++)
scanf("%d",&num[i]);
mp.clear();
for(int k = 1; k <= n; k++)
{
int last = -1;
for(int i = 0; i < gcd[k-1].size(); i++)
{
int u = gcd[k-1][i].first;
int v = gcd[k-1][i].second;
int g = __gcd(u,num[k]);
if(g == last)
continue;
last = g;
gcd[k].push_back(make_pair(g,v));
}
if(num[k] != last)
gcd[k].push_back(make_pair(num[k],k));
for(int i = 0; i < gcd[k].size(); i++)
mp[gcd[k][i].first] += (i + 1 == gcd[k].size()?k+1:gcd[k][i+1].second) - gcd[k][i].second;
}
int q;
scanf("%d",&q);
printf("Case #%d:\n",C++);
while(q--)
{
int s,t,k;
scanf("%d%d",&s,&t);
for(k = 0; k < gcd[t].size(); k++)
if(gcd[t][k].second > s)
break;
printf("%d %I64d\n",gcd[t][k-1].first,mp[gcd[t][k-1].first]);
}
for (int i = 0; i <= n; ++i)
gcd[i].clear();
}
return 0;
}
相关文章推荐
- Web API应用架构设计分析(1)
- 学习《spring 3.x企业应用开发实战》之基于@AspectJ和Schema的AOP
- ios UIImagePickerController 自定义裁剪框 限制比例等
- 循环链表建立、删除
- iOS开发:UITableView的deleteRowsAtIndexPaths方法的使用
- HDOJ 1997 汉诺塔VII(汉诺塔系列,判断排放是否正确)
- Ext.grid.Panel 正确显示日期
- openCV图像形态学
- 【打CF,学算法——二星级】Codeforces 699B - One Bomb (技巧)
- Netbeans自定义折叠代码
- Android 开发之判断当前是否有网络可用
- div居中
- Android Studio快捷键
- 决心书
- 在linux内核编程 不能include <stdio.h>的原因
- 二叉树的镜像
- 在codeBlock中开发vlc的实例
- 【杭电1872】稳定排序
- 转载整理一些网址
- 如何有效应对游戏行业DDoS攻击?