8.22 NOIP 模拟题
编译命令 g++ -o * *.cpp gcc -o * *.c fpc *.pas 编译器版本 g++/gcc 4.9.2 fpc 2.6.2 评测环境 64 位 Linux, 3.3GHZ CPU 评测软件 Lemon 评测方式 忽略行末空格和回车 特别注意：c/c++ 选手使用 printf 输出 64 位整数请使用%lld 1
A 债务
文件名 输入文件 输出文件 时空限制
debt.pas/c/cpp debt.in debt.out 1s 128MB
题目描述
小 G 有一群好朋友，他们经常互相借钱。假如说有三个好朋友 A，B，C。A
欠 B 20 元，B 欠 C 20 元，总债务规模为 20+20=40 元。小 G 是个追求简约的人，
他觉得这样的债务太繁杂了。他认为，上面的债务可以完全等价为 A 欠 C 20 元，
B 既不欠别人，别人也不欠他。这样总债务规模就压缩到了 20 元。
现在给定 n 个人和 m 条债务关系。小 G 想找到一种新的债务方案，使得每个
人欠钱的总数不变，或被欠钱的总数不变（但是对象可以发生变化） ，并且使得总
债务规模最小。
输入格式
输入文件第一行两个数字 n,m，含义如题目所述。
接下来 m 行，每行三个数字 a i ,b i ,c i ，表示 a i 欠 b i 的钱数为 c i 。
注意，数据中关于某两个人 A 和 B 的债务信息可能出现多次，将其累加即可。
如”A 欠 B 20 元”、”A 欠 B 30 元”、”B 欠 A 10 元”，其等价为”A 欠 B 40 元”。
输出格式
输出文件共一行，输出最小的总债务规模。
样例输入 1
5 3
1 2 10
2 3 1
2 4 1
2
样例输出 1
10
样例输入 2
4 3
1 2 1
2 3 1
3 1 1
样例输出 2
0
数据范围
对于 30% 的数据，1 ≤ n ≤ 10，1 ≤ m ≤ 10。
对于 60% 的数据，1 ≤ n ≤ 100, 1 ≤ m ≤ 10 4 。
对于 80% 的数据，1 ≤ n ≤ 10 4 ，1 ≤ m ≤ 10 4 。
对于 100% 的数据，1 ≤ n ≤ 10 6 ，1 ≤ m ≤ 10 6 。
对于所有的数据，保证 1 ≤ a i ,b i ≤ n,0 < c i ≤ 100。
#include<iostream> #include<cstdio> #include<cstring> #define N 1000007 using namespace std; int val[N]; int x,y,z,n,m,ans; inline int read() { int x=0,f=1;char c=getchar(); while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();} while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();} return x*f; } int main() { freopen("debt.in","r",stdin); freopen("debt.out","w",stdout); n=read(),m=read(); for(int i=1;i<=m;i++) { x=read();y=read();z=read(); val[y]+=z;val[x]-=z; } for(int i=1;i<=n;i++) if(val[i]<0) ans+=-val[i]; printf("%d ",ans); fclose(stdin);fclose(stdout); return 0; }
B 排列
文件名 输入文件 输出文件 时空限制
perm.pas/c/cpp perm.in perm.out 1s 128MB
题目描述
小 G 喜欢玩排列。现在他手头有两个 n 的排列。n 的排列是由 0,1,2,...,n−1
这 n 的数字组成的。对于一个排列 p，Order(p) 表示 p 是字典序第 Order(p) 小的
排列（从 0 开始计数） 。对于小于 n! 的非负数 x，Perm(x) 表示字典序第 x 小的
排列。
现在，小 G 想求一下他手头两个排列的和。两个排列 p 和 q 的和为 sum =
Perm((Order(p) + Order(q))%n!)。
输入格式
输入文件第一行一个数字 n，含义如题。
接下来两行，每行 n 个用空格隔开的数字，表示小 G 手头的两个排列。
输出格式
输出一行 n 个数字，用空格隔开，表示两个排列的和。
样例输入 1
2
0 1
1 0
样例输出 1
1 0
4
样例输入 2
3
1 2 0
2 1 0
样例输出 2
1 0 2
数据范围
1、2、3、4 测试点，1 ≤ n ≤ 10。
5、6、7 测试点，1 ≤ n ≤ 5000，保证第二个排列的 Order ≤ 10 5 。
8、9、10 测试点，1 ≤ n ≤ 5000。
/* 可以看出Order perm 互为反函数 康托展开： 一个序列的排名(从0开始计数) = Rank[n]*(n-1)!+Rank[n-1]*(n-2)!+…. 其中Rank[n]表示n位置上的数字在a[1]~a[n]中的排行，并且从0开始计数。 先求出两个序列的康托展开式，相加 但是显然阶乘爆掉 于是我们只加rank数组 如果rank[i]这一位大于等于i，就按i进制进位 证明： ……rank[i] * (i - 1) ! + rank[i + 1] * i !…… 若rank[i] > i 那么可分解为 rank[i]%i * (i - 1)! + rank[i]/i * i! + rank[i + 1]*i! 所以进位显然 mod n!的话，只需要忽略rank[i + 1]即可 可以从n往前找，p从-1向上累计，遇到没有用过的p就标记为用过，同时rank[n]-- 因为num[i + 1..n]求出后能得知num[i]能选哪一些数 也就是num[1...i]有哪一些数 从小到大枚举到对应排名即可 */ #include<iostream> #include<cstdio> #include<cstring> #define N 6010 #define inf 0x3f3f3f3f using namespace std; int num1[N],num2[N],rank[N],rank2[N]; int n,b[N]; int main() { freopen("perm.in","r",stdin); freopen("perm.out","w",stdout); scanf("%d",&n); for(int i=n;i>0;i--) scanf("%d",&num1[i]); for(int i=n;i>0;i--) scanf("%d",&num2[i]); int p; for(int i=1;i<=n;i++) { p=0; for(int j=1;j<i;j++) if(num1[j]<num1[i]) ++p; rank[i]+=p; p=0; for(int j=1;j<i;j++) if(num2[j]<num2[i]) ++p; rank2[i]+=p; rank[i]+=rank2[i]; rank[i+1]+=rank[i]/i; rank[i]%=i; } rank[n]%=n; for(int i=n;i>=1;i--) { p=-1; while(rank[i]>=0) { ++p; if(!b[p]) --rank[i]; }b[p]=1; printf("%d ",p); }return 0; }
C 剪树枝
文件名 输入文件 输出文件 时空限制
tree.pas/c/cpp tree.in tree.out 1s 128MB
题目描述
rzyz 有一棵苹果树。苹果树有 n 个节点（也就是苹果） ，n − 1 条边（也就是
树枝） 。调皮的小 G 爬到苹果树上。他发现这棵苹果树上的苹果有两种：一种是黑
苹果，一种是红苹果。小 G 想要剪掉 k 条树枝，将整棵树分成 k + 1 个部分。他
想要保证每个部分里面有且仅有一个黑苹果。请问他一共有多少种剪树枝的方案？
输入格式
第一行一个数字 n，表示苹果树的节点（苹果）个数。
第二行一共 n − 1 个数字 p 0 ,p 1 ,p 2 ,p 3 ,...,p n−2 ，p i 表示第 i + 1 个节点和 p i 节
点之间有一条边。注意，点的编号是 0 到 n − 1。
第三行一共 n 个数字 x 0 ,x 1 ,x 2 ,x 3 ,...,x n−1 。如果 x i 是 1，表示 i 号节点是黑
苹果；如果 x i 是 0，表示 i 号节点是红苹果。
输出格式
输出一个数字，表示总方案数。 答案对 10 9 + 7 取模。
样例输入 1
3
0 0
0 1 1
样例输出 1
2
6
样例输入 2
6
0 1 1 0 4
1 1 0 0 1 0
样例输出 2
1
样例输入 3
10
0 1 2 1 4 4 4 0 8
0 0 0 1 0 1 1 0 0 1
样例输出 3
27
数据范围
对于 30% 的数据，1 ≤ n ≤ 10。
对于 60% 的数据，1 ≤ n ≤ 100。
对于 80% 的数据，1 ≤ n ≤ 1000。
对于 100% 的数据，1 ≤ n ≤ 10 5 。
对于所有数据点，都有 0 ≤ p i ≤ n − 1，x i = 0 或 x i = 1。
特别地，60% 中、80% 中、100% 中各有一个点，树的形态是一条链。
/* 一道比价杂乱的树形dp，但貌似很经典。 f[i][0][0]代表第i个点是否砍掉，以它为根的子树中有没有黑点。 转移的时候判断当前点是黑是红 黑点：f[u][1][0]=0;f[u][0][0]=0;（这两种情况不成立，不管砍不砍以它为根的子树中一定有黑点） f[u][1][1]=f[v][1][1]+f[v][0][0]的乘积（它被砍了子树中只能有它一个黑点， 所以儿子中要么没有黑点[0][0]要么有黑点被砍了[1][1]） f[u][0][1]不知道暂时赋值为f[u][1][1] 白点：f[u][1][0]这种情况对答案没有贡献设为0，因为它为白点子树仍为白点都不能割； f[u][0][0]=f[v][1][1]+f[v][0][0]的乘积 （它的子树中没有黑点，儿子中要么没有黑点[0][0]要么有黑点被砍了[1][1]） f[u][1][1]=∑f[x][0][1]*（f[v][1][1]+f[v][0][0]），它为白以它为根的子树中有且只有一个黑点， 则只有一个黑点儿子不被砍，其余儿子中要么没有黑点[0][0]要么有黑点被砍了[1][1] 最后答案为f[0][0][1] */ #include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #define mod 1000000007 #define N 100005 using namespace std; long long f[N][2][2]; int tot,n; int head[N],col[N],son[N]; struct node { int to,next; }e[N<<2]; inline int read() { int x=0,f=1;char c=getchar(); while(c>'9'||c<'0'){if(c=='=')f=-1;c=getchar();} while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();} return x*f; } inline void add(int u,int v) { e[++tot].to=v;e[tot].next=head[u];head[u]=tot; } void dfs(int u,int v) { for(int i=head[u];~i;i=e[i].next) { if(e[i].to==v) continue; dfs(e[i].to,u); son[u]=1; } } void dp(int u,int v) { if(!son[u]) { if(col[u]) { f[u][1][1]=1;f[u][1][0]=0; f[u][0][1]=1;f[u][0][0]=0; } else { f[u][1][1]=0;f[u][1][0]=1; f[u][0][1]=0;f[u][0][0]=1; }return; } for(int i=head[u];~i;i=e[i].next) { if(e[i].to==v) continue; dp(e[i].to,u); } if(col[u]) { f[u][1][0]=0;f[u][0][0]=0; long long sum=1; for(int i=head[u];~i;i=e[i].next) { int t=e[i].to; if(t==v)continue; sum=sum*(f[t][0][0]+f[t][1][1])%mod; } f[u][0][1]+=sum;f[u][1][1]+=sum; } else { f[u][1][0]=0; long long sum=1; for(int i=head[u];~i;i=e[i].next) { int t=e[i].to; if(t==v)continue; sum=sum*(f[t][0][0]+f[t][1][1])%mod; } f[u][0][0]+=sum; long long ans=0;sum=1; for(int i=head[u];~i;i=e[i].next,sum=1) { if(e[i].to==v)continue; int tmp=f[e[i].to][0][1]; if(!tmp)continue; for(int j=head[u];~j;j=e[j].next) { int t=e[j].to; if(t==v || j==i) continue; sum=sum*(f[t][0][0]+f[t][1][1])%mod; } ans=(ans+tmp*sum)%mod; } f[u][1][1]=ans%mod;f[u][0][1]=ans%mod; } } int main() { // freopen("tree.in","r",stdin); // freopen("tree.out","w",stdout); /* 手动扩栈 int size = 256 << 16; // 256MB char *p = (char*)malloc(size) + size; __asm__("movl %0, %%esp " :: "r"(p)); */ memset(head,-1,sizeof head); n=read();int x; for(int i=0;i<n-1;i++) { x=read(); add(x,i+1);add(i+1,x); } for(int i=0;i<n;i++) col[i]=read(); dfs(0,-1);dp(0,-1); cout<<f[0][0][1]%mod; return 0; }
100+40+0 T2没想出来。 嗯没错，T3暴力挂了，水了30分链一分没得，这让大数据结构如此擅长水链分的我特别桑心了，，，，，， 唉，愁得我......T3还是没耐心写暴力吧(想到80分暴力貌似不对，写了也不一定调出来)，还是太弱了，myj大佬都开始挨个刷bzoj了.....不怂就是干！晚上吧全排列写5边！！！(被stl坑了一把)。 T3正解一晚上没懂，上午半上午才又找了一个看懂了，唉...... 我就不信这个邪......