2012 Asia Hangzhou/Jinhua Regional Contest
hdu4453 Looploop
好题……虽然有一系列区间翻转区间求和等等平衡树才能玩得起的操作，但翻转的长度和增加的长度都是固定的。
而且题目中给了各种的条件，比如k1<k2≤N..就能有更简单的处理这些操作的办法。
开3个双端队列，第一个存指针后1~k1，第二个存k1+1~k2，第三个存剩下的数字。
那么对于区间翻转和区间求和，只需要修改两个deque的标记。
插入和删除也是只要修改第一个双端队列的头或尾，并且维护后两个队列满足要求
指针左右移也是纸上好好想想就能推出来了，略有些复杂。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <list>
#define pb push_back
#define mp make_pair
#define fi first
#define se second
//#pragma comment(linker, "/STACK:16777216")
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
#define N 1000005
deque<int> Q1 , Q2 , Q3;
int n , m , k1 , k2;
int add , rev;
void insert(int x)
{
int y;
if (rev)
y = Q1.back() , Q1.pop_back() , Q1.push_back(x - add) , Q1.push_back(y);
else y = Q1.front() , Q1.pop_front() , Q1.push_front(x - add) , Q1.push_front(y);
Q3.push_front(Q2.back() + add) , Q2.pop_back();
if (rev)
Q2.push_front(Q1.front()) , Q1.pop_front();
else Q2.push_front(Q1.back()) , Q1.pop_back();
}
void del()
{
int y;
if (rev)
Q1.pop_back() , Q1.push_front(Q2.front()) , Q2.pop_front() , Q2.push_back(Q3.front() - add) , Q3.pop_front();
else
Q1.pop_front() , Q1.push_back(Q2.front()) , Q2.pop_front() , Q2.push_back(Q3.front() - add) , Q3.pop_front();
}
int move(int x)
{
if (x == 1)
{
if (rev)
{
Q1.push_back(Q3.back() - add) , Q3.pop_back();
Q2.push_front(Q1.front()) , Q1.pop_front();
Q3.push_front(Q2.back() + add) , Q2.pop_back();
}
else
{
Q1.push_front(Q3.back() - add) , Q3.pop_back();
Q2.push_front(Q1.back()) , Q1.pop_back();
Q3.push_front(Q2.back() + add) , Q2.pop_back();
}
}
else
{
if (rev)
{
Q3.push_back(Q1.back() + add) , Q1.pop_back();
Q2.push_back(Q3.front() - add) , Q3.pop_front();
Q1.push_front(Q2.front()) , Q2.pop_front();
}
else
{
Q3.push_back(Q1.front() + add) , Q1.pop_front();
Q2.push_back(Q3.front() - add) , Q3.pop_front();
Q1.push_back(Q2.front()) , Q2.pop_front();
}
}
}
void watch()
{
puts("*************");
for (deque<int>::iterator it = Q1.begin() ; it != Q1.end() ; ++ it)
printf("%d " , *it + add); puts("");
for (deque<int>::iterator it = Q2.begin() ; it != Q2.end() ; ++ it)
printf("%d " , *it + add); puts("");
for (deque<int>::iterator it = Q3.begin() ; it != Q3.end() ; ++ it)
printf("%d " , *it); puts("");
puts("*************");
}
int query()
{
if (rev)
return Q1.back() + add;
else return Q1.front() + add;
}
int ca;
void work()
{
int i , j , x , y; char str[10];
while (!Q1.empty()) Q1.pop_back();
while (!Q2.empty()) Q2.pop_back();
while (!Q3.empty()) Q3.pop_back();
scanf("%d%d%d",&m,&k1,&k2);
add = rev = 0;
printf("Case #%d:\n" , ++ ca);
for (i = 1 ; i <= n ; ++ i)
{
scanf("%d",&x);
if (i <= k1) Q1.push_back(x);
else if (i <= k2) Q2.push_back(x);
else Q3.push_back(x);
}
//watch();
while (m --)
{
scanf("%s" , str);
if (*str == 'i')
scanf("%d",&x) , insert(x);
if (*str == 'a')
scanf("%d",&x) , add += x;
if (*str == 'r')
rev ^= 1;
if (*str == 'd')
del();
if (*str == 'm')
scanf("%d",&x) , move(x);
if (*str == 'q')
printf("%d\n" , query());
//watch();
}
}
int main()
{
while (scanf("%d",&n) , n)
work();
return 0;
}
hdu4454 Stealing a Cake
三分
hdu4455 Substrings
好题
hdu4456 Crowd
10000*10000的矩阵，单点增加，询问一个菱形的权值和。
菱形只要转一下就好了(x+y,x-y)在一定的范围内就是一个矩形。
就成了一个20000*20000的矩阵，单点增加，询问一个子矩阵的权值和……
原来想的是树状数组，给id套一个Hash，但是HDU内存有些紧，MLE无解。
没办法就用树套树了，用树状数组维护前缀套一个Treap，Treap询问区间和
空间复杂度就是O(mlogn)了。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <list>
#define pb push_back
#define mp make_pair
#define fi first
#define se second
//#pragma comment(linker, "/STACK:16777216")
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
#define N 1200005
#define M 20005
int n , m , root[M];
struct Treap
{
int priority[N] , child[N][2];
pair<int , int> key[N] ;
int sum[N];
int Size[N] , cnt[N] , nodecnt;
Treap()
{
priority[0] = 0x7FFFFFFF;
}
void clear()
{
nodecnt = 1;
Size[0] = cnt[0] = sum[0] = 0 , key[0] = make_pair(0 , 0);
}
int size(int x)
{
return Size[root[x]];
}
void update(int x)
{
Size[x] = Size[child[x][0]] + Size[child[x][1]] + cnt[x];
sum[x] = sum[child[x][0]] + sum[child[x][1]] + key[x].second * cnt[x];
}
void rotate (int& x , int t)
{
int y = child[x][t];
child[x][t] = child[y][t ^ 1];
child[y][t ^ 1] = x;
update(x) , update(y);
x = y;
}
void insert(int& x , pair<int , int> k)
{
if (x)
{
if (key[x] == k)
++ cnt[x];
else
{
int t = key[x] < k;
insert(child[x][t] , k);
if (priority[child[x][t]] < priority[x])
rotate(x , t);
}
}
else
{
x = nodecnt ++;
key[x] = k;
sum[x] = k.second;
cnt[x] = 1;
priority[x] = rand() << 15 | rand();
child[x][0] = child[x][1] = 0;
}
update(x);
}
void query(int x , int l , int r , int& ans)
{
if(!x) return; int val = key[x].first;
if( l <= val && val <= r)
{
ans += key[x].second * cnt[x];
if(l == -1 << 30)
{
ans += sum[child[x][0]];
query(child[x][1] , l , r, ans);
}
else if(r == 1 << 30)
{
query(child[x][0] , l, r, ans);
ans += sum[child[x][1]];
}
else
{
query(child[x][0] , l , 1 << 30 , ans);
query(child[x][1] , -1 << 30, r , ans);
}
return;
}
if(r < val)
query(child[x][0] , l, r, ans);
if(val < l)
query(child[x][1] , l, r, ans);
}
};
Treap T;
void add(int x , int y , int w)
{
for (int i = x ; i < n + n ; i += i & -i)
T.insert(root[i] , make_pair(y , w));
}
int query(int x , int y , int z)
{
int res = 0 , val;
for (int i = x ; i > 0 ; i -= i & -i)
val = 0 , T.query(root[i] , y , z , val) , res += val;
return res;
}
void work()
{
int i , j , x , y , z , p , q;
memset(root , 0 , sizeof(root));
T.clear();
scanf("%d",&m);
/*for (i = 300000 ; ; -- i)
{
for (j = 2 ; j * j <= i ; ++ j)
if (i % j == 0) break;
if (j * j > i) break;
}cout <<i << endl;*/
while (m --)
{
scanf("%d%d%d%d",&i,&p,&q,&z);
if (i == 1)
{
x = p + q - 1;
y = p - q + n;
add(x , y , z);
}
else
{
x = p + q - 1;
y = p - q + n;
i = x - z , j = y - z , x += z , y += z;
i = max(1 , i) , j = max(1 , j);
x = min(n + n - 1 , x) , y = min(n + n - 1 , y);
printf("%d\n" , query(x , j , y) - query(i - 1 , j , y));
}
}
}
int main()
{
while (scanf("%d",&n) , n)
//int _ ; scanf("%d",&_); while(_--)
work();
return 0;
}
杭州剩下会做的就只有后面四道简单题了。。
金华
hdu4445 Crazy Tank
求出所有存在炮弹落在L1或L2的角度，然后挨个判定就没问题了。
但是这个角度好像不好求。。。
那就挨个枚举角度吧……
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <list>
#define pb push_back
#define mp make_pair
#define fi first
#define se second
//#pragma comment(linker, "/STACK:16777216")
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
#define N 1005
#define M 20005
int n;
double l1 , l2 , r1 , r2 , v[N] , h , pi = acos(-1.) , alpha , g = 9.8;
int cal()
{
int sum = 0 , i ;
for (i = 0 ; i < n ; ++ i)
{
double vx = sin(alpha) * v[i] , vy = cos(alpha) * v[i];
double det = vy * vy + 2 * g * h , t = (-vy + sqrt(det)) / g , x = t * vx;
if (x >= l1 && x <= r1) ++ sum;
if (x >= l2 && x <= r2) return 0;
}
return sum;
}
void work()
{
int i , ans = 0;
scanf("%lf%lf%lf%lf%lf",&h,&l1,&r1,&l2,&r2);
for (i = 0 ; i < n ; ++ i) scanf("%lf",&v[i]);
double eps = pi / 500.0;
for (alpha = 0 ; alpha - 1e-9 <= pi ; alpha += eps)
{
ans = max(ans , cal());
if (ans == n )break;
}
printf("%d\n" , ans);
}
int main()
{
while (scanf("%d",&n) , n)
//int _ ; scanf("%d",&_); while(_--)
work();
return 0;
}
hdu4447 Yuanfang, What Do You Think?
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <queue>
#include <stack>
#include <utility>
#include <map>
#define pb push_back
#define mp make_pair
#define N 1205
using namespace std;
typedef long long ll;
struct poly
{
int drg , a[N];
};
poly p[N];
poly operator / (poly x , poly y)
{
poly ans = {};
ans.drg = x.drg - y.drg;
for (int i = x.drg ;i >= y.drg ; -- i)
if (x.a[i])
{
ans.a[i - y.drg] = x.a[i] / y.a[y.drg];
for (int j = 0 ;j <= y.drg ; ++ j)
x.a[i - y.drg + j] -= ans.a[i - y.drg] * y.a[j];
}
while (ans.drg > 1 && !ans.a[ans.drg]) -- ans.drg;
return ans;
}
bool cmp(int x , int y)
{
if (p[x].drg == p[y].drg)
for (int i = p[x].drg ; i >= 0 ; -- i)
if (abs(p[x].a[i]) == abs(p[y].a[i]))
{
if (p[x].a[i] != p[y].a[i])
return (p[x].a[i] < p[y].a[i]);
}
else return abs(p[x].a[i]) < abs(p[y].a[i]);
return p[x].drg < p[y].drg;
}
void print(poly x)
{
printf("(");
for (int i = x.drg ; i >= 1 ; i --)
if (x.a[i])
{
if (x.a[i] < 0)
putchar('-');
else if (i != x.drg) putchar('+');
if (abs(x.a[i]) != 1)
printf("%d",abs(x.a[i]));
putchar('x');
if (i > 1)
printf("^%d",i);
}
if (x.a[0] > 0) printf("+%d",x.a[0]);
if (x.a[0] < 0) printf("%d",x.a[0]);
printf(")");
}
void init()
{
int i , j;
for (i = 1 ; i <= 1200; i ++)
{
p[i].drg = i , p[i].a[0] = -1 , p[i].a[i] = 1;
for (j = 1 ; j < i ;j ++)
if (i % j == 0)
p[i] = p[i] / p[j];
}
}
int n ;
void work()
{
vector<int> ans;
if (n == 1)
{
puts("x-1");
return;
}
for (int i = 1 ; i <= n ; i ++)
if (n % i == 0)
ans.pb(i);
sort(ans.begin() , ans.end() , cmp);
for (int i = 0 ;i < int(ans.size()); ++ i)
print(p[ans[i]]);
printf("\n");
}
int main()
{
init();
while (cin >> n , n) work();
return 0;
}