[ZJOI2008]骑士

[ZJOI2008]骑士

Time Limit:10000MS  Memory Limit:165536K
Total Submit:62 Accepted:30
Case Time Limit:2000MS

Description

Z国的骑士团是一个很有势力的组织,帮会中汇聚了来自各地的精英。他们劫富济贫,惩恶扬善,受到社会各界的赞扬。
最近发生了一件可怕的事情,邪恶的Y国发动了一场针对Z国的侵略战争。战火绵延五百里,在和平环境中安逸了数百年的Z国又怎能抵挡的住Y国的军 队。于是人们把所有的希望都寄托在了骑士团的身上,就像期待有一个真龙天子的降生,带领正义打败邪恶。
骑士团是肯定具有打败邪恶势力的能力的,但是骑士们互相之间往往有一些矛盾。每个骑士都有且仅有一个自己最厌恶的骑士(当然不是他自己),他是绝 对不会与自己最厌恶的人一同出征的。
战火绵延,人民生灵涂炭,组织起一个骑士军团加入战斗刻不容缓!国王交给了你一个艰巨的任务,从所有的骑士中选出一个骑士军团,使得军团内没有矛 盾的两人(不存在一个骑士与他最痛恨的人一同被选入骑士军团的情况),并且,使得这支骑士军团最具有战斗力。
为了描述战斗力,我们将骑士按照1至N编号,给每名骑士一个战斗力的估计,一个军团的战斗力为所有骑士的战斗力总和。

Input

第一行包含一个正整数N,描述骑士团的人数。
接下来N行,每行两个正整数,按顺序描述每一名骑士的战斗力和他最痛恨的骑士。

Output

应包含一行,包含一个整数,表示你所选出的骑士军团的战斗力。

Sample Input

3
10 2
20 3
30 1

Sample Output

30
【数据规模】
对于30%的测试数据,满足N ≤ 10;
对于60%的测试数据,满足N ≤ 100;
对于80%的测试数据,满足N ≤ 10 000。
对于100%的测试数据,满足N ≤ 1 000 000,每名骑士的战斗力都是不大于 1 000 000的正整数。

Source

啊啊啊。总算A掉了。。
61.187.179.132:8080/JudgeOnline/showproblem
注意到这样形成的图的每个联通块顶上必然是一堆树根组成的环,那么对每个树用树形DP算出选根和不选根的最大收入,然后就变成环上的问题了,就比较简单了。。
自己画的图,好悲剧囧。。

不过感觉我这个方法还不错,就是不是用树形DP去算,而是从底往上更新,每一种状态和规则都写一个类来表示更不容易出错(*^__^*) 。。
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1,maxn=1000000;
struct StateTree
{
long long A[2];
long long operator()(int v)const{return A[v];}
StateTree(int val=0)
{
A[0]=0;A[1]=val;
}
void Renew(const StateTree&son)
{
A[0]+=max(son(0),son(1));
A[1]+=son(0);
}
}Dp[maxn];
struct StateCircle
{
long long A[2][2];
long long operator()(int l,int r){return A[l][r];}
StateCircle(const StateTree&x)
{
memset(A,0,sizeof(A));
A[0][0]=x(0);
A[1][1]=x(1);
}
void Renew(const StateTree&x)
{
long long Ans[2][2];
rep(i,2)
{
Ans[i][0]=max(A[i][0],A[i][1])+x(0);
Ans[i][1]=A[i][0]+x(1);
}
memcpy(A,Ans,sizeof(A));
}
};
int In[maxn],P[maxn],C[maxn],n;
int main()
{
//freopen("in","r",stdin);
scanf("%d",&n);
rep(i,n) scanf("%d %d",C+i,P+i),In[–P[i]]++,Dp[i]=StateTree(C[i]);
queue<int> Q;
rep(i,n)if(!In[i]) Q.push(i);
while(Q.size())
{
int t=Q.front();Q.pop();
Dp[P[t]].Renew(Dp[t]);
if(–In[P[t]]==0) Q.push(P[t]);
}
long long Ans=0;
rep(i,n)if(In[i]>0)
{
In[i]=0;int t=i;
StateCircle now=StateCircle(Dp[t]);
t=P[t];
while(t!=i)
{
now.Renew(Dp[t]);In[t]=0;
t=P[t];
}
Ans+=max(now(1,0),max(now(0,0),now(0,1)));
}
cout<<Ans<<endl;
}

[HAOI2007]受欢迎的牛

61.187.179.132:8080/JudgeOnline/showproblem
额。。这是很老的强联通题目了。刚刚弄懂了Tarjan算法。。就像实践一下。。结果WA了
很久。。最后发现原来是Stack搞错了。。真是悲剧。。
就是求出每个强联通分量,然后如果入度为0的分量只有一个,答案就是这个分量的大小,否则就是0
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stack>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1,maxn=10000;
vector<int> E[maxn];
typedef vector<int>::iterator it;
int ord[maxn],low[maxn],cnt=0,id[maxn],snt=0,n,m;
bool inStack[maxn]={0};
stack<int> Stack;
void dfs(int x)
{
ord[x]=low[x]=cnt++;
inStack[x]=true;Stack.push(x);
for(it i=E[x].begin();i!=E[x].end();++i)
{
if(ord[*i]==-1)
dfs(*i),low[x]=min(low[x],low[*i]);
else if(inStack[*i]) low[x]=min(low[x],ord[*i]);
}
if(ord[x]==low[x])
{
int u;do{u=Stack.top();Stack.pop();id[u]=snt;inStack[u]=false;}while(u!=x);
snt++;
}
}
int Num[maxn]={0},In[maxn]={0};
int main()
{
//freopen("in","r",stdin);
cin>>n>>m;int s,t;
while(m–) scanf("%d %d",&s,&t),E[t-1].pb(s-1);
rep(i,n) ord[i]=-1;
rep(i,n)if(ord[i]==-1) dfs(i);
rep(i,n){
int own=id[i];for(it e=E[i].begin();e!=E[i].end();++e)if(id[*e]!=own) In[id[*e]]++;
Num[own]++;
}
t=0;int ans=0;rep(i,snt)if(In[i]==0) ++t,ans+=Num[i];
if(t>1)cout<<0<<endl;
else cout<<ans<<endl;
}

[ZJOI2008]生日聚会Party

[ZJOI2008]生日聚会Party

Time Limit:10000MS  Memory Limit:165536K
Total Submit:35 Accepted:20
Case Time Limit:1000MS

Description

今天是hidadz小朋友的生日,她邀请了许多朋友来参加她的生日party。 hidadz带着朋友们来到花园中,打算坐成一排玩游戏。为了游戏不至于无聊,就座的方案应满足如下条件:
对于任意连续的一段,男孩与女孩的数目之差不超过k。
很快,小朋友便找到了一种方案坐了下来开始游戏。hidadz的好朋友Susie发现,这样的就座方案其实是很多的,所以大家很快就找到了一种, 那么到底有多少种呢?热爱数学的hidadz和她的朋友们开始思考这个问题……
假设参加party的人中共有n个男孩与m个女孩,你是否能解答Susie和hidadz的疑问呢?由于这个数目可能很多,他们只想知道这个数目 除以12345678的余数。

Input

仅包含一行共3个整数,分别为男孩数目n, 女孩数目m, 常数k。

Output

应包含一行,为题中要求的答案。

Sample Input

1 2 1

Sample Output

1

Hint

对于30%的数据,n , m ≤ 20;
对于100%的数据, n , m ≤ 150,k ≤ 20。

Source

61.187.179.132:8080/JudgeOnline/showproblem
额。就是DP不知道有没有什么公式的。
设dp(i,j,k1,k2)是有i个人,j个男的,右端最多男生比女生多k1个,最多女生比男生多k2个,
跟算最大连续子序列的线段树差不多。。
然后就可以直接DP了。。由于空间不够,要用滚动数组额。
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1,maxn=150+1,maxk=20+1,mod=12345678;
int n,m,k;
int dp[2][maxn][maxk][maxk];
void add(int&a,int x)
{
a+=x;a%=mod;
}
int main()
{
//freopen("in","r",stdin);
cin>>n>>m>>k;
int now=0,old=1;
dp[now][0][0][0]=1;
for(int i=0;i<n+m;++i)
{
swap(now,old);
memset(dp[now],0,sizeof(dp[now]));
for(int j=0;j<=min(i,n);++j)
for(int k1=0;k1<=min(j,k);++k1)
for(int k2=0;k2<=min(i-j,k);++k2)
{
int x=dp[old][j][k1][k2];
if(k1<k&&j<n) add(dp[now][j+1][k1+1][max(k2-1,0)],x);
if(k2<k&&i-j<m) add(dp[now][j][max(k1-1,0)][k2+1],x);
}
}
int ans=0;
rep(k1,k+1)rep(k2,k+1) ans+=dp[now][n][k1][k2],ans%=mod;
cout<<ans<<endl;
}

[HNOI2008]GT考试

61.187.179.132:8080/JudgeOnline/showproblem
额。。水题直接上矩阵乘法。。
由于m比较小实际上没必要用KMP算法算出原始矩阵的,不过我乘机学习了一下KMP算法囧。。
还有就是我昨天一直在调我那个QTREE的程序。。到现在都没有调好。。我真是太菜了555.
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1,maxm=20;
int n,m,mod,A[maxm],next[maxm+1];
struct Mat
{
int M[maxm][maxm];
Mat(){memset(M,0,sizeof(M));}
void operator=(const Mat&o)
{memcpy(M,o.M,sizeof(M));}
int& operator()(int i,int j){return M[i][j];}
Mat operator*(Mat&o)
{
Mat T;
rep(i,m)rep(j,m)rep(k,m)T(i,j)+=M[i][k]*o(k,j),T(i,j)%=mod;
return T;
}
}orig;
Mat Power(int n)
{
if(n==1) return orig;
Mat tmp=Power(n/2);
tmp=tmp*tmp;
if(n&1) tmp=tmp*orig;
return tmp;
}
int main()
{
//freopen("in","r",stdin);
cin>>n>>m>>mod;char x;
rep(i,m) cin>>x,A[i+1]=x-‘0’;
next[1]=0;int t=0;
for(int i=2;i<=m;i++)
{
while(t>0&&A[t+1]!=A[i]) t=next[t];
if(A[t+1]==A[i])++t;
next[i]=t;
}
for(int i=0;i<m;++i)
{
rep(j,10)
{
int t=i;
while(t>0&&A[t+1]!=j)t=next[t];
if(A[t+1]==j)++t;
orig(t,i)++;
}
}
Mat ans=Power(n);int Ans=0;
rep(i,m) Ans+=ans(i,0),Ans%=mod;
cout<<Ans<<endl;
}

[HNOI2008]明明的烦恼

61.187.179.132:8080/JudgeOnline/showproblem
这道题很有意思啊。。就是说一棵树对一些顶点的度数有要求一些没有,让你求这样的树的度数。
我一开始想DP脑子都爆了也没想出来。。最后想到Prufer Code。。然后就明白了。。
Prufer Code是把每棵n点的树对应到长度为n-2,数在1-n之间的数列。。其中一个度数为d的点出现d-1次。。
wiki上的说明:en.wikipedia.org/wiki/Prufer_code
这就是关键。所以问题就转化成求这样的数列的数目了。。这用多重集的公式很方便,就不多说了。。
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1,maxn=1000;
int n;
vector<int> Prime;
typedef vector<int>::iterator it;
int D[maxn],S=0,noLimit=0;
bool isPrime(int x)
{
if(x==2) return true;
if(x/2==0) return false;
for(int i=3;i*i<=x;i+=2)
if(x%i==0) return false;
return true;
}
void getPrime()
{
for(int i=2;i<=n;i++)
if(isPrime(i)) Prime.pb(i);
}
void init()
{
cin>>n;
rep(i,n)
{
cin>>D[i];
if(D[i]>0)
S+=D[i]-1;
else
noLimit++;
}
}
struct BigInt
{
int H[400];
BigInt(){memset(H,0,sizeof(H));}
void mul(int x)
{
rep(i,Prime.size())
while(x%Prime[i]==0) x/=Prime[i],H[i]++;
}
void div(int x)
{
rep(i,Prime.size())
while(x%Prime[i]==0) x/=Prime[i],H[i]–;
}
};
ostream& operator<<(ostream&out,const BigInt&x)
{
static int Hp[100000],last;
memset(Hp,0,sizeof(Hp));Hp[last=0]=1;
rep(i,Prime.size())
{
rep(j,x.H[i])
{
int o=Prime[i],d=0;
rep(k,last+1)
{
d+=Hp[k]*o;
Hp[k]=d%10;
d/=10;
}
while(d) Hp[++last]=d%10,d/=10;
}
}
for(int i=last;i>=0;–i)out<<Hp[i];
return out;
}
void work()
{
if(S>n-2){cout<<0<<endl;return;}
if(noLimit==0&&S<n-2){cout<<0<<endl;return;}
getPrime();
BigInt ans;
for(int i=1;i<=n-2;i++)
ans.mul(i);
rep(i,n)
{
if(D[i]>0)
for(int j=1;j<=D[i]-1;j++)
ans.div(j);
}
for(int i=1;i<=n-2-S;i++)
ans.div(i);
rep(i,n-2-S) ans.mul(noLimit);
cout<<ans<<endl;
}
int main()
{
//freopen("in","r",stdin);
init();
work();
}
PS。。o(︶︿︶)o 这个OJ不能用Java。还得自己写高精度囧。。

TCH 2010-2

额。。这两天我的胃罢工了。。难过4了。。今天抱着别被淘汰就行的心理参加了这个比赛囧。。
比赛前狂灌咖啡。。感觉好多了。。
第一题是水题啊。是水题啊。
第二题真BT啊,真BT。。
一开始我很想啊,这个要分类讨论啊。。然后我脑海里就浮现出200行+的代码和密密麻麻的公式。。囧了一下。。然后想是不是把这个东西摊开来啊。然后写了一会儿。囧掉了。写不出来。。最后没办法了。。突然想到似乎可以bfs。。然后又想如何表示6个面。。纠结了半天。。只好又灌咖啡。。最后豁出去了。干脆开3维数组得了。。结果居然过了。。其实我那个程序好像有个bug总觉得哪个地方搞错了。。不过居然Pass System Test了。真是奇迹额。。
第三题去死吧。去死吧。。
我想到上次USACO的月赛。然后感觉要离散化,然后离散化的时候把角要切出来,之后开个数组递推出每个其他的方格是否是白的?现在感觉好像是这么做的,不过比赛的时候没想到。我的第一反应是绕圈儿。。
因为这个图形本质上好像就是一圈东西囧。。
然后好不容易写了个绕圈的程序。。发现不知道怎么算面积囧。。然后我就撞死了。。
掉Rating了。。

COCI 2010-5

今天两个比赛之间一秒钟都没有。。一连码了5个小时。。脑子小了囧。。
题目见官网。。
www.hsin.hr/coci/
第一题第二题纯水题。。忽略它们吧。。
第三题就是去重么,注意旋转就可以了。。
第四题我写了个dfs。好像写错了。。只有10分囧。。
第五题其实比第四题简单,对于每种jump一起上,那么每个数最多1次,总共也就是nlogn的。。
第六题太BT了。。我看了就晕。。
额。。这次比上次好点。。280分。。
我这个人一做到烦的题目就晕。。

[HNOI2006]花仙子的魔法

囧。。其实这玩意等价于分空间的问题,设维数为d魔法n次为dp(d,n)
则dp(d,n)=dp(d-1,n-1)+dp(d,n-1)
dp(1,i)=2*i
dp(i,0)=1
。。。我一开始是解了个方程。因为我发现结果是d次函数,并且前d+1个值都是2的幂,
所以可以解出系数然后计算囧。。写了半天。真是沙茶到家了555
HNOI2006A了5道。。剩下的就不会做了。。
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1;
const int maxn=100+1,maxd=16;
int n,d;
long long dp[maxd][maxn]={0};
int main()
{
cin>>n>>d;
rep(i,n+1)dp[1][i]=i*2;dp[1][0]=1;
for(int i=2;i<=d;i++)
{
dp[i][0]=1;
for(int j=1;j<=n;j++)
dp[i][j]=dp[i][j-1]+dp[i-1][j-1];
}
cout<<dp[d][n]<<endl;
}

[HNOI2006]公路修建问题

61.187.179.132:8080/JudgeOnline/showproblem
看似水题。。其实不是很水(应该是我太菜了囧。。)
之前想了几种算法都是错的。。就不说了。最后的算法是二分最大值。并用并查集帮忙(*^__^*)
然后首先扫描一遍所有的边,能够建一号公路的就建一号公路。扫描的话如果k没达到,就是false
然后再扫描一遍所有的边,把能建的二号公路也建上去。判断一下是否联通的就OK了。。
Code:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstring>
#include<queue>
#define rep(i,n) for(int i=0;i<n;i++)
#define pb push_back
using namespace std;
const int inf=~0U>>1;
const int maxn=10000;
int n,k,m;
struct Data
{
int s,t,c1,c2;
void operator()(int _s,int _t,int _c1,int _c2)
{s=_s;t=_t;c1=_c1;c2=_c2;}
}E[maxn*2];
int F[maxn];
int find(int x)
{
if(F[x]==x) return x;
return F[x]=find(F[x]);
}
int Union(int i,int j){F[j]=i;}
bool check(int limit)
{
int cnt=n,r=k;
rep(i,n) F[i]=i;
for(Data*i=E;i!=E+m;i++)
{
if(i->c2<=limit)
{
int a=find(i->s),b=find(i->t);
if(a!=b) Union(a,b),–cnt,–r;
}
}
if(r>0) return false;
for(Data*i=E;i!=E+m;i++)
{
if(i->c1<=limit)
{
int a=find(i->s),b=find(i->t);
if(a!=b) Union(a,b),–cnt;
}
if(cnt==1) return true;
}
return false;
}
int main()
{
//freopen("in","r",stdin);
cin>>n>>k>>m;int s,t,c1,c2;
rep(i,m)
{
scanf("%d %d %d %d",&s,&t,&c2,&c1);
E[i](s-1,t-1,c1,c2);
}
int l=0,r=30000;
while(l+1<r)
{
int m=(l+r)/2;
if(check(m))
r=m;
else
l=m;
}
cout<<r<<endl;
}