2015年09月_yysys

原创 BZOJ 2440: [中山市选2011]完全平方数

二分,然后用莫比乌斯做就行了... #include using namespace std; typedef long long LL; const int maxn = 50005; int not_prime[maxn]; int prime[maxn]; int mu[maxn]; int K, p_cnt; void Init() { mu[1] = 1; p_cnt =

2015-09-30 15:06:24 611

原创 DP HDOJ 5492 Find a path

化简公式然后递推.... #include using namespace std; typedef long long LL; const int INF = 0x3f3f3f3f; int dp[33][33][60 * 33]; int g[33][33], n, m; void update(int &x, int y) { if(x > y) x = y; } void w

2015-09-29 20:00:14 416

原创 HDOJ 5483 Nux Walpurgis

先做出最小生成树,然后对每个点dfs判断树边是否可删..... #include using namespace std; typedef long long LL; const int maxn = 3005; const int maxm = 6005; const int INF = 0x3f3f3f3f; struct Edge { int v; Edge *next; }E[

2015-09-29 14:51:43 426

原创 DP Codeforces Round #322 (Div. 2) F. Zublicanes and Mumocrates

树形dp,记录一下根节点的状态即可... #include using namespace std; typedef long long LL; const int maxn = 5005; const int maxm = 10005; const int INF = 0x3f3f3f3f; struct Edge { int v; Edge *next; }E[maxm], *H[

2015-09-28 21:40:39 643

原创 HDOJ 5487 Difference of Languages

记dp[i][j] 代表走到了第一个dfa的i号节点,走到了第二个自动机的j号节点.....然后就可以转移了...... #include using namespace std; typedef long long LL; const int maxn = 1005; struct node { int s1, s2, pre, tran; node(int s1 = 0, int

2015-09-28 15:17:23 448

原创 HDOJ 5469 Antonidas 树分治

经典的树分治,找字符串用hash找就行了.... #include using namespace std; typedef long long LL; #define mp(x, y) make_pair(x, y) const int maxn = 10005; const int maxm = 20005; const int INF = 0x3f3f3f3f; const int m

2015-09-28 13:12:10 428

原创［DP］ HDOJ 5456 Matches Puzzle Game

数位ＤＰ．．． #include using namespace std; typedef long long LL; LL dp[505][2][2][2]; int cnt[] = {6, 2, 5, 5, 4, 5, 6, 3, 7, 6}; int mod; void add(LL& a, LL b) { a += b; if(a > mod) a -= mod; } LL

2015-09-25 13:28:45 394

原创［半平面交］ HDOJ 5462 Manors

很容易就能看出半平面交的公式．．．但是比赛的时候这题看都来不及看．．． #include using namespace std; typedef long long LL; const int maxn = 1005; const double eps = 1e-7; struct Point { double x, y; Point(double x = 0, double y =

2015-09-23 20:00:30 524

原创［网络流］ HDOJ 5457 Hold Your Hand

建出字典树，然后做最小割就行了．．． #include using namespace std; typedef long long LL; const int maxn = 5005; const int maxm = 400005; const int INF = 0x3f3f3f3f; struct Edge { int v, c, next; Edge() {} Edge(i

2015-09-23 18:20:45 511

原创［线段树］ HDOJ Excited Database

把一个矩阵化成３个三角形容斥，然后用等差线段树就可以做了．．． #include using namespace std; typedef long long LL; #define now o, L, R, tree #define lson o << 1, L, mid, tree #define rson o << 1 | 1, mid+1, R, tree #define ls o

2015-09-21 20:13:26 428

原创 hihocoder　#1236 : Scores

高维bitset．．．． #include using namespace std; typedef long long LL; const int maxn = 50005; struct node { int val, id; node(int val = 0, int id = 0) : val(val), id(id) {} bool operator < (const no

2015-09-21 14:28:32 735

原创［动态树］ HDOJ 5467 Clarke and hunger games

按时间点建树，然后递归实现ＬＣＴ的可持久化．．．． #include using namespace std; typedef long long LL; #define mp(x, y) make_pair(x, y) const int maxn = 300005; const int maxm = 300005; const int INF = 0x3f3f3f3f; const in

2015-09-21 10:07:57 540

原创 HDOJ 5458 Stability

由于没有看到图在任意时刻都是连通的．．这一条件．．于是只好打了动态树．． #include using namespace std; typedef long long LL; const int maxn = 30005; const int maxm = 400005; struct Node *null; struct Node { Node *fa, *ch[2]; int s

2015-09-19 18:44:40 429

原创【NTT】 ZOJ 3874 Permutation Graph

推出公式以后分治ntt优化即可。。。 #include using namespace std; typedef long long LL; const int maxn = 100005; const int maxm = 300005; const int mod = 786433; LL a[maxm], b[maxm], c[maxm]; LL f[maxn], dp[maxn

2015-09-18 13:06:22 593

原创【NTT】 HDOJ 5279 YJC plays Minecraft

推出公式以后就可以ntt了。。。 #include using namespace std; typedef long long LL; const int mod = 998244353; const int maxn = 700005; LL dp[100005],dp2[100005]; LL a[maxn], b[maxn], c[maxn], xp[maxn]; LL f[100

2015-09-17 18:21:27 428

原创【矩阵快速幂】 HDOJ 5434 Peace small elephant

对每一列状态压缩，然后矩阵加速即可。。。 #include using namespace std; typedef long long LL; const int mod = 1e9+7; const int maxn = 130; LL mat[maxn][maxn]; LL mid[maxn][maxn]; LL res[maxn][maxn]; int n, m; void ca

2015-09-16 14:04:13 354

原创【树分治】 ZOJ Travel

离线每个询问，然后做树分治。。 #include using namespace std; typedef long long LL; #define lowbit(x) (x&(-x)) #define pii pair #define mp(x, y) make_pair(x, y) const int maxn = 100005; const int maxm = 200005; c

2015-09-16 11:08:54 317

原创 HDOJ 5448 Marisa’s Cake

计算每一条线段对答案的贡献，然后用前缀和优化就行了。。。。 #include using namespace std; typedef long long LL; const int maxn = 100005; const int mod = 1e9+7; struct Point { LL x, y; Point(LL x = 0, LL y = 0) : x(x), y(y) {

2015-09-14 15:53:01 402

原创【NTT】 ZOJ 3899 State Reversing

先找出第二类斯特林数的公式，然后把公式分解成卷积的形式，先做一遍NTT，然后对于每次询问只要用线段树求出有多少个空闲的房间就可以了。 #include using namespace std; typedef long long LL; #define lson o << 1, L, mid #define rson o << 1 | 1, mid+1, R #define ls o <<

2015-09-14 13:11:49 756

原创【DP】 HDOJ 4804 Campus Design

轮廓线dp一下就行了。。和大白书上的轮廓线dp类似。。。 #include using namespace std; typedef long long LL; const int maxn = 105; const int mod = 1e9+7; char g[maxn][11]; int dp[2][1 << 10][22]; int n, m, C, D; void add(in

2015-09-11 19:59:04 307

原创【后缀自动机】 HDOJ 5431 AB String

在sam上做dp统计有多少的串不是a的子串，然后找到长度，dfs找串就行了。。 #include using namespace std; typedef long long LL; const int alpha = 2; const int maxn = 40005; const int maxm = 25; struct node { int len; LL cnt[maxm];

2015-09-06 16:20:28 379

原创【polya】 HDOJ Count the Tetris

套用polya计算出奇偶，旋转，反射16钟情况的方案数。。 import java.util.*; import java.io.*; import java.math.*; import java.lang.*; public class Main { public static void main(String args[]) { Scanner cin = new Scanne

2015-09-03 18:33:38 364

原创【DP】HDOJ 5426 Rikka with Game

把子树合并到根上考虑三种转移。1)切断根与子树的边，建立外交。2)切断根与子树的边，不建立外交。3)不切断根与子树的边。 #include using namespace std; typedef long long LL; const int maxn = 100005; const int maxm = 505; const LL INF = 4e18; struct Edge {

2015-09-02 10:16:07 477

原创 HDOJ 5425 Rikka with Tree II

和上次的上海邀请赛的那个dp优化技巧类似。都是考虑double精度，来减少枚举次数。。。 #include using namespace std; typedef long long LL; const int maxn = 100005; const int maxm = 200005; struct Edge { int v; Edge *next; }*H[maxn], E[m

2015-09-01 14:14:08 320

yysys